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1.  Introduction 


When  a  specimen  with  an  edge  crack  is  subjected  to  fatigue  loading,  amsotropic  and 
inhomogeneous  damage  will  slowly  develop  and  eventually  evolve  into  a  typical 
configuration  in  which  a  crack  of  length  unity  is  lodged  in  a  vanishingly  small  crack-tip 
damage.  Similarly  when  a  notched  specimen  is  fatigue  loaded,  a  crack  completely 
surrounded  by  damage  will  eventually  emerge.  These  two  cases  are  depicted  in  Figs.  1  and 
2  where  t<<!  and  r  (0)  defines  the  shape  of  the  damage  zone. 
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In  order  to  understand  the  many  fundamental  physical  phenomena  associated  with 
the  indicated  crack— damage  evolution,  a  meaningful  description  of  the  so-called  damage 
must  first  be  given.  Then  the  kinetic  relation  associated  with  the  damage  evolution  must 
also  be  fully  explored  and  developed.  Both  issues  have  attracted  the  attention  of  many 
researchers  for  many  years,  but  the  results  have  been  rather  inconclusive.  It  is  partly  for 
this  reason  that  neither  issue  was  considered  in  a  direct  way  in  this  project.  Instead  we 
took  an  extremely  simplistic  approach  by  assuming  that  the  main  characteristics  of  a 
damaged  material  is  that  it  is  softer  than  the  original  material.  This  assumption  reduces 
the  unknown  damage  to  an  inhomogeneity  that  may  be  characterized  by  two  elastic 
constants  which  may  in  turn  be  interpreted  as  damage  parameters. 

When  a  crack  is  interacting  with  an  inhomogeneity,  such  as  the  cases  depicted  in 
Figs.  1  and  2,  the  associated  elasticity  problem  is  well  defined  and  hence  a  full  stress 
analysis  is  possible.  The  availability  of  this  solution  makes  it  feasible  to  link  the  Eshelby 
(configuration)  forces  to  the  inhomogeneitv  geometry  and  moduli,  which  are  used  as  a  gross 
representation  for  the  damage.  Thus,  devising  an  efficient  numerical  procedure  for 
obtaining  the  desired  elasticity  solution  became  a  main  objective  of  this  one— year  project. 
This  objective  was  completed  and  the  relevant  results  may  be  found  in  Chao-Hsun  Chen’s 
thesis  which  was  completed  in  1990  under  the  supervision  of  the  PI  (Appendix  I). 

It  is  perhaps  clear  from  the  general  nature  of  the  problem  that  the  desired  solution 
can  only  be  obtained  by  numerical  means  However,  the  actual  computation  must  be 
preceeded  by  a  skillful  asymptotic  analysis  so  that  the  t-»0  limit  is  analytically  factored 
out  of  the  ensuing  calculation.  This  is  necessary  because  no  numerical  scheme  can  possibly 

—  WO 

handle  the  conflicting  limits  required  by  t->0  (almost  no  inhomogeneitv)  and  r  '  “  -*  a> 
(crack -t.ip  inside  the  inhomogeneity).  Such  an  asymptotic  analysis,  together  with  the 
accompanying  computational  scheme,  was  successfully  accomplished. 

A  more  detailed  description  of  this  accomplishment  is  summarized  in  Section  2. 
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The  above  objective  bypasses  the  need  and  difficulty  of  addressing  the  many 
physical  phenomena  associated  with  a  typical  state  of  damage,  and  yet  the  defining  of  a 
state  of  damage,  together  with  the  associated  kinetics,  remains  to  be  the  most  crucial  issue. 
It  is  with  this  consideration  in  mind  that  we  chose  to  re— examine  the  implications  of  higher 
order  theories  and  found  that  the  v/ell-established  grade— 3  theory  is  the  most  outstanding 
candidate.  This  accomplishment  is  summarized  in  Section  3. 


2.  Numerical  Solution  of  a  Crack  Interacting  with  an  Inhomogeneity 

The  elasticity  problems  considered  are  those  depicted  in  Figs.  1  and  2  where  the 
damaged  material  is  simplistically  replaced  by  a  different  elastic  material,  and  is 
henceforth  termed  an  inhomogeneity.  The  geometry  of  a  typical  problem  is  therefore 
defined  by  a  dimensionless  half-crack  length  1,  an  inhomogeneity  denoted  by  D1  which  in 
term  is  embedded  in  a  medium  denoted  by  D^.  The  associated  elasticity  problems  are 
formulated  in  terms  of  complex  functions. 

Let  (zpZ2)  be  rectangular  Cartesian  coordinates  and  z=z^  +  iz2  the  associated 
complex  variable  in  the  z-plane.  For  plane  problems,  the  displacements  u^z^.z^), 
stresses  rap(zi>z2 )  an^  resultant  force  over  an  arc  R  —  R^  4-  iP-2  may  be  expressed  in 
terms  of  two  complex  functions  W(z)  and  w(z)  viz. 


2/r(uj-t-u2 )  —  k.W(z)  -  zW'(z)  -  w(z) 


iR  =  W(z)  +  z  W"(z)  +  w'(z) 


where 


n  — 


3-4 1/ 

(3-i/)/(l  +  i/) 


plane  stress 
plane  strain 


and  /x  and  u  are,  respectively,  shear  modulus  and  Poisson’s  ratio.  For  the  class  of 
proolems  under  consideration,  the  solution  is  more  conveniently  expressed  in  terms  of  W 
and  another  complex  function  f  defined  by 
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f(z)  —  W(z)  —  z  W'(z)  -  v;(z) 

As  a  convention,  an  additional  subscript  a  is  placed  on  a  parameter  or  variable  to 
indicate  its  region  of  definition  D  Thus,  W  and  f  are  defined  on  D  for  which  the 
elastic  constants  are  v ^  and  Finally,  for  a  two-component  composite  the 

following  composite  parameters  are  important  constants 

_  (1+*2^1  *  _  1 

where  p,  are  associated  with  the  inhomogeneity 

For  a  given  problem  the  four  complex  functions  \Y\,  f^ ,  and  are  obtained 
as  infinite  series.  The  constant  coefficients  are  then  determined  by  employing  the  readily 
available  Fast  Fourier  Transform  Algorithm. 

A )  A  Crack  Lodged  in  a  Tip  Inhomogeneity. 

The  relevant  analysis  and  results  may  be  found  in  Chapter  III  of  Appendix  I.  The 
asymptotic  analysis  is  thoroughly  carried  out  to  the  order  of  e,  and  the  final  set  of 
equations  to  be  solved  numerically  are  given  by  equations  (3.21)  and  (3.22)  of  Appendix  I. 
The  implementation  of  Fast  Fourier  Transform  subroutines  is  illustrated  by  equation 
(3.28). 

Our  scheme  appears  to  be  more  powerful  than  most  of  the  known  techniques  in  that 
solutions  accurate  to  the  order  of  e  can  actually  be  computed  (c.f.  (3.29)).  Moreover, 
nonsvmmetnc  problems  can  be  handled  just  as  easily  as  symmetric  ones.  The  set  of  results 
associated  with  Figs.  3.1  and  3.2  of  Appendix  I  shows  that  the  scheme  is  not  affected  by 
the  slenderness  of  the  inhomogeneity  (a  circula  imhomogeneity  is  the  least  troublesome  in 
every  respect).  To  illustrate  the  efficiency  of  our  scheme  for  nonsymmetnc  configurations, 
a  slender  inhomogeneity  inclined  at  an  angle  u  from  the  main  crack  is  used  as  an  example, 
and  no  convergence  difficulties  were  encountered  in  the  calculation  (Fig.  3). 


7 


(B)  A  Crack  Surround  by  a  Thin  lahomogeneity. 

The  relevant  analysis  and  results  may  be  found  in  section  4.4  of  Appendix  I.  One  of 
the  presumed  properties  of  the  thin  inhomogeneity  is  that  it  shrinks  onto  the  crack  as  e-*0» 
The  most  important  and  new  analytic  conclusion  of  our  investigation  is  the  property 

of  the  solution  for  in  the  immediate  neighborhood  of  a  crack  tip.  It  is  shown  that  the 

—  1  /° 

solution  must  involve  powers  ot  r  where  i  —  0  is  the  crack  tip  which  is  in  D^. 

This  is  in  contrast  to  the  crack-tip  inhomogeneity  solution  where  the  series  is  of  the  form 
of  r-<n  +  !/2) 

This  derived  property  establishes  the  form  of  the  desired  series  solution  upon  w'hich 
the  Fast  Fourier  Transform  Algorithm  is  again  applied.  The  accuracy  and  efficiency  (in 
terms  of  the  number  of  terms  needed)  of  the  scheme  are  checked  by  the  special  case  of  a 
confocal  ellipse  for  which  a  numerically  exact  solution  is  available  (Figs.  4.2.3,  4.2.4  and 
4.3.2  of  Appendix  1). 


(C)  Eshelby  Tensor  and  Eshelby  (Configuration)  Forces. 

The  relevant  analysis  are  earned  out  in  detail  in  Chapter  V  of  Appendix  I.  Let 
f(e^/j)  be  the  strain  energy  density  function,  so  that 

t  n  —  df/de  ... 

ap  '  up 

The  Eshelby  tensor  p^  is  defined  by 

P/3a  =  f 6^  ~  T^7U7,a 

It  is  known  that  the  various  configuration  forces,  which  are  the  generalized  forces 
associated  with  damage  evolution,  are  but  integrated  forms  of  p^.  The  calculation  of  the 


many  configuration  forces  would  be  greatly  facilitated  by  a  suitable  representation  of  p 
Tills  objective  is  realized  in  Chapter  V  of  Appendix  I.  It  is  shown  that 


0a 


Pi2-P2i  =  2(W*'  +  W*  ) 


(P^TPof)  T  *2p22  —  2(z\\  +  w*  ) 


e 


where  w+  and  W*  are  two  new  complex  functions  defined  by 

=  -  J±li(W-(z))2 

w*'(z)  =  -  J~i2W  (z)  w  (z)  . 

Since  both  W  and  w  are  already  m  series  form,  the  desired  convenient  representation  is 
complete. 

3.  Grade-3  Elasticity  and  Surface  Phenomena 

The  issue  of  damage  evolution  was  deliberately  bypassed  in  the  previous  section  for 
the  simple  reason  that  there  has  not  been  to  date  a  complete  and  definitive  continuum 
theory  that  may  be  used  to  study  fatigue  damage  propagation.  At  the  same  time  it  is  clear 
that  damage  leads  to  failure  and  the  creation  of  new  surfaces.  It  is  therefore  desirable  to 
look  for  a  theory  that  actually  includes  surface  tension  as  one  of  its  variables.  The  grade— 3 
elasticity,  which  was  fully  developed  by  hlindlin  in  1965,  has  just  such  a  variable.  A 
detailed  re-examination  of  the  theory,  with  special  emphasis  on  surface  phenomena,  is 
submitted  as  the  second  accomplishment  of  this  one-year  project.  The  main  results  are 
delineated  in  a  manuscript  which  is  attached  as  Appendix  II. 

In  view  of  the  complexity  of  the  theory,  none  of  the  mathematical  deductions  and 
formulae  are  reproduced  in  this  section.  Highlights  of  our  accomplishment  include  the 
specific  determination  of  the  energy  associated  with  a  surface,  with  or  without  the  presence 
of  body  forces  and/or  other  surface  tractions.  These  results  are  presented  in  detail  in  a  self 
contained  manner  in  Sections  6  and  7  of  Appendix  II. 

The  deduction  given  in  Section  4  of  Appendix  II  has  effectively  reduced  the 
complexity  of  the  theory  to  an  extent  that  only  a  series  of  much  simpler  problems  of 
familiar  properties  needs  to  be  solved.  In  view  of  the  exploratory  nature  of  this 
inestigation,  this  accomplishment  is  most  significant  in  that  it  can  be  readily  applied  to 
yield  solutions  to  important  benchmark  problems  upon  which  the  physical  significance  of 


the  theory  could  be  meaningfully  evaluated  Point  load,  cracks  and  notches  are  among  the 
first  ones  to  be  continually  examined  bv  11s. 
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ABSTRACT 


Newiv  engineered  high— performance  composite  matenais  are  very  often  reinforced  bv 
particles,  continuous  or  short  fibers  and  thin  layers.  Cracks  encountered  in  such  matenais 
are  more  often  than  not  affected  by  its  tip  being  iocated  in  one  particular  small  particle  or 
thin  layer  of  a  composite.  The  physical  effect  of  such  apparently  small  geometnc 
alterations  on  the  toughness  of  the  material  is  finite  and  must  be  carefully  examined. 

Fatigue  crack  propagation  usually  leads  to  the  formation  of  a  thin  layer  of  damaged 
material  surrounding  the  propagating  crack.  The  thin  layer  of  damage,  however,  is  known 
to  have  finite  effects  <  ■>  the  various  generalized  Esheibv  forces  that  dnve  the  damaee. 

The  interaction  of  a  crack  and  a  small  inhomogeneitv  in  an  otherwise  homogeneous 
medium  is  studied  in  this  thesis.  Asymptotically  deduced  computer  codes  are  developed 
for  the  purpose  of  computing  any  and  all  physical  quantities  relevant  to  the  aforementioned 
problems. 


CHAPTER  I 


INTRODUCTION 

With  the  advent  of  engineered  multiphase  materials  in  recent  years,  most  notably 
composites  reinforced  by  particles,  continuous  or  short  fibers  and  thin  layers,  there  has 
come  an  interest  in  crack  problems  involving  bodies  of  spatially  varying  material 
properties.  One  particular  aspect  of  the  problem  is  the  alteration  in  the  stress  intensity 
factor  (SIF)  from  its  apparent  value  due  to  either  the  crack  tip  being  lodged  in  a  region 
with  elastic  moduli  differing  from  the  bulk  or  the  complete  crack  being  located  in  one 
particular  phase  of  a  composite.  A  crack  partially  penetrating  a  thin  fiber  and  a  crack 
situated  inside  of  a  thin  layer  are  but  two  examples.  It  is  noted  that  the  word  thin  is 
appended  to  stress  the  very  often  encounted  physical  situation.  The  results  presented  in 
this  thesis  have  direct  bearing  on  the  understanding  of  the  class  of  problems. 

Fracture  toughness  enhancement  has  been  observed  in  a  number  of  ceramic  systems 
containing  particles  which  undergo  a  transformation  of  the  martensite  type  (T.K.  Gupta 
etal  1978).  The  high  stresses  in  the  vicinity  of  a  macroscopic  crack  induce  a  transformation 
of  the  particles  and  thereby  alter  the  crack  tip  stress  field.  The  associated  SIF  could  be 
approximated  by  that  for  a  crack  tip  lodged  in  a  thin  inhomogeneity  with  moduli  softer 
than  the  bulk. 

It  is  now  a  common  knowledge  that  there  exists  a  fracture  process  region  or  damage 
zone  near  a  crack  tip  where  fracture  process  such  as  nucleation  of  voids  or  microcracks  and 
their  coalescence  take  place  and  the  usual  continuum  theory  does  not  apply.  The  damage 
zone  is  usually  small  compared  with  the  length  of  the  crack  and  may  be  approximated  by  a 
thin  inhomogeneity  with  softer  moduli. 

The  size  and  shape  of  a  process  region  change  as  it  moves  along  with  the  crack  tip 


under  fatigue  loading  conditions.  The  net  result  is  that  of  a  crack  surrounded  by  an  active 
crack-tip  process  region  together  with  an  inactive  wake,  the  damage  region  created  and  leit 
behind  by  the  traversing  active  region.  This  is  just  the  crack-layer  configuration 
investigated  by  Chuduovsky,  Moet  and  Botsis  (1987).  If  we  approximate  the  damaged 
region  by  an  elastic  material  with  softer  moduli,  the  crack— layer  configuration  is  just  that 
of  a  crack  lying  inside  of  a  thin  inhomogeneity  embedded  in  an  otherwise  homogeneous 
medium.  The  change  of  the  crack-layer  configuration  leads  to  a  number  of  identifiable 
energe-release  rates  which  cannot  be  determined  withour  a  careful  stress  analysis.  The 
results  presented  in  this  thesis  provide  effective  means  to  perform  the  needed  calculation. 

With  the  above  problems  and  the  attending  importance  and  applications  in  mind,  we 
direct  our  attention  in  this  thesis  to  the  following  two  specific  classes  of  problems: 

Problem  I  .  Disjointed  Inhomogeneities 

A  crack  of  length  2  in  an  infinite  plane  is  assumed  to  have  its  tips  separately  lodged 
in  vanishingly  small  inhomogeneities  of  size  e  (e  <  <  1),  Fig.  1.1. 

Problem  II  .  Single  Inhomogeneitv 

A  crack  of  length  2  in  an  infinite  plane  is  wholly  embeded  inside  of  an  inhomogeneity 
of  vanishingly  small  thickness  e  (e  <  <  1),  Fig.  1.2 

The  choice  of  the  basic  configuration,  a  straight  crack  in  an  infinite  plane,  enables  us 
to  remove  the  geometric  and  loading  complications,  which  may  be  handled  by  conventional 
means,  from  the  new  and  essential  asymptotic  deduction  as  well  as  the  appending 
numerical  scheme.  In  fact,  the  ultimate  objective  is  to  incorporate  the 


asymptotic/numerical  result  of  the  thesis  into  a  general  situation  for  practical  applications. 

Both  the  inhomogeneity  and  the  infinite  medium  are  assumed  to  be  homogeneous  and 
isotropic  in  this  thesis.  The  analysis  may  be  straightforwardly  extended  to  cases  where  the 
thin  inhomogeneity  is  anisotropic,  but  detailed  calculations  are  not  included  in  this  thesis. 
Similar  approach  could  be  applied  to  the  situation  where  the  inhomogeneity  is  actually 
inhomogeneous.  This  latter  consideration  requires  extensive  analysis  and  is  not  considered. 


The  solutions  to  the  two  classes  of  problems  depend,  among  other  parameters,  on  the 
size  of  the  inhomogeneity,  i.e.,  the  small  parameter  «.  It  is  clear  from  the  general  nature  of 
the  problems  that  no  analytic  and  explicit  solutions  can  be  expected.  The  implementation 
of  a  computational  scheme,  however,  must  be  preceded  by  an  asymptotic  analysis  to  get  nd 
of  the  e,  as  no  numerical  scheme  can  possibly  handle  the  conflicting  limits  required  by  e-o 
(almost  no  inhomogeneity)  and  r  -*  ®  (crack  tip  inside  of  the  vanishingly  small/thin 
inhomogeneity).  Such  an  asymptotic  analysis,  together  with  the  attending  numerical 
scheme,  represents  the  main  original  contribution  of  this  thesis.  The  final  product  are  two 
computer  codes  which,  together  with  other  codes  that  may  be  developed  to  accommodate 
geometric  and  loading  conditions,  may  be  immediately  adapted  for  application. 

To  ensure  the  correctness  and  accuracy  of  the  computer  codes,  a  number  of 
benchmark  problems  are  also  considered.  Some  of  them  are  also  new  and  original  but  the 
main  motivation  was  for  the  purpose  of  validating  the  anticipated  computer  codes. 

For  problem  I,  both  the  inhomogeneity  and  the  medium  contain  a  traction-free 
boundary  and  the  analytical  structure  of  the  solution  follow  directly  from  that  of  the 
well-known  crack— tip  field  representation.  The  standard  Fast  Fourier  Transform  routines 
are  used  to  determine  the  unknown  coefficients  in  the  series  representation,  and  the 
capability  of  generating  the  solution  to  the  order  of  e  is  established.  The  benchmark 
problem  is  that  of  a  semi-infinite  crack  penetrating  a  circular  inhomogeneity,  a 
numerically  exact  solution  obtained  by  Steif(1987). 

For  problem  II,  the  medium  contains  no  portion  of  the  traction— free  boundary  and, 
as  a  result,  the  analytical  structure  of  the  solution  for  the  medium  in  the  vicinity  of  a  tip  is 
not  clear.  For  this  reason,  the  case  of  a  crack  embeded  inside  of  a  confocal  elliptic 
inhomogeneity  is  introduced  as  a  benchmark  problem  and  studied  in  detail.  The  confocal 
gemetry  accommodates  a  Fourier  series  representation  but  the  convergence  of  the  series 
becomes  extremely  slow  as  the  inhomogeneity  becomes  vanishingly  thin.  Nevertheless  the 
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e-o  iimit  is  established  bv  extrapolation  and  a  numencailv  exact  bencnmarx  is  established. 
The  aforementioned  analytic  structure,  however,  is  still  unknown.  This  important 
information  is  revealed  by  a  detailed  asymptotic  analysis. 

The  required  asymptotic  analysis  parallels  to  that  used  in  thin  airfoii  analysis  but  is 
much  more  involved,  as  both  the  "airfoil"  and  the  outside  medium  are  each  governed  by 
two  complex  functions.  Moreover,  there  is  a  square— root  singularity  inside  the  "airfoil". 
The  analysis  is  carried  out  systematically  and  the  needed  analytical  structure  is  extracted 
from  the  inner  expansion  of  the  outer  expansion.  This  property  is  used  to  construct  the 
needed  series  solution  whose  coefficients  are  again  determined  by  employing  the  readilv 
available  Fast  Fourier  Transform  routines.  As  an  unexpected  byproduct,  an  approximate 
but  explicit  solution  is  also  obtained  for  the  confocai  crack  problem.  This  result,  together 
with  the  numerically  exact  Fourier  series  solution,  is  used  to  validate  the  final  computer 
code. 

Numerical  results  obtained  from  the  computer  codes  developed  for  the  two  problems 
are  further  validated  by  the  following  intuitive  considerations.  For  problem  II,  the  center 
portion  of  the  inhomogeneity  is  expected  to  have  small  effect  on  the  solution  as  e-o.  If  so 
the  solution  for  Problem  II  should  tend  to  that  for  Problem  I  in  terms  of  a  geometrically 
obvious  parameter.  Such  a  tendency  does  appear  to  exist.  Extending  the  iength  of  the 
inhomogeneity  of  problem  II  would  result  an  increase  in  SIF  and  its  vallue  would  always  be 
bounded  from  above  by  that  for  a  crack  embeded  in  an  infinite  layer  of  vanishing 
thickness.  Our  numerical  results  also  conform  with  this  intuitive  observation. 

Chapter  2  summarizes  the  formulation  in  terms  of  a  complex  variable.  The  rest  of 
the  exposition  follows  approximately  the  order  of  outlining  delinedted  in  this  chapter. 


CHAPTER  II 


.u 


COMPLEX  VARIABLE  FORMULATION 


Let  (zpZj)  be  rectangular  Cartesian  coordinates  and  z  =  z^  +  iz2  the  associated 
complex  variable  in  the  z-plane.  For  plane  problems,  the  displacements  u  (z^^),  stresses 
Ta0  (zrz2)  anc*  resu^tant  force  over  an  arc  R  =  +  iR2  may  be  expressed  in  terms  of 

two  complex  functions  W(z)  and  w(z),  viz., 


2**(Uj  +  iu0)  =  *W(z)  —  zW'(z)  -  w(z)  , 


where 


iR  =  W(z)  +  zW'(z)  +  w(z; 


K 


3  -  \v  plane  stress 
(3-i/)/(l+i/)  plane  strain 


(2.2 ) 

(2.3) 


and  /i  and  v  are,  respectively,  shear  modulus  and  Poisson’s  ratio. 

For  regions  containing  a  portion  of  the  real  axis  along  which  displacement  and  trac¬ 
tion  are  continuous,  the  function  w(z)  may  be  expressed  in  terms  of  W(z)  and  a  new 
function  f(z)  as  follows  (England,  1971): 


w(z)  =  W(z)  -  zW'(z)  —  f(z)  . 
Using  the  above,  we  obtain  from  (2.1)  and  (2.2) 


(2.4) 


iR  =  W(z)  +  W(z)  +  (z— z)W'(z)  -  f(z)  ,  (2-5) 


2/t(u1  +  iu2)  =  (k+1)W(z)  -  W(z)  +  W(z)  +  (z-z)W'(z)  -  £(i). 

(2.6) 

A  traction-free  crack  of  length  2  is  located  on  the  real  axis  with  |  z^  (  <1.  The  crack  is 


partially  or  wholly  embedded  in  an  inhomogeneity,  denoted  by  D,,  which  m  turn  is 

embedded  in  an  infinite  medium  denoted  by  D0)  Fig.  2.2a.  An  additional  suDsenpt  a 

wiil  be  placed  on  a  parameter  or  variable  to  indicate  its  region  of  definition.  Thus,  u  ,  v^, 

\  ,  W  and  f  are  defined  for  region  D  The  crack  is  located  in  D..  and  the  associated 
a  a  a  °  a  l 

traction-free  condition  may  be  integrated  once  to  become  (c.f.  (2.5)) 

Wj(x)  +  W*(x)  -  f^(x)  =  0  ( j x |  <  1 )  (2.7) 

where  the  notation  F±(x)  =  F(x  ±  iO)  has  been  used. 

The  infinite  medium  is  loaded  at  infinity  bv  r  a  —  it  a  so  that 

*  -  ap  ap 


=  W z  ,  fj  =  fz  as  z  -  o 


where 


W  =  4^1l  +  '  f  =  \  (*ii  “  ffr>)  +  i<ri2  - 


2  '  ll  ~  22'  "r  12 


<r  =  2W-/=«r22-i<r12 


(2.10) 


is  another  parameter  to  be  used  in  the  sequel. 

The  interface  between  and  D2  is  denoted  by  C  and  is  defined  by 

Ci  z  =  z  . 


(2.11 ) 


It  is  assumed  that  C  is  perfectly  bonded  so  that  traction  and  displacement  are  con¬ 
tinuous  along  C.  In  particular,  the  traction  continuity  condition  may  be  integrated  once 
to  become  a  continuity  condition  in  R  (c.f.  2.5).  The  two  conditions  are 

Wjl.p  +  wt(sc)  +  (*c-ic)wj(.c)  -  yic)  =  w2(zc)  +  w2(ic) 

+  (Vic)^-f2(2c), 


(2.12) 


Wl(zc)  =  ?W2(z2>  t  7*  [W2(  V  +  W2^c^  +  <zc"'c)W2^c)  "  rV  V 


(2.13) 

where  the  R-continuity,  (2.12),  has  been  used  in  simplifying  the  displacement  continuity 
conditions,  (2.13),  and 


7  ”  (1+«|)m2 


(2.14) 


are  two  composite  parameters.  A  discussion  of  composite  parameters  may  be  found  in 
(Dundurs,  1969).  Equations  (2.7),  (2.8),  (2.12)  and  (2.13)  constitute  the  governing  condi¬ 
tions  for  the  solution  of  the  desired  problem. 


We  shall  be  dealing  with  vanishingly  small  inhomogeneities  and  the  following  two 
cases  will  be  considered. 

A.  Disjointed  Inhomogeneities. 

The  crack  tips  are  separately  embeded  in  a  vanishingly  small  inhomogeneity,  Fig. 
1.1a.  The  interface  around  z=l  is  defined  by 


C:  (Rightj 


0<s<s 


o’ 


(2.15) 


where  s  is  an  arc  parameter  and  e<<x  a  small  parameter.  The  interface  around  z=— 1  is 
assumed  to  satisfy  the  condition  C(Left]=  -C(Right). 


B.  Single  Inhomogeneity 

The  crack  is  wholly  embedded  in  a  vanishingly  thin  inhomogeneity  Fig.  1.2a.  The 
interface  C  is  given  by 
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,  x  ±  iey*(x)  |x|  <  1 

C:  z  =  z  (x;e)  =  < 

1  e|{(s)  -r  ids'll  I xj  ~  1 


(2.16) 


where  s  is  a  conveniently  chosen  arc  parameter  that  may  be  expressed  in  terms  of  x,  and 
e<<l  is  again  a  small  parameter. 

The  solutions  to  the  associated  boundary  value  problems  must  therefore  depend  on 
the  size  of  the  inhomogeneity  e.  We  shall  use  the  genenc  symbol  F(z;e)  to  indicate  the 
dependence  of  F(z)  on  e.  The  objective  of  this  thesis  is  to  obtain  the  asymptotic  expansion 
of  F(z;e)  as  e  -  0.  We  note  that  the  solution  characterized  by  u(z)  =  u^  +  iu^  satisfies  the 
condition 

u( — z )  =  — u(  z )  2.17) 


for  both  cases. 
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CHAPTER  IH 

A  CRACK  WITH  TIPS  LODGED  IN  VANISHINGLY  SMALL  INHOMOGENEITIES 

3.1  Introduction 

For  the  case  of  a  finite  crack  with  tips  lodged  in  vanishingly  small  inhomogeneities, 
the  asymptotic  limit  is  just  the  solution  for  a  semi-infinite  cracir  lodged  in  an 
inhomogeneity  of  finite  size.  The  benchmark  problem  is  that  of  a  semi— infinite  crack 
penetrating  a  circular  inhomogeneity  which  was  solved  by  Stief(1987).  NoT:circular 
inhomogeneity  cases  have  been  considered  by  Hutchinson  ^1986).  Dimension— analysis 
considerations  indicate  that  the  asymptotic  limit  can  only  depend  on  the  shape  of  the 
inhomogeneity,  as  far  as  geometric  dependence  is  concerned.  This  dependence  is  fully 
accounted  for  in  our  calculation  via  the  use  of  the  readily  available  Fast  Founer  Transform 
Algorithm.  Moreover,  the  size  effect  may  be  obtained  by  including  an  additional  term  in 
the  asymptotic  expansion. 

A  recap  of  the  formulation,  together  with  the  introduction  of  a  boundary— layer 
complex  variable  (  ,  is  given  in  section  3.2.  Outer  expansion,  which  is  valid  away  from  the 
crack  tips,  is  presented  in  section  3.3  in  terms  of  z.  Inner  expansion  of  the  outer  expansion 
is  then  used  as  a  guide  to  establish  the  inner  expansion,  which  is  presented  in  section  3.4. 
The  application  of  the  Fast  Founer  Transform  Algorithm  to  the  system  of  interface 
boundary  conditions  is  discussed  in  section  3.5.  Some  numerical  results  are  presented  in 
section  3.6  mainly  for  the  purpose  of  illustrating  the  efficency  of  the  numencal  scheme. 


3.2  Formulation 

Using  (2.15)  and  (2.5)  ,  the  resultant— tree  condition  along  the  crack  become  (c.t. 
Fig. 1.1a): 


W*(x;s)  +  W*(x;e)  -  f^(x;e)  =  0  for  |x|  s  1 

(3.1) 

W2(x;t)  +  W2(x;e)  -  f^(x;e)  =  fl  for  |x|  <  1 

(3.2) 

The  resultant-continuity  and  displacement-continuity  conditions  along  C  are  : 


w2( V)  +  w2^c:£^  +  ^zc"'c^W 2  (  *  c  1  £  }  ~  y "c;£) 

=  W1(2c:e )  +  Wi(2cif)  +  (2c-ac)Wi  (  2  c  ;  f  “  fi(2cu). 

W1(z(.,e)  =  7W2(zc;e) 

+  7*[  W2(ze-<)  +  W2(ic;€)  +  (*c-ic)W2  (  =  c  .  <  )  -f2(V«)] 


(3.3) 


(3.4) 


where  (3.3)  has  been  used  in  deducing  (3.4).  The  loading  condition  at  infinity,  (2.8)  ,  is 

W2(z,e)  =  Wz,  f 2 ( z ; e )  =  fz,  as  z  -  ®.  (3-5i 

As  e  -  0  the  conditions  (3.1),(3.3)  and  (3.4)  disappear  altogether  .  Borrowing  the 
terminology  used  in  thin-airfoil  analysis  (Van  Dyke  1975)  we  shall  call  the  attending 
solution  the  outer  expansion.  Such  an  expansion  is  not  valid  near  |z|  =  1.  Inner 
expansions  near  the  crack  tips  must  be  constructed.  For  this  purpose  a  boundary  layer 
complex  variable  (  is  introduced  for  the  neighborhood  containing  z  —  1  .  It  is  defined  by 

/  i*  z  -  1 
(  =  P  e  f  =  — — 


(3.6) 


.0 


and  the  interface  C  [nghtj  may  oe  conveniently  written  m  the  form 


z  -  i  -  to  i  f)  ,  ‘  =  =  p f)  e 


■f 


1  3.1 


where  the  single  function  p  iv)? j  may  be  used  to  define  the  vanisningiy  smail  inhomogeneitv 
(Fig. 1.1b).  It  is  assumed  that 

and  (3-8> 

so  that  p^j  charactenzes  the  maximum  dimension  of  the  inhomogeneity.  The  quantites  p^ 
anc^  ^m  are  ^ey  tealures  °f  function  p^if).  Finally,  the  square— root  character  of  a 
crack  tip  suggests  that  the  required  asymptotic  sequence  must  involve  powers  of 

3.3  Outer  Expansion  ( a  ^  *  1 ) 

The  outer  expansion  for  W2  and  f2  are  gonvered  by  (3.2)  and  (3.5).  They  are 


L(z;e)  -  fz  + 


n=l 


Att(«)  An(f) 


2-1 


z_+J_l 
e  p 

L 


(3.9) 


2W2(z;e)  -  f2(z;e)  ~ 


,  2  .  ,1/2 

<r(z  -1)  ' 


CD 

1  +  X 

n=l 


Ba(0 


Bn<«> 


z-  1  " 

□ 

z  +  1 

L«'nJ 

(3.10) 


where  a  is  given  by  (2.10)  and  W  and  f  by  (2.9),  and  the  coefficients  A  (e)  and  BQ(e)  have 
representations  of  the  forms 


A „(<)=’>'  A  <<W)m 
1  t'm4  am'  1 
m=0 


B  (e)  =  S  B  (e 
nv  ’  t'm4  am' 

m=0 


1/2 


m 


(3.11) 
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We  shall  be  needing  the  inner  expansions  of  the  outer  expansion.  These  resuits,  after 
normalized  by  an  as  yet  undefined  factor  a*  j y'S,  are 


f2(1  +  e(iO  ,/2  a  ®, 

- - +  «/2 -  >  A.U'/pJ 

//A  /  /  n0  m 


&  ' 

tr  *  m  m 


+  f3/2'  I  An2^m)  n  +  — A10  '  — +' 

n=l 


(3.12) 


2W2(l+e(;e)  -  f2(l  +  e<;e) 
<r*/  fl 


’  1/2«  )1/2  1  +  ’£  Bn0«/om)  n 

n=l 


+*«  )1/2  :>>nlu/o-n  - 


n=l 


+<3/v/2  y  B„,«//.mrn + 4-  b. 


n2vw,m' 

n=l 


1  ®  0 

+— «  )  i+y-r 

4  — -  (n 

L  n  1  > 


(3.13) 


The  above  expansions  reduce  to  the  ordinary  crack  solution  as  t  -  0  ana  the  Laurent 
series  with  unknown  coefficients  are  needed  to  account  for  the  as  yet  unknown  effects  of  the 
vanighmgly  small  inhomogeneities  located  at  z  =  *  1.  Similar  situations  have  been 
encountered  in  crack  branching  (Wu,  1978)  and  crack  tip  contact  (Wu,l982).  That  the 
Laurent  series  at  a  =  *  1  have  same  coefficients  is  a  result  of  (2.17).  The  sigulanties  are 
contained  in  and  have  no  effects  on  the  analyticity  of  Wj  and  {2  which  are  defined  in 


3.4  Inner 


nons  Near  *  =  1. 


For  a  generic  function  F(z;e),  the  associated  inner  expansion  is  defined  by 


?•((;«).  F(l+«<;0  / 


(3.14) 


Moreover, 


4  mm 


(3.15) 


It  follows  that 


Y  f,M). 


(3.16) 


UJ 

-  <?(<;«)  ~  ]>>1/2)n  l2WL«)  -  ■ 


(3.17) 


where  the  right-hand  side  terms  are  just  the  terms  defined  in  (3.12)  and  (3.13). 


The  functions  f^((;e)  and  W*((;e)  must  satisfy  (3.2)  in  the  appropriately 
transformed  form.  This  condition  is  identically  satisfied  if 


0) 

2w;((;0-^(i;.)~  y  (e1/2)n  [jw;tt(0  -  fjn(i) 


(3.19) 


f10«)  =  3f /  ~~  .  2Wj0(o-^0U)=o. 


I.  ak,n-I  M)k  ■ 
k=0 


wi.co-S.m-d'T. 


iL2bk,n-l  ’ 

k=0 


(3.20) 


where  the  last  two  expressions  are  for  n  >  1.  All  the  unknown  coefficients  must  now  be 
determined  to  satisfy  (3.3)  and  (3.4)  which  now  become 

"Liu  +  w2„y>  +■  «c-  yw;0uc)-g!ci 


=w;n(y+w;n(yM(c-yw«ln(.c)-^c) , 

WIn«c>  =  3Wlnlic>  +  3*[  WlnUc)  +  W^) 

+<<c-yw2n(cc)-r2n(y  - 


(3.21) 


(3.22) 
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The  expansions  defined  by  (3.16)  (3.17)  (3.12)  and  (3.13)  are  valid  for  the  region 
outside  the  inhomogeneity,  Fig  1.1b,  and  the  unknown  coefficients  involved  in  (3. 12), (3. 13) 
and  (3.20)  are  determined  by  applying  the  standard  Fast  Fourier  Transform  (FFT) 
algorithm  to  (3.21)  and  (3.22). 

Before  implementing  the  FFT  calculation,  the  normalizing  factor  c  should  be 
properly  chosen  for  specific  cases.  This  is  explained  in  the  following  section. 


3.5  Fast  Fourier  Transform  Algorithm  And  Numerical  Procedure. 

It  is  convenient  to  consider  shear  loading  and  axial  loading  (^^1^22^  seParate^Jr- 
Let  us  first  consider  the  case  ^2=^'  a22  anc*  ^rst  ^  sets  °*  un^nowT1 

sequences  become: 


4,(0  =  ■n  1  /  »*.  2w;0(o  -  4,10=0 

4,(0  =  74,10.  2w*0«)  -  q0to=o 


(3.23) 


t  r2i(o=X  An0((/pmr" 


n=l 


2w*j(O-f*21(0  =  2((  )1/2 


© 


1  +  L 

n=l 


,-n 


*11^  'I  anO^^M^ 

n=l 


2wi,(o-4(o  =41/2,r  2b„o((/<’M)° 


Jim<  n0^lHW 
n=0 


(3.24) 
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Each  set  must  satisfy  the  continuity  conditions  (3.21)  and  (3.22).  We  note  that 
(3.21)  is  the  continuity  in  resultant  force.  The  first  set  (3.23)  satisfies  the  conditions 
identically.  For  the  other  sets  we  truncate  the  sums  at  n=N  so  that  a  total  of  4N+2 
unknown  coefficients  are  involved.  Since  the  standard  Fast  Founer  Transform  subroutines 
FFTCF  and  FFTCB  are  written  for  functions  of  a  real  variable,  complex  functions  of  the 
complex  variable  (  —  p  typified  by  a  generic  symbol  G((  ),  are  represented  as 

C  C  v# 

follows: 


G«c) 


N  N 

Tcn^=T 


n=l 


n=l 


M 


-M 


')'c  .“'t'j'c 

i  mn  i  L*  mn 


m=0 


m=0 


(3.28) 


where  M  =  500  has  been  used  in  all  calculations.  The  truncation  number  N  is  determined 
by  a  convenient  convergence  criteria. 

We  normalize  the  physical  stress  intensity  factors  Kj  and  Kjj  by  a22^*  anb  tbe  resu^t 
is  : 


(1)  (1) 

K:  -iK2  =(Kj-iK  n)l*22fi 
~  b00  +  fl^b01  +  eb02  + 


(3.29) 


It  is  clear  from  (3.25)  and  (3.27)  that  bg^  =  0  if  ^ 

For  a  pure  shear  loading  condition,  =  0  and  we  choose  a  —  ~iai2  ' 

Equations  (3.23)  -  (3.26)  remain  unchanged  and  (3.27)  is  replaced  by 


(3.30) 


00 


The  associated  stress  intensity  factors  are  given  by 


(2)  (2)  ,  ^ 

K2  +  i  Kx  =  (  Kn  +  i  Kj  )  /  <rl2y? 


b00  +  fl/2b01  +  eb02  + 


(3.31) 


It  is  clear  from  the  solution  procedure  that  the  solutions  for  (3.24)  and  (3.26)  are 
identical  in  form  for  both  loading  conditions.  The  solution  for  (3.25)  is  determined  by 
either  (3.27)  or  (3.30)  Thus,  bg^^  is  never  zero  in  (3.31).  On  the  other  hand,  bgg  and  bg2  for 
both  cases  are  identical.  The  actual  field  variables,  however,  must  be  modified  by  the 
factor  j*  and  hence  are  completely  different  for  the  two  cases. 


3.6  Results  and  Discussion 

As  we  have  mentioned  in  the  introduction  that  the  e  -*  0  limit  depends  only  on  the 
shape  of  the  inhomogeneity,  as  far  as  geometric  dependence  is  concerned.  To  illustrate  this 
dependence  and  also  to  test  the  efficiency  of  our  computer  code,  the  class  of  problems 
depicted  in  Fig.  3.1  is  considered.  The  asymptotic  configuration  is  that  of  a  semi-infinite 
crack  with  a  tip  inhomogeneity  of  thickness  2  and  length  2  +  a+  b.  We  mention  in  passing 
that  the  actural  physical  dimensions  are  2e  *  (2  +  a  +  b)  e.  The  material  properties  are 
fixed  by  =  0.2  and  /^//*2  =  0-5  80  the  parametric  study  is  purely  geomotrical. 

Moreover,  onl>  mode— I  loading  is  considered  in  the  illustration. 

The  benchmark  situation  of  a  circular  inhomogeneity  is  recovered  by  setting  a  =  b  = 
0  and  the  associated  normalized  SIF  is  0.645  (Steif  1987).  It  is  anticipated  that  increasing 
a  would  lead  to  an  decrease  in  SIF  while  increasing  b  would  actually  intensify  the 
associated  SIF  (Hutchinson  1986).  Every  method  has  its  limitations  and  the  present  code 
cannot  be  expected  to  function  for  cases  where  a  and  b  are  very  much  greater  than  1. 
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Considered  as  a  function  of  a  and  b,  the  normalized  SIF,  K  (a,  b),  is  expected  to  approach 
to  certain  asymptotic  limits  very  rapidly.  K  (m,  b)  is  just  the  solution  for  problem  II  and 
K  foD,  oj)  ~  0.74  which  is  extrapolated  from  the  results  given  in  (Hilton  and  Sih,  1970). 

Tne  pertinent  results  are  protted  in  Fig.  3.2  as  a  function  of  b  with  a  as  a  parameter. 
The  a  —  cd  curve  is  produced  by  our  second  computer  code  which  will  be  developed  in 
Chapter  IV.  It  is  seen  that  the  results  conform  with  all  the  expected  trends. 
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CHAPTER  IV 

A  CRACK  IN  A  THIN  INHOMOGENEITY 


4.1  Introduction 

The  present  problem  differs  from  that  of  Chapter  III  in  that  the  complete  crack  is 
wholly  surrounded  wy  the  inhomogeneity.  As  a  result,  the  infinite  medium  contains  no 
portion  of  the  traction— free  boundary  ..nd  hence  the  analytical  structure  of  the  solution  for 
the  medium  in  the  vicinity  of  a  crack  tip  is  not  known.  For  this  reason  and  a’  >  the 
purpose  of  setting  up  a  benchmark  for  validation,  the  problem  of  a  crack  imbeded  insiae  01 
a  confocal  elliptic  inhomogeneity  is  first  considered  in  detail. 

A  Fourier  series  solution  is  first  constructed  in  section  4.2.  The  analysis  is  preceeded 
by  a  brief  discussion  of  the  associated  antiplane  shear  problem  for  which  the  exact  solution 
is  presented.  The  exact  asymptotic  limits  for  very  large  and  very  small  inhomogeneities 
are  deduced  from  the  exact  solution.  Their  exact  dependence  on  size  and  moduli  serves  as 
a  guidence  for  the  desired  plane  problem  solution. 

The  Fourier  series  solution  obtained  for  the  plane  problem  is  considered  to  be 
numerically  exact,  as  the  convergence  can  be  easily  checked  by  varying  the  number  of 
terms  included  in  the  actual  computation.  The  convergence  of  the  senes,  however, 
becomes  extremely  slow  as  the  inhomogeneity  becomes  vanishingly  thin.  The  asymptotic 
limit  for  the  SIF,  though,  is  easily  extrapolated  from  the  numerical  results. 

The  series  solution  sheds  no  light  on  the  desired  analytical  structure  of  the  solution. 
A  complete  asymptotic  analysis  is  therefore  carried  out  in  section  4.3.  The  needed 
analytical  structure  is  revealed  by  the  inner  expansion  of  the  outer  expansion.  As  an 
unexpected  useful  byproduct,  an  approximate  but  explicit  solution  is  obtained  for  the 
confocal  situation.  This  completes  the  establishment  of  the  benchmark  solution. 


The  general  <  ase  is  finally  presented  in  section  4.4  where  the  tmckness  of  the 
inhomogeneity  and  the  nose  thickness,  i.e.,  the  length  of  the  inhomogeneity  minus  *he 
crack  length,  are  assumed  to  have  the  same  order  of  magnitude  e.  We  note  that  tor  the 
confocal  situation  the  nose  thickness  is  of  the  order  of  e4"  while  the  innomogeneity  thic  me's 
is  of  the  order  e.  In  this  regard  and  in  the  context  of  the  final  boundary  layer 
computation,  the  confocal  situation  is  even  more  difficult  than  our  general  case,  as  the 
former  requires  the  satisfying  of  boundary  conditions  specified  on  a  parabola  in  the 
boundary— layer  variable.  Neverthless,  the  analytic  structure  revealed  by  the  benchmark 
analysis  is  the  key  to  the  success  of  the  general  representation.  The  unknown  constents 
involved  in  the  senes  solution  are  again  determined  by  the  Fast  Founer  Transform 
Algonthm. 


23 


4.2  Series  Solution 

4.2.1  Introduction 

A  family  of  confocal  ellipses  may  be  characterized  by  a  single  parameter  p  >  1.  In 
the  limit  as  p  •*  1  the  ellipse  degenerates  into  a  straight  line  of  length  2.  It  tends  to  e 
C’-cle  of  infinite  radius  as  p  ■*  ®.  The  geometry  of  the  problem  is  fixed  by  the  crack  ( p  = 
l)  and  the  size  of  an  elliptic  inhomogeneity  (p  =  p  ).  Both  the  mhomogeneity  and  the 
infinite  medium  are  assumed  to  be  -omogeneou-  and  isotropic. 

Section  4.2.2  summarizes  the  formulation  in  a  transformed  complex  plane.  The 
exact  solution  for  the  anti— plane  shear  case  is  presented  in  Section  4.2.3.  Explicit 
asymptotic  limits  for  large  and  small  inhomogeneities  are  extracted  from  the  exact  formula. 
Plane  problems  are  dealt  with  in  Section  4.2.4.  The  case  of  equal  shear  modulus  and 
unequal  Poisson’s  ratio  is  solved  exactly,  and  the  general  case  is  presented  as  r  senes 
solution. 

A  large  number  of  references  on  inhomogeneity  problems  may  be  found  in  Mura 
(1982,1988),  but  we  have  not  found  any  reference  dealing  with  the  consideration  of  a  crack 
in  a  vanishingly  small  inhomogeneity.  The  closest  situation  is  the  one  given  by  Warren 
(1983),  who  considered  the  edge  dislocation  inside  an  elliptical  inhomogeneity,  including 
vanishingly  small  inhomogeneities. 

Numerical  results  for  the  plane  problems  are  presented  only  for  plane  strain  arid 
Mode— I  conditions.  Parametric  dependence  of  SIF  on  the  size  of  the  inhomogeneity  is 
discussed  in  detail  in  Section  4.2.5  for  the  range  where  the  inhomogeneity  is  softer  than  the 


medium. 
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1.2.2  Formulation 

Let  (z^Zj)  be  rectangular  Cartesian  coordinates  and  z  =  z.,  +  iz0  the  associated 
complex  variable  in  the  z-plane.  A  crack  of  length  2  in  the  z-plane  is  mapped  onto  a 
unit  circle  in  a  new  complex  (-plane  via  the  mapping  function 


=  m(()  =  2  (f  +  f] 


(4.2.1) 


where  (  =  +  i^  —  P^  (Fig-  4.2.1b)  and  p  >  1.  The  image  of  the  circle  (  —  (Q  — 


i  6 

pQe  is  the  ellipse  (Fig.  4.2.1a) 

(zT/a)2  +  (z2/b)2  =  1 

where 


1 

l 

,  i 

' 

l 

a  -  2 

r°  v 

•  b"2 

{'*-4 

(4.2.2) 


(4.2.3) 


The  infinite  z-plane  is  now  conveniently  divided  into  two  regions  and  Dj  by  the 


single  parameter  pQ,  viz., 


Di:  1  <  P  <  .  Dr  P  >  Pn 


(4.2.4) 


We  shall  call  the  inhomogeneity  and  the  medium.  The  two  regions  are  of 

different  elastic  materials  characterized  by  shear  moduli  p  and  plane— elasticity  constants 


K 

a 


3-4i/fl  plane  strain 

.  (3-l/a)/(l-bi/a)  plane  stress 


(4.2.5) 


where  v  are  Poisson's  ratios.  The  infinite  plane  is  loaded  at  infinity  by 


Anti-plane:  as  |z|  -•  ®  , 

Piane:  Ta0=  %0  as  !zl  -  «  - 


(4.2.6) 

(4.2.7) 


where  r  •  are  the  stress  components. 


For  the  anti  piane  problem,  the  displacement  u^z^Zj),  stresses  r3a(zpz?i  ana 
resultant  force  along  an  arc  may  be  expressed  in  terms  of  a  smgie  complex  function 
F(z).  We  have 


»3-j[*«) +  ♦«)). 

(4.2.8) 

r3l  ~ir32  =  M*'(0/m'(0  , 

(4.2.9) 

R3  =  -i$  [♦(()  + ♦(«]  , 

(4.2.10) 

where  $(()  =  F(m(())  and  m(()  is  the  mapping  function  (4.2.1). 

For  the  plane  problem,  the  displacements  ua(zpz2)>  stresses  ra^(zi>zo)  and  resul¬ 
tant  force  +  1R2  =  R  along  an  arc  may  be  expressed  in  terms  of  two  complex  func- 


tions  W(z)  and  w(z). 

We  have 

2m(u1  +  iu2) 

=  «Q(0-— M  0'(0-M  , 

(4.2.11) 

1R  =  0(0  + 

■ m(^  -Fto-lio. 

(4.2.12) 

m'  (0 


where  0(0  =  W(m(0)  and  u(()  =  w(m(()). 

The  complex  functions  must  be  determined  for  the  two  regions  and  D2  sub¬ 
jected  to  the  loading  conditions,  (4.2.6)  and  (4.2.7),  and  the  continuity  conditions  along  the 
interface  boundary  characterized  by  p q.  The  crack  surface  is  assumed  to  be  traction  free. 
The  traction-free  condition  may  be  integrated  along  the  crack  to  become  a  resultant— free 
condition.  Similarly,  traction  continuity  along  the  interface  may  be  integrated  to  become  a 
resultant  continuity  condition.  The  integrated  forms  of  these  conditions  will  be  used  in  the 
calculations  to  follow. 

We  shall  place  a  subscript  a  on  a  complex  function  to  indicate  its  region  of  defini¬ 
tion.  For  example,  Ffl(z)  and  *a(()  are  defined  for  region  D  . 


4.2.3  Anti— Plant*  Shear 


The  problem  may  be  most  conveniently  solved  in  the  plane.  The  integrated  trac¬ 
tion-free,  integrated  traction  continuity,  displacement  contiruty  and  loading  conditions 
are 


-  *1(e1#)  =  0  , 

(4.2.13) 

(4.2.14) 

*l«o>  +  *l«o>  =  #2(40)  +  ♦,«„)  , 

(4.2.15) 

#2  ~  $  C  as  ®  , 

(4.2.16) 

where  (  =  p  el$  and 

0  0 

A  1  1  ,  \ 

*“5^2^31  ~  1<r32' 

(4.2.17) 

The  solution  is 

*j(0  =  A<  +  ^  (i  <  id  <  i>0) , 

(4.2.18) 

♦2(0  =  +  [(1+/.JJA - p>l]  >  (|(|  >  pa) , 

where 

(4.2.19) 

A  =  2#/  K+1)  +  ^  (Oj-l)  • 

(4.2.20) 

The  stress  intensity  factor  may  be  readily  determined.  It  is  convenient  to  normalize  the 
SIF  Kjjj  by  the  factor  ,  and  the  result  is 


k,  „  A  =  ill  l 

l 0  “2J  *32a  “2  ,j+i  t  +  ^  ,*_i 

1  +  "2  ',o+1 


The  following  limits  may  be  easily  obtained 


Kj®,-± 


Mi  Ai 

=  2  —  /  1,  —  , 

**2  [  >>J 


(4.2.21) 


(4.2.22) 


(4.2.23) 


These  two  limits  are  plotted  in  Fig.  4.2.2.  The  sign  of  dK ~j dp  is  governed  by  (1  — 

j  0 

M1/m2)-  Thus 


K  fl  ill  ^  K  "'l  k-  f  'll  ..  "1  <  , 

3[  ’,2J  >  3(  °'^j  >  'f  ^  ? 


(4.2.24) 


It  is  noted  that  the  bounds  are  of  practical  significance  for  the  cases  where  the  inhomo¬ 
geneity  is  softer  than  the  medium. 

A  very  slender  inhomogeneity  may  be  defined  by  p  =  1+e  where  e  <<  1.  The 

o 

exact  result  (4.2.21)  may  be  used  to  obtain 


Mi  Mi  Mi 

K3  1+e,  —  ~  —  1+e  1 - +  •  •  • 

3  m2  ^2  ^2 


(4.2.25) 


provided  that  -*  0  as  e  -»  0.  It  is,  therefore,  clear  that  such  a  two— term  asymp¬ 

totic  expansion  is  valid  only  for  the  cases  where  the  inhomogeneity  is  either  softer  or  slight¬ 
ly  harder  than  the  matrix. 
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4.2.4  Plane  Prahlem 

Since  the  plane  problems  cannot  be  solved  exactly,  we  begin  by  constructing  a  senes 
solution  in  the  (—plane.  The  traction— free  condition  on  the  unit  circle  enables  us  to  intro¬ 
duce  the  stress  continuation  (England,  1971) 

^(0  =  -Q1(l/n~  —  ■  (4.2.26) 

®'(0 

The  function  fl^(()  is  now  extended  to  the  region  1  /pQ  <  p  <  p  and  the  traaion-free 
condition  on  p  =  1  is  identically  satisfied.  The  conditions  (4.2.7)  are  met  if 

n2(0=Q<  +  oj|]  (4.2.27) 

«2(0  =  «4  +  o[^]  K |  -  CD  ,  (4.2.28) 

where 

fl  =  g(an  +  <r22)  ,  «  =  5  [(<^22  “  ^n)  +  l2ffi2  •  (4-2-29) 

The  integrated  form  of  the  traction  continuity  condition  along  (  may  be  obtained 
from  (4.2.12),  i.e., 

W  +  -_ ~o)  fl2(g  +  gg 
m'(g 

=  +  Mg—  f  (4.2.30) 

.°J  1  NJ  “'(g 

where  (4.2.26)  has  been  applied.  Continuity  in  displacements  along  (  yields 
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mU'J 


"A(y +  fli 


l 


m{(Q)-m 


*{(Q 


m'((J 


which,  after  applying  (4.2.30),  may  be  reduced  to 


a1(g  =  Tfl2(C0)  +  7,n2«.)  + 


m((  )  -  - 

=^W-“2  «o> 


m'((0) 


where 

(i+^l 
T  =  (T+^ 


r  ^ 


are  two  composite  parameters  (Dundurs,  1969).  We  note  that  (4.2.30)  may 
from  (4.2.31)  via  the  relation 


(4.5)  =  Letting  ^  =  1  and 


Before  proceeding,  we  shall  first  consider  the  special  case  ^  ^ 


(4.2.31) 


(4.2.32) 


(4.2.33) 
be  deduced 


(4.2.34) 
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i)  Exact  Solution  for  Equai  Shear  Modulus 

For  this  case,  the  composite  parameters  defined  by  (4.2.33)  become 


_  1+*2 

^  7o  1+^  ’ 


7*  =  0 


and  (4.2.32)  becomes 

<W  =  to«2U0)  , 

which  serves  as  an  analytic  continuation  of  the  two  functions  fi.(()  and  ^(O- 
the  substitution  fl^(()  =  we  obtain  from  (4.2.30) 

f  > 

>o; 


*2(/W  + 


m<U 


m  ( eJ/?0) 


Wi0]-r0n(i0»j;y 


=  (To-ijiyg  - 


where 


m((  )  -  ml/?  ) 

M(g  =  ■  ■°4j_ — 

m/(g) 


p20(p2Q-')((l  -  P2n) 


c0«l  - 


It  follows  from  (4.2.37)  that  the  function  H(()  defined  by 

t0o2U/pJ)  ■ 

H(0  =  ‘ 


<  I  >  o 


«2(^/g  +  _mU1  w  -  70M(on2(gg 


m'(p  W 


(  K  I  <  Pr) 


(4.2.35) 

(4.2.36) 
Making 

(4.2.37) 

(4.2.38) 


(4.2.39) 


is  holomorphic  in  the  whole  (-plane.  Moreover,  its  properties  at  (  =  0  and 
governed  by  the  right-hand  side  of  (4.2.39).  The  complete  solution  is 


<d  are 
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H(0  = 


-1  +  7 


1 

^0 


a<  +  +  [1  +  70(Pq-1)!^}  i 


Q2(0  =  fl(- 


a  + 


5/J2 

0 


1  +  T0(Pq- !) 


(4.2.40) 


(4.2.41) 


»j(0='H(/.J/f)  +  ToM(/.J/J)OJ(C)-  [m(p'/J)/m'(C)]  0j«)  ,  (4.2.42) 

0j(0  =  7o0j(0  •  (4.2.43) 

Let  us  use  the  factors  and  to  normalize  the  SIF’s  KT  and  KTT,  and  wnte 


I 


"II* 


Ki  —  Kj/<t22>/?  ,  Kj  —  Kjj /ai2^*  • 


(4.2.44) 


The  following  explicit  results  are  readily  obtained 


TqI^+1)  4-7X-D1  70(1-T0)(^-1)  TU 

2I1  +  70(pJ-l)|  2(1  +  70(PJ-1)I  '22  ' 


1 


+ 


70Pp 

70^o-1) 


There  are  the  following  exact  limits 


l+*« 

lim  Kj  and  K2  =  ,0=  *, 


lim  K|  =  j  (1  +  7#)  -  j  (1  —  7q)  111 

P0*m  22 

lim  K0  =  1  . 

V* 


(4.2.45) 


(4.2.46) 


(4.2.47) 


(4.2.48) 


(4.2.49) 
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ii)  Series  Solution 

The  complex  functions  and  together  with  their  regions  of  definition, 

admit  the  following  senes  representations 

CD 

fll(0  =  X  "TI1[Anin  +  V_nj  '  I4-2-50) 

n=l 


ho = +  X  '’o+°iy  ■ 

n=l 


(4.2.51) 


,(o  =  .<  +  y 


Am* 

n=l 


(4.2.52' 


where  n  =  1,3,5,...  The  factors  p*  Q  and  Pg+n  are  included  for  convenience  and  pose 
no  restrictions  on  the  validity  of  the  series  representations.  They  are,  nevertheless, 
conceived  from  ihe  fact  that  ufl  along  |  ( j  =  pQ  must  be  of  the  order  of  pQ  as  pQ  -  a>. 

Substituting  the  above  into  (4.2.31),  setting  (  —  p  e1^,  and  equating  coefficients  of 

in#  , 

e  to  zero,  we  obtain 


— n 


i~L 


An  +  (n+2) 


1-1 


n-f  2 


1  + 


n 


1 

+  — 


*1 

1  4-  .  .  .  L 

a  i  *  - 

Ti  1  D 

K0 

an+2  +  p2  *2 

n  a  n+2 

po 

'  r\ 


for  n  =  1 


(4.2.53) 
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and 


' 

’ 

' 

' 

1 

*1 1-  — 

A1  + 

1.  — — 

A  1 

A1  T  ~ 

1 - — 

p2 

°. 

Pl 

P2 

v Q 

"O 

P0 

+  ^i 


-B  -b  1--B1 

Pi  1  1  p\  1 


^2 


+  —  u 

p2 

Ko 


*1  ~  “3n 
P 


1 

? 

o 


A„  +  — 

2  P2 
y0 


a,  + 


* 

A  1  11 

An-2  +  T5 

1-i 
_  2 

n-2 

0 

^  o 

-  fin  3 - 

"2  »2o 


(4.2.54) 


' 

• 

1  -  — 

(n-2lB  ,  +  —  B  -  b  +  —  bn  n 

plo 

n-2  p2 

K  ■  n-2  2  n  n  2  n-2 

po  po 

for  n  =  3 

for  n  =  5,7,9,. 


(4.2.55) 


The  relations  derived  from  (4.2.30)  are  obtained  by  the  substitution  (4.2.34),  i.e., 

[Letting  n^=  ^2  =  l  and  —  *2  =  -1  *n  (4-2. 53), (4. 2. 54), (4. 2. 55))  . 

(4.2.56) 


The  infinite  system  (4.2.53 ) — (4.2.56)  are  truncated  and  the  resulting  finite  system  is 
inverted  numerically.  The  SIF’s 

KI-iKn  =  2(r)1/2n|(l)  (4.2.57) 


are  then  normalized  in  accordance  with  (4.2.44).  Values  of  plane— strain  are  plotted  in 

Fig.  4.2.3  for  the  case  —  i/j  =  0.2  and  =  0.  The  numerical  results  approach  to  a 
limit  very  rapidly  as  pQ  -  ®.  The  analytic  expression  for  this  limit  is  determined  in  the 
next  sub-section. 
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The  convergence  of  the  numerical  scheme  becomes  extremely  slow  as  pQ  -  1-  For 
this  reason,  values  of  K1  for  fixed  values  of  are  plotted  as  functions  of  1  j p20  in 

Fig.  4.2.4.  It  is  seen  that  all  curves  tend  to  finite  limits  as  pQ  -  1-  The  pQ  =  1  curve 
indicated  in  Fig.  4.2.3  is  extrapolated  from  Fig.  4.2.4. 


iii)  Very  Large  Inhomogeneitv 

As  we  have  indicated  before  that  the  factors  p^±n  built  in  (4.2.50)— (4.2.52)  are  con- 

o 

ceived  from  the  fact  that  the  displacements  ufl  along  |  (|  =  pQ  must  be  of  the  order  of 
pQ  as  pQ  -  <d.  Thus  the  series  representation  may  also  be  interpreted  as  an  asymptotic 
expansion  for  large  p  .  In  this  interpretation,  however,  the  four  sets  of  constants  must  be 
re— expanded  in  powers  of  p'Q2,  i.e., 

(A,a,B,b)a  =  (  )m  +  i  ()„!+••  •  <4-2-58> 


Substituting  the  above  into  (4.2.53) — (4.2.56)  and  equating  to  zero  the  coefficients  of  powers 
of  p'Q2,  we  obtain  an  infinite  system  of  infinitely  many  equations  for  the  determination  of 
the  coefficients  (  )  .  The  asymptotic  limit  for  the  case  of  a  very  large  inhomogeneity  is 

thus  governed  by  the  coefficients  (  )  . 

The  system  governing  (  )  actually  decouples  into  finite  systems  and  the  first  clus¬ 
ter  of  equations  are 


,-s  a  ,M1  p  _M1- 

+  3A30  a10  +  "2B10  ~  ' 


(4.2.59) 


*10  +  3*30  a10  B10  “  u  • 


(4.2.60) 


*iAio 


l10 


__*1 

p2  bio  M2 


(*2-i)n, 


(4.2.61) 


L10  +  *10  610  ”  20 


(4.2.62) 
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*1A30  +  Mo  '®10  63q)  “  U  ' 

A30  +  B10  ~  b30  =  0  ’ 
which  may  be  explicitly  solved  to  yield 

A10  =  jr2 <*2+1w  si-l-  +  2^  ■ 

[Mi  \h 

B10  =  b30  =  1  5/+  ^  '2  +  1 

A10  ~  ~B1(T  A10  “  5 

A30  =  0  ’  b10  =  2(A10  ~  • 

In  fact,  the  second  cluster  of  equations  yields 

A50  =  a30  =  B30  =  b50  =  0  ’ 

The  asymptotic  limit  for  is  merely 

»l(0  ~  [a104  +  a10  +  0 

and 

Q{(!)  ~  (A10  -  a1Q)  +  0 


4.2.63) 

(4.2.64) 

(4.2.65) 

(4.2.66) 

(4.2.671 

(4.2.68) 


Equations  (4.2.44),  (4.2.57)  and  the  above  lead  to  the  following  explicit  asymptotic  limits 
for  very  large  elliptic  inhomogeneity: 


Kl"'5^(1+'l2> 

K2'^(1+‘2>/  1+^'2  («-7°) 

Equation  (4.2.69)  is  in  perfect  agreement  with  the  numerical  asymptotic  limit  given  in  Fig. 
4.2.3. 

The  series  solution  provides  us  with  a  complete  family  of  numerically  exact 
bencnmark  results.  Still,  the  resuits  yield  no  useful  analytic  information  concerning  the 
behavior  of  the  solution  for  the  medium  in  the  vicinity  of  the  crack  tips.  This  important 
information  will  be  deduced  from  the  asymptotic  analysis  of  section  4.3. 
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ig-  4.2.1  A  crack  in  a  confocail  elliptic  inhomogeneity 

a)  Configuration  in  the  physical  z-plane 

b)  The  image  of  of  the  arcie  (  =  pQe^is  ^e  e^Pse 

and  the  unit  circle  is  the  crack  in  the  new  complex 
(-plane  via  mapping  function 
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4.3  Asymptotic  Solution  of  a  Crack  in  a  Vanishingly 
Thin  Elliptic  Tnhomogeneity 

4.3.1  Formulation 

Attention  is  now  turned  to  the  specific  situation  where  the  parameter  pQ  used  in 
(4.2.3)  is  almost  equal  to  1,  i.e., 


1 

r  ii 

,  1 

1 

a“2 

■k° 

4- 

o 

'  b  ~  2 

P°  V 

v  J  "  (4.3.1) 

where  a  and  b  are  the  major  and  minor  axes  of  an  ellipse  with  focal  points  located  at 


±1.  The  ellipse  is  thin  if  e  <<  1.  We  have 

a  =  1  +  ^  (e2  -  f3  4-  •  •  • )  i  b  =  e  -  ^  (e2  -  f3  4-  •  •  •)  (4-3-2) 

f(e)  =  b/a«0(e)  (4-3-3) 

where  6  is  another  convenient  parameter.  Moreover,  the  radius  of  curvature  of  the  ellipse 
at  z  =  ±a  is 

Pa  =  X  ~  ^  '  (4-3-4) 

The  ellipse  is  assumed  to  be  the  interface  boundary,  viz., 

C:  z  =  z  =  x  ±  i  ^  (a2  -  x2)1/2  ,  (x2  <  a2)  (4.3.5) 

C  a 

Thus,  as  e  -»  0,  z  has  the  expansion 

zQ  ~  x  ±  i£(l-x2)*/2  ±  i 63  ^  (1-X2)"1/2  +  •  •  •  (4-3-6) 


which  is  valid  for  x2  <  1.  It  follows  that  the  expansion  cannot  be  used  to  satisfy  the  con¬ 
tinuity  conditions  near  z  =  ±a.  The  associated  expansion  is  termed  the  outer  expansion 
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which  will  be  presented  in  Section  4.3.3.  The  terminology  has  its  origin  m  thin  airfoil 
theory  which  is  well-known  (Van  Dyke,  1975). 

To  remedy  the  shortcoming  of  the  outer  expansion,  inner  expansions  will  be  construc¬ 
ted  in  Section  4.3.4.  In  view  of  the  symmetry,  we  shall  concentrate  on  the  region  near  z  = 
a.  The  scaling  factor  for  the  needed  boundary— layer  complex  variable  is  dictated  by 
(4.3.2).  This  complex  variable  is  defined  by  (Fig.  4.3.1b) 

<  =  'TTi  (4-3.7) 


The  appropriate  portion  of  C  is  now  given  by  the  expansion: 


C:  (  ~  (c  ~  i  *  i2 


(1-0 --(1-0  + 

4 


1/2 


where  (  <  1  and  the  leading  term  is  just  the  parabola: 

<c"  (0  =  (W)  +  i2jl  (-»<t ;<+■), 

,-2 


"  [cos  £)  e1^  (-1  <  t  <  ,)  . 


(4.3.3) 


(4.3.9) 

(4.3.10) 


4.3.2  Outer  E^nninn 

The  representation  (4.3.6)  will  be  used  to  satisfy  (2.12)  and  (2.13).  In  view  of  (4.3.3) 
and  (4.3.6),  we  choose  5n  as  the  asymptotic  sequence  and  seek  the  solution  of  a  generic 
unknown  function  Ffz)  in  the  form 

© 

F(z)  ~F(  £  SnFn{z) 

n=0 


(4.3.11) 


where  F  stands  for  anv  of  the  four  unknown  functions  W  and  f  .  The  value  of  the 

a  u 

function  F(z  )  may  be  computed  from  the  scheme: 


F(*c)“  X ff“Fn(*C)  ~  I  ■*"[f’*W*^(l-x:)1/2Fn*(x).+  --] 

(4.3.12) 

where  (4.3.6)  has  been  used  for  z  . 

The  system  of  governing  conditions  (2.7),  (2.12),  (2.13)  and  (2.8)  are 

in  accordance  with  (4.3.1)  and  (4.3.2).  The  J°— terms  are: 

now  expanded 

W*„(x)  +  W'0(x)  -  gx)  =  0 

(4.3.13) 

W^„(x)  +  W'0(x)  -  (J0(x)  =  W^q(x)  +  W*Q(x)  -  (J„(x)  . 

(4.3.14) 

W^q(x)  =  rW^Q(x)  +  7'[w^0(x)  +  W’„(x)  -  f^0(x)]  , 

W2Q(z)  =  Wz  ,  f2o(z)  =  k  as  z  -»  «  , 

(4.3.15) 

(4.3.16) 

where  the  first  three  conditions  apply  to  the  interval  |x|  <  1.  It  follows  from  (4.3.13) 
that  the  left-hand  side  of  (4.3.14)  vanishes  as  well.  It  merely  implies  that  when  an 
interface  is  asymptotically  near  a  traction-free  boundary  it  is  itself  asymptotically 

traction— free.  This  is  essentially  the  nature  of  the  iteration  associated  with  the  outer 

expansion.  At  the  same  time,  it  is  clear  that  the  iterative  mechanism  can  not  be  valid  near 

—1/2 

z  =  ±a  where  the  interface  traction  is  asymptotically  large  in  (z-1)  '  as  z  =  a  -*  1. 

In  any  case,  the  solution  to  the  above  equations  is 

fl0M  =  1f20(z)  =  rf*  • 

(4.3.17) 

2W10W  -  f10w  =  t[2W20(z)  -  f20(«)]  =  T-xX(z)  , 

where 

(4.3.18) 
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The  6— terms  of  the  system  of  equations  are  then  derived.  They  mav  be  further  sim¬ 
plified  by  the  use  of  (4.3.17)  and  (4.3.18).  The  results  are 

Wj^x)  +  W^(x)  -  f^(x)  =  0,  (4.3.20) 

W^x(x)  +  W^(x)  -  f^(x)  =  (7-l)j^(<t-?)x  ±  (f-hT)i(l— x2)1/2  , 

1  (4.3.21) 

Wll(x)  =  T^^x)  +  ±  (£-h?)i(  1— x2)1  /2j  ,  (4.3.22) 

W2l(z)  and  f21(z)  “  0  as  1  m  •  (4.3.23) 

Combining  the  ±  forms  of  (4.3.21),  we  obtain 

f2lW  "  ^l(x)  =  i2(7-l)(f+?)(l-x2)1/2  ,  (4.3.24) 

[2Wn(x)  -  f2  i(x)]  +  [2W2i(x)  “  f2i(x)]  =  2( <r — ^ )( 7 — 1  )x  , 

(4.3.25) 

which,  together  with  (4.3.23),  may  be  solved  for  f21  and  W21  in  terms  of  a  Cauchy  and  a 
Hilbert  problem.  The  solutions  are: 

f21(z)  =  (7-l)(f+I)[X(z)-zj,  (4.3.26) 

2W2i(z)  -  f21(z)  =  -(7-l)(,-j)(X(z)-zj  ,  (4.3.27) 

where  z  is  a  homogeneous  solution  of  (4.3.24)  and  X(z)  a  homogeneous  solution  of  (4.3.25). 
They  are  included  to  meet  the  conditions  at  infinity  Substituting  the  above  into  (4.3.22) 
and  (4.3.20),  we  get 

fn(*)  =  7[(7  -27*)('— *)-(7— l)(f+?)]»,  (4-3.28) 

2Wn(z)  -  fu(z)  =  7l(T-l+27*)(f+f)  -  (7~l)(^)lX(z)  .  (4.3.29) 


The  62— terms  of  the  system  of  governing  conditions  are  obtained  in  a  similar  man¬ 
ner.  They  are  further  simplified  by  the  explicit  lower  order  solutions  and  the  results  are: 

Wl2(x)  +  Wi2W  ~  ^2(x)  =  0  ,  (4.3.30) 

+  w^2(x)  -  fj2(x)  =  (r— i  )[(2T-i)f  -  (27+3)?jx 

*  ( 7-1)1(27-1  )f  +  (27-3)?]i(l-x2)1/2  -  . . ,  (4.3.31) 

±i(l 

W*2(x)  =  7W22(x)  -  77^(f+?)x  *  ( s — <j )i(  1 — x2)1  ^2| 

+  77*  2(7-1  )(f-?)x  ¥  2(7-l)(f+?)i(l-x2)1/2 - ^7»  . 

*i(l— x2)1'2-* 

(4.3.32) 

w22(z)  and  f22(*)-0  as  z  «  «  .  (4.3.33) 

Combining  the  ±  forms  of  (4.3,31)  again  leads  to  a  Cauchy  and  Hilbert  problem,  and  the 
solutions  are 

*22(z)  =  ( 7  1 )[( 2 7 — 1  )f  +  (2 7— 3) ?][X(z) — z]  -  (7-l)?X"l(z)  , 

(4.3.34) 

w22(z)  -  t22(z)  =  (7-l)l(27+3)f  -  ( 2 7 — 1  )f] (X( z ) — z]  +  C(7-l)?X-l(z)  , 

(4.3.35) 

where  the  last  term  of  (4.3.35)  is  a  homogeneous  solution  and  C  an  arbitrary  constant. 
The  other  two  functions  may  again  be  determined  from  (4.3.32)  and  (4.3.30).  They  are 

fl2(z)  =  2|~7(7-1)[3?  -  (27— l)f]  -  77*(f+f)  +  277*(7-l)(f-f)Jz  , 

(4.3.36) 

2W12(z)  ~  ^2(z)  =  [(C-l)(7-l)  -  27*]7ffX"‘(z) 

+  2^7(7  1)[3?  -  (27-I )f  -  77'(,-<01  “  277*(7— l)(f+?)|x(z)  . 


(4.3.37) 


The  procedure  could  be  continued  to  include  as  many  terms  as  we  wish  but  the  expli¬ 
cit  3— term  expansion  is  sufficient  to  reveal  the  most  important  characters  of  the  desired 
solution.  To  unveil  these  properties  and  also  to  facilitate  matching,  we  need  the  inner 
expansion  of  the  outer  expansion.  This  is  accomplished  by  expressing  z  in  terms  of  (  via 
(4.3.7)  and  then  expanding.  The  results  are 

f!  ~  ?f  +  S'  7l(l-7)(f+?)  -  ( I-7-27* )(<r-5))  +  £  ■  +  0(£2)  , 

n=l 


(4.3.38) 


2W1  ~  -  S'  7‘K1/2  +  7?l(7“l)C  -  (7-1+27*)!  -jji 

+  <in  Y,  “jr ' +  a{S1)-  <4-3'39) 

n=2  ' 

<2  «  <  +  i  (l-7)(f+I)  +  (1-7)5  T/I  +  £  +  • 

(4.3.40) 

<D  ,  \ 

2W2  -  ^2  ^  ^  <r(  l  +  (7— l)(<r— +  (7—1)05  ^  •  +  0(0. 

(4.3.41) 

where  the  generic  symbol  (  )n  denotes  an  expression  depending  on  the  explicit  forms  of 
the  outer  expansion,  (4.3.11).  Some  of  them  may  involve  arbitrary  constants  such  as  C 
given  in  (4.3.35).  These  constants  can  only  be  determined  from  matching. 

The  above  expressions  represent  the  properties  of  inner  expansions,  which  will  be 
derived  in  the  following  section,  for  (  -*  ®.  It  follows  from  (4.3.38)  and  (4.3.39)  that 


=  Holomorphic  function  of  ( 


as  4  -  m 


f 


•  1/2 

2W^  —  f^  =  (  '  [Holomorphic  function  of  (] 


(4.3.42) 


a  property  that  can  be  directly  deducted  from  the  traction-free  condition.  On  the  other 
hand,  (4.3.40)  and  (4.3.41)  indicate  that  both  f ^  and  Wj  involve  powers  of  (  ^  as  ( 
-  ®.  This  property  cannot  be  conceived  directly  from  any  of  the  governing  conditions  and 
will  be  needed  in  the  series  representation  to  be  constructed  in  the  sequel. 


4.3.3  Inner  Expansion 

The  independent  variable  is  the  boundary  layer  variable  (  defined  by  (4.3.7). 
Guided  by  (4.3.38)-(4.3.41),  we  shall  seek  the  expansion  of  a  generic  unknown  function 
F(z)  in  the  form  (c.f.(4.3.11)) 


j“F;(0  (4.3.43) 

n=0 

where  F  stands  for  any  one  of  the  four  unknown!)  Wfl  and  ffl.  Using  F*(()  to  denote 
F *((J)  for  simplicity,  we  obtain  from  (2.7),  (2.12)  and  (2.13)  the  governing  system: 

W;*(()  +  W"«)  -  =  0,  ((  <  0)  (4.3.44) 


w;(g  +  w;(y  +  (<c-yw;'(y-f*2(?c) 

=  w*(cc)  +  w;<y  +  «c  -  ?c)wyy  -  fj(?c). 

(4.3.45) 

W2«c)  =  \  Wl<  V  -  f  [Wl«c>  +  wl(?c) 

+  uc-y^X)-fl(y], 


(4.3.46) 
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where  (  is  defined  by  (4.3.8).  As  (  -  »,  the  functions  must  be  matched  with  the  inner 
expansions  of  the  outer  expansions  ^4.3.38)— (4.3.41). 

The  6n— order  conditions  may  be  derived  by  using  the  following  substitution: 

F'ty  ~  X  *np;cy  *  I  '"fay  ’ *!i  ?  d-{)3/2p;-(y  +  •  ■  ■] 

(4.3.47) 


where  (  is  defined  by  (4.3.9).  For  the  6° -terms,  the  choice 


w;0«)  -  jfloto  =  7w;0(o  =  ^0(o  =  l  f 

(4.3.48) 

satisfies  all  the  conditions  and  yields  no  useful  information.  The  if— order  system  governs 

the  desired  asymptotic  limit  (£  -•  0)  and  the  subsequent  discussions  will  be  concentrated 

on  the  solution  of  the  system. 

We  begin  by  writing 

f^U)  =  T [(l-Tl(f-t-I)  -  (l-7-2?’)(«-*)]  +  Ifhj (0  , 

(4.3.49) 

2W*1(0-f^1(0  =  ry/2[i  +  H1(()] , 

(4.3.50) 

&(()  =  (l-7)(f+i)  +  Aj(0  , 

(4.3.51) 

W^(0  =  \  (l-7)[({+!)  -  (»->)]  +  §  (U2  +  #H,(0 

(4.3.52) 

Comparing  the  above  with  (4.3.38)— (4.3.41),  we  conclude  that  the  explicitly  listed  terms 
are  matched  and 


t\,(0  .  Hfl(0  =  0  as  (-*«. 


(4.3.53) 
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Since  the  geometry  is  symmetric  and  the  loading  symmetry  is  maintained  by  the  factor  j 
defined  by  (2.10),  we  conclude  from  the  above  that  the  functions  h^  and  Hfl,  signified  by 
a  generic  symbol  F,  must  satisfy  the  condition 


F(0  =  F(?)  for  hfl  and  . 


(4.3.54) 


Finally,  (4.3.44)  is  identically  satisfied  if  h^  and  are  holomorphic  in  in  the 

C— plane.  We  must  now  determine  h  and  H  to  meet  the  two  remaining  conditions 

na 

(4.3.45)  and  (4.3.46). 

Before  proceeding,  we  establish  a  few  identities  associated  with  the  function  ( Q 
given  by  (4.3.9).  Since 


=  1+iJf , 


we 


may  express  if  in  terms  of 

if  =  +(U  =  2-<V2 


(4.3.55) 


(4.3.56) 


where  4*(()  is  holomorphic  in  the  ((=(^+i(2)-plarre  with  a  cut  along  the  axis,  the 
crack  boundary.  We  have 


*[*!(o]  =  <1/2 . 

1  4 

m  =  ~v*  “ 

1  ^o  "  ^0  1 

*0 


(4.3.57) 

(4.3.58) 

(4.3.59) 


Substituting  (4.3.49),  (4.3.50)  and  (4.3.52)  into  (4.3.46)  and  using  (4.3.54),  (4.3.56) 
and  (4.3.59),  we  find  that  the  condition  is  satisfied  if 


55 


Hj(i)  =  =  \  [<1/2H,(0  +  hj 


+  i^2H]U)  +t(()H 


where 


+  2i(()Hi[i!(0] 


(4.3.60) 

(4.3.61) 


Substituting  the  above  into  (4.3.45)  and  using  (4.3.56)  and  (4.3.57),  we  find  that  the 
condition  is  satisfied  if 


h2«)  =  =  h2«)  =  + 1  (i-7 -?•)  ttj 

+  «C)H1p(o]  +h1(*I(o]  -(T+7‘)  5  -»!«)] 

. 

+  **(l-<1/2]  +  2C1/2H'(0  +  2h'(0  •  .  (4.3.62) 

The  above  two  expressions  H^H^.h^}  and  h^H^h^}  may  be  taken  as  two  calcu¬ 
lating  machines  which  convert  the  input  h^  and  into  the  products  and  Hj  in 
such  a  way  that  the  continuity  conditions  are  identically  satisfied.  For  a  pair  of  chosen 
and  hj  satisfying  (4.3.53),  Hj  and  calculated  from  (4.3.61)  and  (4.3.62)  also  satisfy 
(4.3.53).  The  only  condition  that  is  left  to  be  enforced  is  that  hfl  and  Hfl  must  be  holo- 
morphic  in  in  the  (-piane  (c.f.  ig.  4.3.1b). 
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Let  us  begin  by  setting  H,  =  =  0.  Then, 


h;{o,o>  =  t  JUidl ill 

2  1(77  c1/2[2_41/2]- 


which  have  poles  at  (  =  4,  i.e., 


h;{o,o;  =  p*0(o  = 


((-■») 


hj{o,o} » p;0(o  =  ^  +  ((-4) 


(4.3.63) 

(4.3.64) 


(4.3.65) 

V-x.3.66) 


Thus,  =  h^  =  0  cannot  be  the  correct  choice.  To  obtain  the  complete  solution,  we 
write 


H1(O  =  R10(O 

h«>  =  rio«> 

Hj(0  =  ^{0,0}  +  [Rj0(O  -  P*o«)] 
h2(0  =  hj{0.0}  +  |r20(O  -  P20(()j 


(4.3.67) 


where  RflQ  and  rflQ  are  holomorphic  in  D  .  The  following  senes  are  assumed 


L10 


(0=r  V  —In — 


n=l  4 


■I 


10  L*  (/  x'Q  ' 
n=l  " 

B 


(4.3.68) 


r20^  ~  £  7^7?  ’  R2Q~Yd7^' 


“  C 

n=i  s 


which  conform  with  the  requirement  at  (  =  ®  characterized  by  (4.3.38) — (4.3.41 ).  The 

constants  a  ,  b  ,  A  and  B  must  be  chosen  in  such  a  way  that  the  substitutions  (c.£. 
n  a  n  n 

(4.3.49)— (4.3.52)) 


2Wi  _fi;  7»(1/2R.0iO 

*2:  "[r20^  ~  p20'^ 
W;:»[R2o(()~P;0(O 


(4.3.69) 


satisfy  the  two  continuity  conditions  (4.3.45)  and  (4.3.46).  Since  all  terms  are  bounded, 
the  standard  Fast  Fourier  Transform  algorithm  may  be  applied.  This  provides  a  numerical 
scheme  for  the  direct  computation  of  the  needed  asymptotic  solution  as  opposed  to  the 
numerically  extrapolated  result  of  section  4.2  . 

Let  Kj  and  Kjj  be  the  SlF’s.  They  may  be  normalized  by  the  factors  s ^ and 
al2Ji  ,  i.e., 

K1  =  KI^22^’  ’  K2  =  (4.3.70) 

where  and  K2  depend  on  the  small  parameter  S.  As  6  -*  0,  they  tend  to  finite 
mans  which  may  be  obtained  from  (4.3.50).  The  results  are 


where 


K1  =  7[l  +^(0)' 

M  - +1) 

k2  =  7[i  +  h1(o) 

1  1 

» 

11 

11 

1 

►— * 

*  b 

hi(°)  =  Ri0(°)  = 

y  n 

£,  (H)“ 

(4.3.71) 


(4.3.72) 


The  above  analysis  establishes  the  basis  for  the  general  solution  which  wiil  be  discussed  in 


section  4.4. 


4.3.4  An  Approximate  Analytic  Solution 

It  appears  that  an  approximate  but  explicit  solution  may  be  constructed  for  the 
special  confocal  geometry.  This  is  accomplished  in  the  following  manner.  Another  way  of 
removing  the  poles  defined  by  (4.3.65)  and  (4.3.66)  is  to  define  H1  and  hj  as  follows 


h1(o  =  p;1(o  =  p1(o=— +-^-> 

C-4  «-4)‘ 

hi(0  =  pJi(0  =  Pj(0  =  —  +  — ^-5 

((^)2 


(4.3.73) 

(4.3.74) 


where  P^’ ^12’ P11  anc*  ^12  are  constants  and  both  expressions  are  holomorphic  in  D, 
in  the  (-plane. 

Substituting  the  above  into  (4.3.60)  and  (4.3.62)  and  requiring  that  the  calculated 
and  h2  are  free  of  poles  at  (  =  4,  we  obtain 


and  hence 


where 


8tV(1+H.(0)1 

2Pj(o  +  Pi«)  +  — i  1  ^.0, 

(4.3.75) 

(2— OPjU)  -  P1(()-4[2P'(()  +  Pi(o]  =  0  , 

(4.3.76) 

Pj(0  =  —  f-  —  +  -i-j 

(4.3.77) 

4-»*  l  (-4  «-4)2j 

PjU)-  —  [-  —  4--2-J 

(4.3.78) 

r  =  [1  +  Hj(o>] . 

(4.3.79) 

The  newly  generated  Hj  and  hj  are  now  free  of  poles  at  (  =  4,  but  new  poles  are  gener¬ 
ated  at  (  =  16.  This  is  so  because  H^('l|2(())  now  include  powers  of 


1 
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m-i  (-is 


as  4  -  16 


Removing  the  poles  at  4=16  will  lead  to  new  poles  at  4'  =  36. 

The  process  may  be  continued  for  a  chosen  number  of  times  so  that 


;  4.3.801 


hi  -  ?;nm  - 1  v«  ■ 


(4.3.81) 


where 


hi  =  p!n<0  =  p"(0  • 


(4.3.82) 


p  (4)  =  V  - !H5__ 

m=l  U-(2n)2J 


(4.3.83) 


m=l  LH2n)2J 


(4.3.84) 


which  may  be  expressed  in  terms  of  and  PQ_^(0  ^a  recurrence  formulas.  The 

derivation  is  straightforward  but  lengthy  and  is  not  included  in  the  thesis. 

The  functions  generated  by  substituting  (4.3.81)  and  (4.3.82)  into  (4.3.60)  and 
(4.3.62)  have  poles  at  4'  =  [2(N+l)j2,  i.e., 


h;(p;n«).  Pi*n(o}  =  p;N«) 

h2{PIN^^.  PijjCO}  =  P2n(0 


as  M2(N-fl)r 


(4.3.85) 


where  and  p^  have  forms  similar  to  those  given  by  (4.3.81)  and  (4.3.82).  Thus, 

(4.3.81)  and  (4.3.82),  together  with  the  generated  Hj  and  hj,  cannot  be  the  correct 
solution.  The  complete  solution  must  be  of  the  form 
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H^O  =  PlN^  +  RlN^  1 
^(0  =  pJnU)  +  r:NU)  . 

(4.3.86) 

hj(o  =  h*{p;n,p;n}  +  [r2N(c)-p2N(o], 

hj(0  =  h2{PiN'  p1n)  +  [r2N^  "  P2N^]  ' 

where  R  ^  and  play  the  roles  of  R^q  and  rflQ  in  (4.3.67)  and  must  be  determined 
in  a  similar  manner.  In  fact,  the  series  (4.3.68)  may  be  taken  as  the  solution  for  the  new 
unknowns.  The  input  to  this  problem  is  provided  by  and  p^  which  have  poles  of 

order  2N  at  (  =  (2(N+1)|  .  Noting  that  the  nearest  interface  point  is  at  4  =  1,  we 
conclude  that  the  contributions  of  R^  and  rfl^  are  perhaps  small  for  large  N.  This, 
aowever,  is  not  proven  in  the  thesis.1  Taking  the  statement  for  granted,  we  have 

hx(o«p;n(o,  ^(OsPiVo. 

H2(0  .  H*j{p;n,  P;n}  -  P‘N«)  ,  (4.3.87) 

82(0  =  h2jpiN.  Pin}  -  p2N^)  ■ 

For  a  two-term  approximation,  and  p1  are  given  by  (4.3.77)  and  (4.3.78).  A 
lengthy  computation  yields  the  following  expressions  for  P£  and  Pj 


'The  recurrence  formulas  associated  with  (4.3.83)  and  (4.3.84),  together  with  the  explicit 
forms  of  (4.3.85),  are  needed  for  such  a  proof. 
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p2(0  =  --rp2(0^5(?=JTI  [- — -i-T’" 


4(2  +  <r*)7*l 

1 

1-7’ 

TF^ 

+ 


'64(1— t-t*)  +  48tV*]  1  +  3847*(ll/*-3) _ l_ 


1+1 
\3„* 


(3-16)  7*  1 

1-7’  //_„ 


1-7^ 

i 


((-16)' 


((— 16)v 


((—16) 

where  T  is  defined  by  (4.3.79).  It  follows  from  the  approximation 


HjfO)  «  Pt(0>  +  Pj(«) 


that 


rV* 


Hj(0) 

i+tqw =  Vi-TlW'') ' 

+  2rr-Tl'-J 


_  1  pf2+g-j  ( i-T-r-1 

16— 3<r  |_  7*h7v 


(4.3.89) 


(4.3.90) 


(4.3.91) 


The  approximate  expressions  for  the  asymptotic  values  of  and  (4.3.71)  are 


37(1-7*) 


3-7T* 


57(1-7  *) 


(4.3.92) 


(4.3.93) 
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4.3.5  R-suIts  And  Pisa^on 

We  shall  restrict  our  attention  to  the  range  0  <  <  1  t0  reflect  our  interest  in 

the  physical  situation  of  crack  damage  interaction.  The  confocal  nature  of  the  problem 
allows  a  series  solution  in  a  transformed  plane.  Such  a  series  solution  was  obtained  in 
section  4.2.  The  convergence  of  the  needed  matrix  inversion  becomes  extremely  slow  for  e 
<  0.05.  Nevertheless,  the  e  -*  0  limit  was  accurately  extrapolated  from  the  extensive 
numerical  results.  This  limit  is  reproduced  for  for  the  case  =  0.2  in 

Fig. 4. 3. 2  where  the  numerical  results  generated  by  (4.3.68)  is  also  shown  to  be  in 
agreement  with  the  series  solution.  To  substantiate  the  accuracy  of  our  two— term 
approximate  solution,  (4.3.92)  and  (4.3.93),  Kj  is  also  plotted  in  Fig.4.3.2.  The 
agreement  is  almost  perfect. 

Our  formulation  indicates  that  the  solution  to  the  problem  depends  explicitly  on  the 
two  composite  parameters  7  and  7*  defined  by  (2.14).  It  is  therefore  desirable  to 
present  the  SIF’s  in  terms  of  them.  The  choice  of  the  two  composite  parameters  in  a  typi¬ 
cal  two— phase  problem  is  not  unique  (Dundurs,  1969).  The  parameters  7  and  7*  appear 
naturally  in  our  problem.  Moreover,  exact  solution  may  be  obtained  for  the  case  7*  =  0, 
(4.2.47).  Still,  the  range  of  the  two  parameters  may  be  easily  determined  from  Dundurs1 
discussion.  In  terms  of  7  and  7*,  Dundurs'  a  and  0  parameters  are 


which  may  be  inverted  to  yield 

7  =  iTe 1  7*  =  ■  t4-3-95) 

Dundurs1  a—0  plane,  together  with  the  admissible  ranges  of  a  and  0,  are  reproduced  in 
Fig.  4.3.3.  Also  plotted  in  the  figure  are  the  constant  7  and  7*  lines.  Thus  the  ranges  of 
7  and  7*  may  be  straightforwardly  calculated,  i.e., 
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j  (l- ^7*)  <  7  <  2(1— 27*)  for  0<7’ <5 

0  <  7  <  2(1-27*)  for  \  <  7’  <  5 


(4.3.96) 


which  covers  the  range  of  interest  set  by  0  <  n^ln ^  i  1-  Results  pertaining  to  (4.3.92)  and 
(4.3.93)  are  plotted  in  Fig. 4. 3.4  as  functions  of  7  with  7*  as  a  parameter.  It  is  felt  that 
results  for  all  inhomogeneity  problems  should  be  presented  in  terms  of  7  and  7*  in  the 
fashion  characterized  by  Fig. 4. 3.4. 


* 


an 
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oo  S’  Z 


Fig.  4  3.3 


Admissible  regj 
constants 


4.4  A  Crack  in  a  Thin  Inhnrnop;eneitv  of  Arbitrary  Shai 
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e§ 


4.4.1  .Introduction 

Equiped  with  the  benchmark  results  of  section  4.2  and  the  analytic  information 
revealed  by  the  asymptotic  analysis  of  section  4.3,  we  are  now  ready  to  construct  the 
asymptotic  solution  for  problem  II  where  the  inhomogeneity  is  thin  but  is  otherwise 
arbitrary.  The  geometric  configuration  is  examplified  by  the  right-hand  side  of  (2.16).  It 
is  noted  that  the  nose  dimension  implied  by  the  second  expression  of  (2.16)  is  of  the  order 
of  the  thickness  ,  e  ,  of  the  mhomogeneity.  This  choice  is  made  to  reflect  the  physically 
observed  damage-zone  shape  reported  by  Chudnovskv,  Moet  and  Botsis  (1987). 
Mathematically,  though,  the  situation  is  in  contrast  to  that  of  a  confocal  geometry  where 
the  nose  dimension  is  of  the  order  of  (4.3.2)  and  (4.3.7).  In  this  connection  and  also  in 
terms  of  mathematical  considerations,  the  confocal  situation  is  even  harder  to  solve  than 
the  general  case,  as  the  former  requires  the  satisfying  of  boundary  condition  on  a  parabola 
in  the  boundary  layer  variable. 

Our  final  product  is  an  asymptotically  accurate  series  solution  which  is  valid  in  the 
neighbourd  of  a  crack  tip  for  both  the  inhomogeneity  and  the  medium.  As  a  result,  the 
solution  may  be  manupilated  to  yield  any  and  all  desired  physical  quantities  of  physical 
importance.  For  example,  the  various  generalized  Eshelby  forces  along  the  interface 
boundary  may  be  straightforwardly  computed.  These  quantities  will  be  discussed  in 
Chapter  V. 


69 


4.4.2  Formulation 

We  now  turn  to  the  situation  depicted  in  Fig  1.2.  Geometrically  speaking,  it  is  the 
elasticity  counterpart  of  the  thin  airfoil  problem  (Von  Dyke,  1975).  The  proolem  is  of 
course  much  more  involved  as  it  involves  four  complex  functions.  The  geometric 
description  is  given  by  (2.16)  where,  in  general,  y*(x )  =  y±(x,e)  may  be  assumed  to  have  a 
given  asymptotic  expansion.  We  are  mainly  concerned  with  the  first  order  solution  and 
such  a  generalization  is  not  pursued  in  our  analysis.  The  mathematical  problem  consists  of 
(2.7),  (2.12),  (2.13)  and  (2.8),  and  is  recapitulated  as  follows: 


W^(x,e)  +  W*(x,f )  -  f^(x,e)  =  0  for  |x|  <  1,  (4-4.1) 

W2(zc;e)  +  W2(zc;e)  +  (z~lc)^2  (  *  c  ;  O  -  fj(i c;«) 

=  W^je)  +  W1(zc;e)  +  (zc~ic)W '  (  z  c  ’•  £  )  "  VV'O.  (4A2) 

Wl(V£)  =  7W2(zcie) 

+  7*[  W2(zcif)  +  w2(zcif)  +  (zc~^c^W2  ^  2  c  ‘  e  ^  ~¥2c;£> 

(4.4.3) 

W2(z;e)  =  Wz,  ,  f2(z;e)  =  fz  as  z  -*  ®.  (4.4.4) 


4.4.3  Outer  Expansion, 

We  seek  the  expansion  of  a  generic  unknown  function  F(z)  in  the  form. 

F (•;«)  ~  F0(x)  +  i,(,)Fj(.)  +  <2(<)F2(«)  + .  (4-4'5; 

where  6q(«)  is  an  asymptotic  sequence  that  may  depend  on  zc(x;e).  The  value  of  the 


:o 


function  F(zc;e)  may  be  computed  from  the  scheme 
F(V,<)~E*n(e)Fn!zc) 

~  S  in(E)[F*(x)  .  iey*(x)Fu*(x)  +  ....]  (4-«) 

where  (2.16)  has  been  used  for  a  .  I  n(x)  =  Fq(x  ±  i0). 

The  system  of  equations  (4.4.1)  -  (4.4.4)  are  now  expanded  inaccordance  with  (4.4.5) 
and  (4.4.6).  ^he  0(1)  -  terms  are: 


wio<x>  +  wioW  -  finW  =  0  ■ 

(4-4.7) 

W^0(x)  +  W^x)  -  qQ(x)  =  W^Q(x)  +  W*Q(x)  -  f^0(x)  , 

(4.4.8) 

W^0(x)  =  ,W^(x)  +  T*[w^0(x)  +  W’0(x)  -  f^0(x)]  , 

W2q(z)  =  Wz  ,  f2Q(z)  =  fz  as  z  -  □>  , 

(4.4.9) 

(4.4.10) 

where  |x|  <  1.  As  e  -*  0  the  interface  boundary  is  asymptotically  close  to  the  traction— free 
boundary.  This  is  formally  implied  by  (4.4.7)  and  (4.4.8).  The  vanishing  oi  the  left-hand 

side  of  (4.4.8),  together  with  (4.4.10),  leads  to  the  ordinaarv  crack  solution,  viz. 

f23(x)  =  tl  , 

'4.4.11) 

2W2q(s)  -£2q(z)  =  ffX(z)  , 

(4.4.12) 

where 

x(«)  =  (i.7-i;1/2 

(4.4.13) 

Equations  (4.4.7)  and  (4.4.9)  are  identically  satisfied  by 


(4.4.14  j 


ficAz)  =  T*20^  =  ?fz  ’ 

2W1q(z)  “  f10(z)  =  ?[2W20^  ~  f20^2' 

Unlike  the  situation  described  in  chapter  3.  the  analysis  for  the  present  case  depends 
heavily  on  the  shape  z  (x;e),  (2.16).  For  example,  if  z  is  an  ellipse  then  the 

C  w 

2 

boundary -layer  variable  should  be  scaled  by  e  and  the  required  asymptotic  sequence  is 
just  £n(e)=en.  On  the  other  hand,  if  y*(x)  =  l  in  (2.16)  the  required  asymptotic  sequence 
may  involve  e^ne.  For  this  reason,  we  shall  restnt  our  attention  to  the  determination  of  a 
one— term  approximation  for  a  round-nosed  inhomogeneity  with  (3.6)  as  the  assocciatea 
boundary-layer  complex  variable.. 

To  the  first  order  of  approximation  the  inner  expansion  of  the  outer  expansion  near  z 
=  1  is 


=  'TffX(z) 


(4.4.15) 


f2(l  +  e(;e)  ~  f  +  o(e^2)  . 

(4.4.16) 

2W2(l  +  c(;e )  -  f2(l+«<;e)  ~  e1/2/2  <1/2  +  •  •  • 

(4.4.17) 

yi+efte)  ~  +  o(e1//2)  • 

(4.4.18) 

2W1(l  +  e(;<)  -  f^(l  +  f(ie)  -  &  s  C'2  + 

(4.4.19) 

4.4.4  Inner  Expansion  Near  »  =  1 

For  a  generic  function  F(z;e),  the  associated  inner  expansion  is  defined  by 

f*U;')-f;(0  +  «1/2fJ(0+ . 

where  <^(«)  =  is  determined  from  (4.4.16) — (4.4.19).  Thus 


(4.4.20) 


fI(C;O-t:0(O  +  e1/2fJ1(O+  •••■ 

(4.4.21) 

2W*U;«j -£^((i«)~  faw*  (0-f* 0(0 

+  ‘1/2[2Wj(0-£j1(o]  + . 

(4.4.22) 

%(■•')  ~  $o«)  +  <1/2  ^(0  + . 

(4.4.23) 

+  '1/JPn(0-i51(0|  + . 

(4.4.24) 

In  the  (—plane,  is  infinite  and  is  a  semi-infinite  strip  Fig  1.2.  Equation  (4.4.1) 
requires  f^((;e)  to  be  holomorphic  and  2W*— f^  double  valued  aiong  the  negative  axis.  The 
continuity  conditions  are  just  (4.4.2)  (4.4.3).  Their  conditions  at  (  =  ®  are  governed  by 
(4.4.16)— (4.4.19).  The  choice 

<H)(0  =  7^0(O  =  7*  ■  (4.4.25) 

2W;o(0  -  fj„(0  =  2W*2q(0  -  fJ0(O  =  0,  (4.4.26) 

meets  all  conditions  and  yields  no  useful  information. 

Before  proceeding,  we  give  a  more  precise  descnption  of  the  interface  boundary  z=z c 
in  terms  of  the  boundary-layer  variable  (  defined  by  (3.6).  It  is 

z  -1 

(  =  (c  =  -S - =  pc(f)  e1'  (  x  <  jf  <  x)  (4.4.27) 


and 


Pc  (?)  e1*  -  -  <r  as  f ->  *  X  . 


(4.4.28) 


.Moreover,  the  scaling  of  the  e  is  chosen  in  such  a  way  that 


~3 


Max  /?c(p)cosf  =  1 


(4.4.29) 


It  follows  that 


^c 

T77 


<  1 


(4.4.30) 


Finally,  the  minimum  value  of  p^if)  is  denoted  by  p  ,  i.e., 


Pm  =  <  1  . 


(4.4.31 ' 


The  desired  asymptotic  limit  is  governed  by  f^.,  W^,  ant*  For  the  region 
Dj  in  the  (  -  plane,  the  following  series  are  assumed  : 


*21 


—  =  A  {(/p  ) 
^  n=l 


-n/2 


(4.4.32) 


2Wn(0-^(0 


V? 


=  2(1/2 


1  +  ]?  WJ 


-n/2 


n=l 


(4.4.33) 


where  a*  is  defined  by 

/  =  ^22  for  'l2  ~  °  (4.4.34) 

-1<r12  ffU  =  <r22  =  0 

so  that  ulif*  =1  for  both  cases  (ci.  (2.10)  for  definition  of  a).  It  is  noted  that  the  above 
satisfy  the  conditions  at  (  =  ®  required  by  (4.4.16)  and  (4.4.17).  The  form  of  the  series 
follows  from  the  analytic  structure  revealed  by  the  benchmark  analysis. 

We  employ  the  property  (4.4.30)  to  set  up  the  series  solution  for  region  D^.  It  is 
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(4.4.35) 

(4.4.36) 


so  that  the  conditions  at  (  =  <b  ,  (4.4.18)  and  (4.4.19)  ,  are  identically  satisfied. 

The  constants  involved  in  the  series  are  then  determined  from  the  interface  conditions 
that  have  forms  similar  to  (4.3.45)  and  (4.3.46)  via  the  Fast  Founer  Transform  Algorithm. 
The  associated  stress-intensity  factors  are  : 


K*1)  -  i  K< 1 l>  i  (  K,  -  i  Kn)/.nfi  ~  7  b„  ,  <4-4-38> 

K<2>  +  i  KJ2>  i  (  K,j  +  i  K[) / -  7  b„  ,  <4-4-39> 

where  bg  is  the  same  complex  number  solved  from  the  same  mathematical  equations  which, 
via  the  scaling  factor  a* ,  govern  the  solutions  to  two  different  physical  problems. 


4.4.5  Results  and  Discussion 

The  computer  code  is  used  to  compute  the  K  (®,  b)  curve  discussed  in  section  3.6  and 
presented  in  Fig.  3.2.  Other  relevant  comments  may  be  found  in  section  3.6. 


75 


CHAPTER  V 

ESHELBY  TENSOR  AND  THE  ASSOCIATED  ENERGY  RELEASE  RATES 
5.1  Plane  Elasticity  —  The  r  -  Problem 

We  begin  by  reviewing  the  well-established  results  of  plane  elasticity.  The 
exposition  will  be  delineated  in  such  a  way  that  parallel  development  in  the  manupilation 
of  Eshelby  tensor  in  a  homogeneous  region  will  be  demonstrated.  Most  of  the  results  can 
be  found  in  Eshelby’s  publications  which  will  be  cited  at  the  proper  places.  The  purpose  is 
to  summarize  all  the  known  results  and  the  many  still  not  fully  utilized  concepts  advanced 
by  Eshelby  in  a  unified  setting. 

Let  O  -  (zp  Z2,  z^)  be  rectangular  Cartesian  coordinates  a*  d  let  (i^,  lj,  1,3)  be  the 
associated  unit  vectors.  The  two  dimensional  Kroncker  delta  and  alternator  are  denoted 
by  SaQ  and  e^  A  typical  two  dimensional  region  in  (zp  Z2)  will  be  denoted  by  R,  its 
boundary  by  dR,  and  the  outward  unit  normal  to  dR  by  n.  The  displacement,  strain  and 
stress  are  denoted  by  ufll  and  ,  respectively.  The  strain-displacement  and 
stress— strain  relations  are: 

1/2(ua,0+ufla)  (51) 

\0=dil  77  (5’2) 


where  f  is  the  strain-energy  :ensity  function,  and  n  and  X  the  Lame1  constants. 

We  shall  first  restrict  ourselves  to  the  case  of  constant  /1  and  A  and  consider  the 
r-problem  defined  by  : 

=  0  ,  (5-4) 


7SM 


7St  =  0' 

Tfan0=Ta  ' 


(5-5) 

(5.6) 


2 

where  7“  is  the  Lapiacian  and  are  prescribed  on  dR.  There  are  the  following  associated 


relations  : 


’St  =  0' 

72u  =  -  ( 

a  /i  77  .a 

_  l-f-2u 

U,7~  M  *7/33  eaa,/?’ 


7  «  =  0  , 


where 


"  “  Ca03  Ua,0  ■ 


Equilibrium  considerations  lead  to  the  necessary  conditions 


/ 

*/JR 


T  d<r  =  0  , 

/>  ’ 


(  e  ~  z  T  dff  =  0  , 

J9R  7“3  7  a 


and  the  identity 


I  r  da 

JR  7‘ 


-  f  . 

7 


T  d<r  , 
1 


(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 


(5.12) 

(5.13) 


(5.14) 


which  is  commonly  used  to  define  average  stress  in  terms  of  boundary  traction.  The  trace 
of  the  above  may  be  related  to  the  trace  of 


JR.^da=l/2/jR(V^+uflna)  d„, 


(5.15) 


which  is  commonly  used  to  define  average  strain  in  terms  of  boundary  displacement. 
Multiplying  (5.4)  by  u^  and  integrating  over  R,  we  obtain 


J 

■*  t * 

u  T  d <r  =  r  /jU_  ^da 

7  a  *'R  aw 

(5.16) 

The  trace  of  the  above  is  Clapeyron’s  Theorem  : 

J 

f  T  u  d<r  =  2  1  f  da 

9R  a  a  JR 

(5.17) 

Forming  the  skew-symmetric  part  of  (5.16),  we  obtain 

e/jouTadff—  e^oU  r  da 
/^R  ap  3  a  p  ap  3  a, 7  a0 

(5.18) 

which  must  hold  for  all  R. 

Finally,  in  terms  of  the  complex  variable  z  =  +  izj  and 

W(z)  and  w(z),  there  are  the  following  results  (England,  1971)  : 

two  holomorphic  functions 

2/i  (uj  +  iug)  =  kW  —  zW  '  -  w 

(5.19) 

rH  +  r22  =  2  (W'  4-  VT) 

(5.20) 

^rll  “  r22^  +  l2  r12  =  “2  (zW  '  '  +  w  '  ) 

(5.21) 

(Tx  +  ;.T2)  dir  =  -i  d(  W  +  zW“  +  w  j 

(5.22) 

^z1^2  ~  Z2^V  ^zl^l  z2^2^ 

=  -d  zzW'  -t-  zw  —  J wdz  —  i  Im  j 2 z W'd 1 

(5.23) 

[^2U1  ^lu2^  J<~  *^lul  +  ^2U2^  ^ 

= [n W  -  zW  '  -  w)  d(W  +  z  W'  +  w) 


(5.24) 
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.vhere 


I  (3  —  4i/)  plane  strain 

*■  (3  -  u)  /  (1  +  u)  plane  stress 


(5.25) 


and  v  Poisson’s  ratio.  The  imaginary  part  of  (5.23)  and  real  part  of  (5.24)  are  not 
commonly  used  in  studying  elasticity  problems.  These  well-known  results  are  recorded 
here  because,  as  we  shall  see,  the  Eshelby  tensor  permits  a  representation  that  is  identical 
in  form  to  (5.19)— (5.24).  In  particular,  the  full  form  of  (5.23)  and  (5.24)  are  associated 
with  the  now  well  known  conservation  integrals. 


5.2  The  p  —  Problem  and  r*  —  Problem 

We  shall  discuss  in  this  section  a  mathematical  problem,  the  p-problem,  governing 

the  solution  of  a  tensor  field  p  n  which  is  not  symmetric.  The  p-problem  is  defined  by 

PM  -  0  ,  <««> 


P/M"0' 


'  ea/Bpa<r0' 


(5.27) 

(5.28) 

(5.29) 


p/3a  ap~  Pa  ’ 

where  are  prescribed  on  dR.  The  problem  is  well  posed,  as  there  are  four  unknowns  and 
four  equations.  Indeed,  the  situation  becomes  apparent,  if  we  convert  the  p-problem  into  a 
r* -problem  via  the  transformation: 

\0  =  eaJ3  =e@3  paJ  •  (5'30) 


T*  =  e  ,,  P.  . 
a  a!3  A 


(5.31) 


'9 


we  shall  call  r*^a  pseudo-stress  field. 

Substituting  the  above  into  (5.26)— (5.29),  we  find  that  the  equations  governing 

are  identical  in  form  to  (5.3)— (5.6).  It  follows  that  a  r*— problem  is  a  r— problem.  As  a 

result,  holomorphic  functions  W  (z)  and  w*(z)  may  be  obtained  to  defined  the 

pseudo-stress  field  r*  a  . 

ap 

Since  no  "kinematic  variables"  are  involved  in  the  p— problem,  and  hence  the 
t  —problem,  we  may  even  stretch  the  above  mathematical  analogy  further  by  introducing 


pseudo  strain  rate  e  *a  via  the  inverse  of 
ap 


(5.32) 


where  n  and  1  are  pseudo  viscosity  coefficients.  Since  the  field  is  "Compatible"  by 

(5.30)  and  (5.28),  the  psseudo  strain  rate  e^may  be  integrated  to  yield  a  pseudo  velocity 

field  u*  via 
a 

(5-33) 


To  obtain  the  solution  to  a  mathematical  p— problem  we  merely  express  the 
associated  r  —problem  solution  in  terms  of  W*(z)  and  w*(z),  and  ‘hen  make  :he 
substitutions  : 


012 

~P21 

=  r  *  4-  r  * 
rll  +  r22  ’ 

(5.34) 

(012 

+  p21)  +  i  2  p22  _  (r11  -  r22)  +  i  2 

* 

r21  • 

(pll  _  P22^’ 

(5.35) 

<P1 

+  iP. 

2)  dir  =  i  (Tj  +  i  T*)  d<r  , 

(5.36) 

(*1 

P2- 

z2  Pl)  +  i(zi  Pi+*2  P2)i 

dir 

=  -i 

(Z1  T2  ~  z2  Tl^  +  i  (zi  Ti 

+  z2 

Tl) 

d<r. 

(5.37) 

The  asterisked  versions  of  (5.12)  and  (5.13)  yield 


SO 


!  P  da  =  0  ,  ]  z  P  da  =  0  ,  /c 

JdR  a  JdR  a  a  (5'38) 

and  the  asterisked  version  of  the  trace  of  (5.14)  is 

/}Re«fl3IaP^d'  =  JRe^3!J^da  (5-39) 


The  byproduct  of  the  solution  to  a  r*— problem  is  a  pseudo  velocity  field  u*  which 
may  be  expressed  in  terms  of  W*  and  w*  via  the  assumed  fi*  and  A*.  We  may  therefore 
introduce  yet  another  pseudo  velocity  field  vfl  to  pair  with  the  tensor  .  It  is  defined  as 
follows  : 


v 

a 


Qa0  3  n0  ' 


Ua  =  ea03  V0  • 


(5.40) 


5.3  Eshelbv  Tensor  Associated  with  Ar  —  Problem 

For  the  sake  of  argument  and  convenience,  let  us  assume  that  the  strain— energy 
density  function  f  introduced  in  (5.2)  depends  explicitly  on  place  and  tune,  i.e.,  f  =  f(efl^  ; 
i^\  t).  Using  Eshelby’s  language,  we  may  regard  f(...;  z^,  t)  as  a  calculating  machine  which 
works  ou*  the  value  of  f  when  are  inserted.  The  explicit  dependence  on  zfl  and  t  then 
means  that  the  machine  characteristics  are  continuous  functions  of  place  and  time.  We 
introduce  the  notation: 


1  f 

TT 

1  f 

TT 


1  e  a  and  z  constant 
a/3  a 

d  f 

TT 

7 


(5.41) 


and  zfl(a^7)  constant 


(5.42) 
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so  that 


d  f  _  d  f  •  H  { 
&  t  d  £  Z,(aS  '  d  t 


d  f  d  f 

~Tz  ~1  e 

7 


‘a/«A7  + 


(5.43) 


(5.44) 


Let  be  the  Eshelby  tensor  (1951,  1954,  1956,  1970).  In  terms  of  our  notation,  the 
tensor  is  defined  by 

PA>  =  f  -  rA  (5.45) 


Using  (5.1)— (5.4)  and  (5.44),  we  may  easily  shows  that 
n  -  ^ 

p0a,0--rT  • 


(5.46) 


where  the  right-hand  side  is  zero  if  f  is  homogeneous  in  space.  We  now  show  that  for  such 
a  case  the  Eshelby  tensor  field  is  the  solution  to  a  p  —  problm  defined  by  (5.26)— (5.29). 
Equation  (5.26)  is  identically  satisfied  by  (5.45).  Equation  (5.27)  follows  from  (5.46)  in  a 
homogeneous  region.  Using  (5.1)  and  (5.2),  we  find  from  (5.45)  that 


ea03pa0=-(i+^‘'u7>7  ■  (5-47) 

It  follows  from  (5.7),  (5.9)  and  the  above  that  (5.28)  is  also  satisfied.  The  necessary 
conditions  (5.36)  axe  the  three  well-known  conservation  integrals  (Budiansky  and 
Rice, 1973;  Chadwick,  1975;  Golebiewska— Herrmann.  1981,1982;  Knowles  and  Sternberg, 
..972,  Rice,  1968).  This  completes  the  assertion  that  the  Eshelby— tensor  field  is  the 
solution  to  a  p  —  problem.  Finally,  using  the  identity 


a03  ?a0~  ea03  T 


U  Q 

l3 


=  -  e 


a0  3  r; 


(5.48) 


we  obtain  from  (5.39)  that 
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i*  * 

JSR  +  3  “a  V d'  = 0 


(5.49) 


which  is  another  one  of  the  conservation  integrals2 

We  shall  say  that  a  p  —  problem  and  a  r  —  problem  are  associated  problems  in  a 
homogeneous  region  if  p^  is  the  Eshelby  tensor  associated  with  the  r  —  problem.  We 
recall  that  the  solution  to  a  r  -  problem  may  be  expressed  in  terms  of  (W,  w),  and  the 
solution  to  a  p  —  problem  may  be  obtained  via  the  solutionof  a  r  —  problem  expressed  in 
terms  of  (W*.  w*).  A  tedious  but  straightforward  calculation  leads  to  the  following  simple 


relations  for  a  perfect  associations  : 


(5.50) 


^  (2)  =  2  W'(2)  W'(Z) 


(5.51) 


5.4  Eshdby— Tensor  fores  in  Ela3tostatics 


The  many  important  implications  of  Eshelby  tensor  are  discovered  by  Eshelby  via  the 
consideration  of  the  total  energy  of  a  system.  In  elastostatics,  the  total  energy  EtQt  is  the 
potential  energy.  Using  Eshelby’s  machine  description  f  (efl^  i  z7>  0.  the  total  energy  at 


time  t  is 


=  /  f  da  -  J 


T  u.  da 

m  a  a 


(5.52) 


It  follows  from  the  equilibrium  conditions  (5.2)  and  the  notation  (5.43)  that 

“  JR  -TV  da 


(5.53) 


2The  first  two  equations  of  (5.38)  are  commonly  referred  to  as  are  even  more  commonly 

referred  to  crack-tip  contours  (Rice, 1968).  In  view  of  the  resultant  force  appearance,  it  is 
perhaps  more  appropriate  to  call  them  F  -  integrals.  The  last  of  (5.38)  is  commonly  called 

the  M-integral. 


S3 


which  wiii  be  examined  in  conjunction  with  (5.46),  i.e., 

P/M-rv-  (554) 

The  discussion  cannot  be  continued  without  first  specifying  a  machine  f(  •  •  z^,  t). 

The  machines  that  have  been  considered  by  Esheiby  and  many  other  researchers  are  mostiy 

piecewiseiy  uniform  in  space.  To  fix  ideas,  let  us  introduce  a  curve  or  contour  C  defined  by 

0  :  z  =x  (t)  (5.55) 

a  a 

and  hence  vfl=  xfl(t)  defines  the  velocity  of  C.  If  the  machine  f(  •  •  •;  z^,  t)  is  uniform  in 
space  on  both  sides  of  C,  then  for  a  region  R  containing  C 


E.o,W  =  J. 


da  =  -  !  (f+  -  (~)  v  n  da 
Jr>x  1  a  a 


(5.56) 


where  f*  are  the  uniform  machines  cn  the  ±  n  sides  of  C.  Since  f(  •  •  •;  z  t)  is  a  step 

tsi  U 

function  in  space,  the  right-hand  side  of  (5.56)  is 


-J>+-r)van 


r  u  . 

L -TT\il  ■ 

K  a 


(5.57) 


and  hence 


E.  .(t)  =  -  I  (  P+  -  Prt  )vdf  ,  (5.58) 

totv  '  Jq  v  a  a  '  a 

which  was  first  obtained  by  Esheiby  (1951)  and  is  the  essence  of  the  physical  significance  of 
the  Esheiby  tensor.  If  C  is  the  boundary  5R  of  a  region  R  and  vfl  =  x^  is  linear  in  x^, 
(5.58)  gives  the  energy  release  rates  associated  with  the  "hypothetical",  in  Budiansky  and 
Rice's  words  (1973),  translation,  rotation  and  expansion  of  R.  Esheiby  identified  (5.58) 
with  the  interface  force  associated  with  a  martenstic  transformation. 
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The  various  energy— release  rates  associated  with  (5.58)  play  the  roles  of  generalized 
forces  and  hence  are  of  fundamental  importance  in  analyzing  fatigue  crack-<iamage 
propagation.  If  the  crack-damage  description  could  be  approximated  by  the  two  classes  of 
problems  examined  in  this  thesis,  then  analytic  means  must  be  provided  to  perform  the 
nontrival  calculations  implied  by  (5.58).  Thi.»  was  our  motivation  and  we  have  completed 
our  self-imposed  goal.  The  following  flow  chart  indicates  the  sequence  of  steps  that  may 
be  straightforwardly  performed  to  complete  (5.58)  • 


Series 

W(  z  ) 

via  (5.50) 

i  K (•) 

* 

solution 

w  (  z  ) 

rmr! 

w  (  z  ) 

T  a  0 

(5.34)  | 

P 

!  Energy 

(5.37)'”  | 

|  ^  a  0 

a 

1  1  |  Release  Rates 

>  -  - - 

The  linking  of  the  energy-release  rates  to  the  damage  evolution  kinetics  is  the  intent  of  our 
next  objective. 


CHAPTER  VI 


SUMMARY 

Our  interest  in  toughness  of  composites  reinforced  by  particles,  fibers  and  layers,  and 
generalized  forces  associated  with  crack-damage  propagations  have  led  us  to  consider  the 
tw<?  classes  of  elasticity  problems  stipulated  in  Chapter  I.  Two  computer  codes,  one  for 
each  of  the  two  cases,  have  been  developed  and  are  ready  for  adaptation.  The 
mathematical  complication  associated  with  the  small  geometric  parameter  has  been 
effectively  removed.  The  remaining  arbitrariness  associated  with  a  general  physical 
problem  may  be  handled  by  one  or  several  of  the  conventional  analytic  and/or  numerical 
methods. 

The  computer  codes  effectively  give  the  coefficients  of  a  series  solution  which  is  valid 
in  the  neighbourhood  of  a  crack  tip.  It  follows  that  they  may  be  manupilated  to  yield  any 
and  all  physical  quantities  of  significance.  In  particular  the  various  Eshelby  forces  along 
the  interface  boundary  may  be  easily  computed  by  following  the  flow  chart  provided  in 
Chapter  V. 

The  efficiency  of  the  two  codes  are  illustrated  and  the  class  of  problems  depicted  in 
Fig.  3.1  is  discussed  in  section  3.6.  Since  we  are  most  interested  in  the  situation  where  the 
inhomogeneity  is  softer  than  the  medium,  Mode-I  SIF’s  for  five  cases  are  summarized  in 
Fig.  6.1,  where  the  normalized  SIF  is  plotted  against  ^or  v\  =  v2  ~  ^  *s  Seen 

that  the  effect  of  an  inhomogeneity  is  strictly  of  a  shielding  nature  for  the  range  ^  <  1- 
Moreover,  the  exact  value  of  the  SIF  for  a  specific  configuration  falls  in  the  rather  narrow 
lens— shaped  region  bounded  by  the  very  large  inhomogeneity  limit,  (4.2.69)^  and  the 
straight  line  K  =  which  is  the  exact  SIF  associated  with  a  semi— infinite  crack 

perpendicular  to  a  two— phase  boundary  (Hilton  and  Sih,  1971).  We  conclude  the  summary 
with  the  following  rather  unexpected  conjecture. 

A  conjecture:  Let  a  crack  be  surrounded  by  a  doubly  symmetric  inhomogeneity  of 


moduli  and  =  v  which,  in  turn,  is  embedded  in  an  infinite  medium  of  moduli  ^  and 
i/j  =  v.  The  inhomogeneity  may  consist  of  two  disjointed  inhomogeneities  surrounding  the 
tips.  The  normalized  Mode— I  SIF  is  a  function  of  the  two  composite  parameters 


(1+*2^1  _  **1 

*  _  1  ^1  X  _1  X 

^  1+*^  ^  1+*  /*2 

and  is  denoted  byK^yj*).  It  satisfies 


7  <  (7.7*)  < 


7  (7+I-7*) 
2(7+7*)  (1-27*) 


for  0  <  7  =  ^1/^2  ~  **  nght-hand  side  of  the  above  is  deduced  from  (4.2.69)  and 
equalities  hold  for  7  =  =  0  and  1. 

The  curves  summarized  in  fig.  6.1  apparently  satisfy  the  conjecture.  Curves  (l)  and 
(5)  are  respectively  the  right-hand  and  left-hand  sides  of  (6.1).  Curve  (4)  is  the  thin 
confocal  ellipse  and  (3)  the  circular  inhomogeneity.  Curve  (2)  corresponds  to  the 
configuration  of  Fig.  3.1  with  a  =  0.5  and  b  =  2.0. 
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ABSTRACT 

Cohesive  elasticity  is  the  grade— 3  theory  of  elasticity  developed  by  Mindlin  in  1965. 
has  a  modulus  of  cohesion  which  gives  rise  to  surface-tension.  The  concept  of  adhesion 
is  introduced,  and  interfacial  energies  and  energy  of  adhesion  are  defined.  The  interfacial 
energy  solution  may  also  be  used  to  define  a  grain  boundary  energy.  Also  presented  are  the 
thin  film  energy  and  the  concept  of  an  interface-phase.  The  stretching  of  a  thin  film  is 
analyzed  in  deta.1;  and  it  is  found  that  the  apparent  Young’s  modulus  obtained  from  a  film 
is  higher  than  that  obtained  from  a  plate. 
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1.  Introduction 

Every  continuum  field  theory  contains  only  a  limited  amount  of  physics,  and  even 
the  limited  amount  of  physics  is  very  often  only  of  a  phenomenological  nature  in  the  sense 
that  the  various  constants  employed  in  the  theory  are  closely  related  to  but  not  directly 
derived  from  more  fundamental  physical  parameters.  It  could  be  argued  that  Yeung’s 
modulus  should  be  directly  derived  from  more  fundamental  bonding  forces,  but  the 
"correct"  Young’s  modulus  has  always  been  directly  measured  from  a  tension  test.  It  is 
hard  to  define  the  exact  physics  content  of  the  theory  of  elasticity,  but  it  is  perhaps 
reasonable  to  say  that  most  of  the  so-called  micromechanics  phenomena  we  are  so  eager  to 
understand  are  not  derivable  from  the  limited  amount  of  physics  that  was  built  in 
elasticity.  When  the  implications  of  a  continuum  theory  have  been  exhaustively  revealed 
by  the  mathematic  solutions  to  the  associated  field  equations,  additional  mathematical 
manipulations  become  redundant  even  though  new  solutions  will  always  be  of  usefulness. 

In  terms  of  a  lattice  theory,  elasticity  incorporates  only  the  nearest  neighbor 
interactions,  and  that  is  it.  The  theory  does  not  have  an  intrinsic  length  scale  and,  as  a 
result,  a  30— cm  slab  behaves  the  same  as  a  40— /an  film,  and  there  is  no  difference  between 
a  microcrack  and  a  geological  fault.  An  intrinsic  length  scale  appears  when  the  forces 
between  particles  are  extended  to  include  first,  second  and  n— th  neighbor  interactions. 
Moreover,  an  initial,  homogeneous,  self— equilibrating  stress  will  lead  to  surface  tension 
(Toupin  and  Gazis,  1964;  Gazis  and  Willis,  1964).  The  surface-tension  solution  is  of  a 
boundary— layer  type,  and  the  associated  decay  constant  has  been  estimated  from 
electron-diffraction  data  obtained  by  Gerner,  MacRae  and  Hartman  (1961). 

The  continuum  version  of  the  n— th  neighbor-interaction  lattice  theory  is  the 
so-called  grade— n  theory.  If  the  strain  energy— density  is  assumed  to  depend  on  the 
rotation— gradient,  in  addition  to  the  strain,  there  results  the  couple  stress  theory.  A 
complete  grade— 2  theory  depends  on  all  eighteen  components  of  the  strain-gradient,  and 
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the  additional  stress  quantity  is  sometimes  termed  the  double  stress  (Toupin,  1962).  The 
inclusion  of  yet  a  third  gradient,  which  has  thirty  independent  components,  leads  to  a 
grade— 3  theory  and  the  associated  new  stress  quantity  is  sometimes  termed  the  tnple  stress 
(Mindlin,  1965).  There  is  also  a  very  general  theory  which  includes  strain-gradients  of  any 
order  (Green  and  Rivlin,  1964). 

Couple  stress  and  higher  order  gradient  theories  were  popular  topics  of  research  in 
the  sixties,  and  Mindlin’s  work  are  most  noteworthy  in  that  his  goal  was  specifically 
targeted  at  understanding  the  effect  of  microstructures  on  the  failure  of  solids  (Mindlin, 
1962,  1964,  1965a,  1965b,  1968;  Mindlin  and  Eshel,  1968).  In  particular,  and  to  the  best  of 
our  knowledge,  his  grade— 3  theory  (Mindlin,  1965)  is  the  only  continuum  theory  that  is 
fully  developed  and  has  the  capability  of  defining  surface  phenomena  via  nontnvial 
displacement  fields.  The  most  important  concept  of  this  theory  is  a  new  constant  called 
the  modulus  of  cohesion  which  is  essentially  an  initial,  homogeneous,  self— equilibrating 
triple  stress.  It  is  for  this  most  important  constant  that  we  have  coined  the  name  cohesive 
elasticity  to  stress  the  significance  of  the  theory. 

The  availability  of  self— equilibrating  states  is  fully  exploited  in  this  paper  In 
addition  to  surface  free  energy,  which  was  considered  by  Mindlin,  we  have  introduced  the 
concept  of  an  adhesive  joint  to  define  various  interfacial  energies  and  energy  of  adhesion. 
When  suitably  defined,  an  interfacial  energy  becomes  the  grain  boundary  energy.  Also 
introduced  are  the  thin  film  energy  and  the  concept  of  an  interface-phase  The  stretching 
of  a  thin  film  is  analyzed  at  the  end.  It  is  shown  that  the  apparent  Young’s  modulus 
obtained  from  a  film  is  higher  than  that  obtained  from  a  slab. 

The  basic  equations  are  first  recapitulated  in  a  dimensionless  form.  Stress  function 
representations  are  then  presented  in  terms  of  solutions  to  thirteen  uncoupled  second  order 
equations.  A  complete  boundary-layer  solution  is  obtained  for  a  regular  boundary,  and 
several  significant,  outstanding  problems  are  outlined  at  the  end. 
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2.  Notation 

We  consider  vectors  in  three-dimensional  Euclidean  space.  Such  vectors  wiil  be 
denoted  by  lower-case  bold-face  letters  (a,  b,  n,  etc.).  Let  e2<  e^  denote  the  unit 
vectors  in  the  directions  of  the  three  coordinate  axes  z^,  Zj,  z^  of  a  rectangular 
right-handed  cartesian  system  with  origin  0.  The  position  vector  of  a  point  P  relative  to 
0  is  denoted  by  a. 

In  addition  to  scalar  product  a«b  and  vector  product  axb  of  the  two  vectors  a  and 
b,  the  dyadic  product  is  written  a«b  which,  in  case  of  no  confusion,  is  also  written  ab.  A 
second-order  cartesian  tensor  T  can  be  expressed  as  a  linear  combination  of  unit  dyads, 
so  that  it  takes  the  form 

T  =  T  e®e.  (2.1 ) 

ij  i  J 

where  the  summation  convention  is  employed.  As  a  rule,  we  shall  use  bold-face  capitals  to 
denote  cartesian  tensors  of  second  (and  higher)  order.  Moreover,  an  over-scnpt  is  used  to 
denote  the  order.  Thus,  a  cartesian  tensor  of  order  n  can  be  expressed  in  components  as 


n 


n  indi  c  es  n  factors 


n  m  n  m 

The  inner  product  A«B  and  outer  product  A®B  are  illustrated  by 


A*B  =  A-jBj^e-»e,K  , 

(23) 

2  2 

(2.4) 

A«B  =  A.jBy  e-«ej«e^»e^  . 

n  m 

In  addition,  the  rule  of  contraction  A*  B  of  tensors  is  illustrated  by 


(a1a2)*(blb2b3)  =  (vMa2#lVb3  ’ 

(ala2a3)*(blb2)  =  ai(a2«bi)(a3*b2)  , 

(ala2a3)*(bib2b3)  =  (a1.b1)(a2*b2)(a3.b3)  , 

where,  for  example,  (a^a2)  =  a^®a2.  The  dyadic  product  indicator  ®  will,  for 
convenience,  not  be  explicitly  shown  in  places  where  no  apparent  confusion  will  be  resulted 
from  the  omission. 

The  theory  developed  by  Mindlin  involves  a  material  length-scale  L  which  is 
assumed  to  be  much  smaller  than  the  geometnc  length-scale  L  associated  with  a  problem. 
The  ratio  f  =  l/L  plays  a  key  role  in  the  theory,  as  well  as  in  the  ensuing  asymptotic 
analyses  and  results.  It  is  for  this  reason  that  we  have  chosen  to  recapitulate  the  theory  in 
a  dimensionless  form. 

Throughout  this  paper  the  length-scale  L  is  prefixed  so  that,  e.g., 

Dim.  *  =  L  *, 

and  hence  z-  are  dimensionless.  Let  n  be  the  shear  modulus.  The  force-scale  used 
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throughout  this  paper  is  just  L  fi.  Also,  the  notation  Dim.(  )  is  used  to  mean  dimensional 
form  of  (  ),  instead  of  dimension  of  (  ). 

3.  Governing  Equations 

Let  u(*)  be  the  dimensionless  displacement  vector.  Following  the  original  work  of 
Mindlin  (1965),  the  dimensionless  strain  energy-density,  W,  is  assumed  to  be  a  function 
of  three  displacement-gradient  tensors: 

W  =  W(E,E,E)  (3.1  M 

where 

E  =  *  (Vu  +  uV),  Eij  =  \  (u;  a  +  Uj -)  ,  (3  2) 


(2.5  : 

(2.6  i 
(2.71 


3*Multiply  W  by  (i  to  recover  its  dimension. 
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E  =  ?Vu 


E  =  m  u 


Eijk  Uk,ij  ’ 
Eikjl  =  ul,ijk  ' 


(3.3) 
(3.4). 4 


and  V  is  the  gradient  operator  in  the  dimensionless  b.  The  symmetry  properties  of  the 
above  are  seif— explanatory. 

The  variation  of  the  total  strain  energy,  in  a  (dimensionless)  volume  V,  with 
arbitrary  variation  of  u,  is 

f  C  2  3  4 

Wdi/=J  [T*(V^u)  +  T.(VVJu)  +  T*[W6v)\dv  (3.5) 

V  V 

where  dv  =  dz^dz^dz^  and 

T  =  dW/dE  ,  T  =  dW/dE  ,  T  =  dVJ/dE  .  (3. 6). 3 


The  right-hand  side  of  (3.5)  may  be  reduced  to  include  a  surface  integral  by  application  of 
the  chain  rule  of  differentiation  and  the  divergence  theorem.  In  doing  so,  however, 
additional  caution  must  be  taken  to  recognize  the  fact  that  Vtfu,  together  with  other 
similar  terms,  is  not  independent  of  S\x  on  the  (dimensionless)  surface  <?V  of  the  volume 
V. 


A  number  of  surface  operations  and  surface  operators  are  needed  in  completing  the 
desired  reduction.  Let  a  be  the  unit  normal  to  dV  and  pointing  a  .vav  from  V.  The 
following  are  applicable  and  defined  on  dV : 


V£u  =  ad  <5u  +  V°£u 
n 

(3.7) 

d  e  n.V 
n 

(3.8) 

V°  e  (I  —  nn)«V 

(3.9) 

$  1  _  n-i_2 

4  Multiply  E  by  L  L  to  recover  its  dimension. 

*  n  n-2 

5  Multiply  T  by  jiL  to  recover  its  dimension. 


~y 


7’  e  mVu«n)  -  V°  = 


I  +I 

R1  R2 


n-  r 


(3.10) 


where  I  is  the  unit  tensor  and  (1/R^  4-  I/R2)  the  mean  curvature  of  the  surface.  It  is 
noted  that  for  a  flat  surface  V’  is  just  the  negative  of  V°,  the  plane  gradient  operator. 
Using  the  above  relations,  we  obtain  from  (3.5) 

Wdi/  =  -J  jv*  [t- V»T  +  V*  [v®t]  ]  |  •  Svdv  +  J  [t»*u  +  t*dn<5u 


V 


dV 


2  2  -i 

+  t»d  £u|  da 
n  J 


(3.1i; 


where 


t  =  n  •  [t  -  V.T  +  (  W)*t]  +  7'.  [v’»  (h»t)]  +  V’.jn*  iT-V*T 

-(?0n).[(nn).T]}  , 

t  =  (iui).[t-V«t]  +  n*[?*.(n.T)]  +  V’* [(nn)*T]  , 


t  =  (mm)*f  , 


(3.12) 

(3.13) 
(3.14)*S 6 


are  the  generalized  surface  traction  vectors  associated  with  the  generalized  surface 
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kinematic  vectors  u,  dnu  and  dnu,  respectively.  The  Generalized  Principal  of 
Stationary  Potential  Energy  and  stress— equation  of  equilibrium  are: 


S  J  Wdi/  =  J  ft^udi'  +  J  f  t»<hi+t«dn<hi+t*dn<hi  da  (3.15 

V  V  dV 


V# 


2  3  4  1 

T-V*T+(VV)*T 


+  f  =  0 


(3.16) 


*  n  n 

6  Multiply  t  by  /zL  to  recover  its  dimension. 


For  homogeneous  and  isotropic  materials,  the  following  2nd  degree  polynomial  W 
was  deduced  by  Mindlin: 

W=  [e..E..  +  E-.E.-l  +  e2/3E.r. 

ij  ij  2  ( /c  1 )  u  jjJ  Ha  njj 

+  f^aiEijjEikk+a2EiikEkj,+Q3EiikEjjk  +  a^E^  +  aE^E^ 

+  +  /34EiijkEllkj+/35EiijkE]lik+  tfE^E^ 

+  ^EijklEjkli]  +  e2[^EUEiikk^2EijEijkk+^EijEkkij]  (3-17)’7 8 


where  e  =  if L, 

’  3 — 4i/  (general  and  plane  strain) 

K  =  •  8 

3— v  (plane  stress  of  elasticity  )* 

T+iJ 


(3.18) 


and  v  Poisson’s  ratio.  All  the  remaining  constants  are  expressed  in  dimensionless  form 
via  n  and  a  material  length-scale  t,  which  is  assumed  to  be  small  so  that  e<<l.  We 
shall  refer  to  e  as  the  small  parameter  from  time  to  time.  For  example,  as  the  small 
parameter  tends  to  zero  (3.17)  tends  to  the  strain  energy— density  of  the  theory  of  linear 
elasticity.  All  important  mathematical  and  physical  implications  of  the  assumed  W  may 
be  found  in  Mindlin’s  original  work  and  are  amplified  by  the  explicit  introduction  of  L 
and  e.  Finally,  the  material  constants  introduced  by  Mindlin  are 

V  b0'  bn'  Cn}  =  M 4’  ^7n}  ' 


7  The  dimensional  forms  of  the  various  terms  are  illustrated  by,  e.g.,  t26  E-... 


o  njj 


* 

8  The  meaning  of  plane  stress  in  the  new  theory  should  be  reexamined. 
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From  (3.6)  and  (3.17)  follow  the  constitutive  equations: 

Tpq  =  [2Epq  +  fel  Eiiipq]  +  (![nEiij/p,  +  1'2Epq,i 

+  5^%,+%)  ■ 

Tpqr  ~  ^  “t  (Epii^qr+Eqii  V  +  5a2:  Eiipl5qr+2Em^qp+Enq'ipr^ 

+  2°,Biir*pq  +  2*(Epqr+°(  %+Erpq)  , 


(3.19) 


(3.20) 


Pqrs  M3  ^iEujjWs  +  H^ku^jkpqrs  8  ^[[Eiijk+Eiikj]  ^jkpqrs  +  2Ejsii^jpqr] 
+  3  ^Eiisj^jpqr  +  3  ^Eiijs^jpqr  +  2^6Epqrs  +  3  ^Eqrsp  +  Erspq  +  Espqr^} 

+  e"  Vpqrs+  H  Vijpqrs+  37jBtS  W]  +  l^pqrs  1  ^ 


(3.21) 


wheere  6-^  is  the  Kronecker  delta  and 


*ijkl  "  Vkl  +  ^ik^jl  +  ^il^jk  - 

^ijklmn  ^ik^jl^mn  +  ^ik^jm^ln  +  ^il^jm^kn 


(3.22) 

(3.23) 


kk 

The  symmetry  of  T  follows  from  that  of  E,  and  the  very  last  term  of  (3.21)  was 
identified  by  Mindlin  as  a  cohesive  force  which  gives  rise  to  surface  tension.  As  we  shall 

see  that  it  also  gives  rise  to  interface  phase  and  interphase  interface  energies. 

* 

Let  W  [u]  be  the  total  elastic  energy.  Then  (3.17)— (3.21)  lead  to  the  conclusion 


1  =  J  Wdv  =  j  e20^  J  V2V*u  d^  +  j  J  (t*E  +  T*E  +  T*E]di/  (3.24) 


that 
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where  the  appearance  of  the  first  integral  is  a  result  of  the  fact  that  the  0  —term  in  W  is 

o 

linear  in  E-.-.  The  second  volume  integral  may  be  converted  into  a  surface  integral  by 
following  the  steps  from  (3.6)— (3.11),  and  the  net  outcome  is  the  Generalized  Clapeyron’s 
Theorem: 

2  J  Wd a  —  J  f«uda  +  e2/3^  J  V2V*udi/  4-  J  [  t*u  +  t«dnu  +  t*d  u]da  .  (3.25) 


Finally,  the  total  potential  energy  U  of  the  system  is 


*  *  i'*  i'*  ro  1  2  2  1 

U  =  W  —  I  f*udi/ —  I  I  Uu+Ud  u+t*d  u  da 
•>  J  n  n  J 


(3.26) 


0  1  2 


It  follows  from  the  above  that  if  there  are  no  external  forces,  i.e.  if  t,  t,  t  and  f 
are  identically  zero,  the  total  potential  energy  is  merely 

U  =  W*  =  \  e20a  J  V2V*u^eW  =  j  e20o  J  dnV* Je^da  (3.27) 

°  v  0  dv  n 

where  u  =  u  ,  if  exists,  is  a  self— equilibrating  state.  The  total  potential  energy  of  the 

system,  (3.26),  is  the  free  energy,  either  the  Gibbs  free  enerrgy  or  the  Helmholtz  free 

energy.  When  a  traction— free  surface  is  the  source  of  the  disturbance,  the  free  energy  is 

(3.27).  It  follows  that  the  dimensionless  surface  free  energy  per  unit  surface  area  F  is 

s 

given  by 


Fs  =  \  *l0  O.J 

s  £  on 


(3.28) 


the  dimensional  form  of  which  is  termed  the  surface  tension  by  Mindlin.  We  shall  see  that 
(3.26)  can  be  used  to  define  many  surface  related  energies. 

The  displacement -equation  of  equilibrium  deduced  from  (3. 16)— (3 .2 1 )  is 

K.+  1  2  2  2  2 

D  VV*u  -  D  D  VxVxu  +  f=  0 

K— 1  tl  12  21  22 


(3.29) 
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where 


■  *  =1'2'  '•-1-* 


(3.30) 


2  2  1/2 
2a  i  =  a’-7  *  [(a-?,)  -4/3’l 

2A  3  3 


(3.31) 


(3.32) 


in  which  A=l,2;  and 


a  =  2(q  +q  +a  +a  -fa  )  ,  a1  =  2(a  -fa  )  , 

V  1  2  3  4  S'  3  4 

0  =  2(0+0+0+0+0+0+0j)  ,  0'  =  2(0 +06) 

7=  7  +72+73  • 


(3.33) 


As  the  small  parameter  e  tends  to  zero,  (3.29)  reduces  to  the  displacement— equation  of 
equilibrium  of  the  theory  of  elasticity. 

As  it  was  pointed  out  by  Mindlin  that  the  conditions  for  positive  W  do  not  include 

relations  between  a(or  a')  and  0,  7  (or  0',!^)  and  hence,  supply  no  indications  of  the 

character,  real  or  complex,  of  the  four  constants  a In  all  the  specific  solutions 

considered  in  the  sequel,  the  a ^  will  be  treated  as  if  they  were  real  and  positive;  but 

complex  a ^  are  equally  admissible.  The  character,  real  or  complex,  of  the  four  constants 

a dictates  the  behavior  of  the  field  variables.  On  the  other  hand,  for  the  most  part,  the 

final  products  of  the  problems  to  be  analyzed  are  energy  expressions  of  one  form  or  another. 

These  energies  turn  out  to  be  indpendent  of  the  character  of  a In  any  case  the  two 

pairs  of  constants  (a  ,a  )  and  (a  ,a  )  are  assumed  to  be  positive  throughout  this  paper. 

11  12'  21  22/ 

In  case  they  are  actually  complex,  the  following  substitution  is  presumed: 

a  =  a  ,  Real  a  >  0  ;  (3.34) 

12  11  u 

a  =  a  ,  Real  a  >  0  ;  (3.35) 

22  21  21 
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which,  of  course,  satisfy  (3.31)  and  (3.32). 

Equations  (3.5),  (3.11),  together  with  the  equilibrium  equation  (3.16),  yield  the 
identity 


J  [T*(Vfa)  +  T*(Wfu)  +  T*(m«u)]di/ 

V 

J  ft&idi'  +  J  [t«<hi+t*<Mu+t«dn£u]da  ,  (3.36) 

V  av 


which  may  be  used  to  establish  certain  useful  integral  identities.  The  result  of  setting 
£u=c,  an  arbitrary  constant  vector,  is 


c« 


f  f° 

J  fdi/  +  J  tda 

V  av 


=  0, 


(3.37) 


and  the  vanishing  of  the  resultant  force  is  a  necessary  condition  for  (3.16)  to  have  a 
solution.  The  result  of  setting  <fai=<i/xs,  where  u  is  an  arbitrary  constant  vector,  is 


u  • 


r  r  °  1 

J  8  x  fdi/  +  J  (*  x  t  +  n  x  t)da 


dV 


=  0, 


(3.38) 


and,  hence,  the  vanishing  of  the  moment  becomes  another  necessary  condition.  The  rest  of 
12  1  2 

t  and  t,  i.e.  t«n  and  the  full  t,  are  self-equilibrating  and,  hence,  may  be  arbitrarily 
specified. 

It  is  possible  to  obtain  additional  integral  indentities  from  (3.36)  which  may  become 
useful  in  establishing  the  expected  new  Eshelby  (1951,  56,  70)  tensors  and  the  associated 
conservation  laws  (Knowles  and  Sternberg,  1972;  Budiansky  and  Rice,  1973),  but  the 
possibility  is  not  pursued  further  in  this  paper. 

The  stress  equation  of  equilibrium,  (3.16),  may  be  written 

7*Tc+f  =  0  (3.39) 


where 


Tc  a  T-V.T  +  V.(».T) 


(3,40) 
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may  be  considered  as  the  Cauchy  stress  tensor  with  the  associated  Cauchy  stress  vector  tc 
defined  by 


tC  =  n.Tc  =  t-r  •{...} 


(3-41) 


where  V’*{...}  may  be  identified  with  the  terms  given  in  (3.12).  It  is  clear  that 

j  fdi/  +  J  tcda  =  0  ,  (3.42) 


dV 


J  8  x  fdi/  +  J  s  x  tcda  =  0  . 


(3.43) 


V 


dV 


The  following  surface  integral  indentity  will  enable  us  to  show  that  (3.37),  (3.38)  are 
equivalent  to  the  above. 

Let  A  be  an  aribtrary  second  order  tensor.  We  have 

J  V°.Ada  =  J  V°«|n(n.A)  +  [A-n(n.A)]J  da 
dv  dv 

=  J  (V°.n)(n.A)da  +  V°.[A-n(n.A))  da  (3.44) 

dV  dv 

where  V°  is  the  plane  gradient  operator  defined  by  (3.9).  For  a  smooth  closed  surface  dV, 
the  last  surface  integral  vanishes.  It  follows  that 


J  [(V°«n)(n.A)  -  ?°.A)da  =  f  V’»Ada  =  0, 
dV  dV 


(3.45) 


which  is  the  desired  identity. 

In  the  usual  continuum  theory  of  thermodynamics,  the  work  term  is  axiomatically 
defined  to  be 

J  n*Tc*uda  (3.46) 

av 

even  if  Tc  is  assumed  to  depend  on  7u,  ?7u  and  VVVu.  This  led  Dunn  and  Sernn  (1984) 
to  conclude  that  "a  troubling  aspect  of  all  higher-grade  models  is  that  they  are  in  gneeral 
incompatible  with  the  usual  continuum  theory  of  thermodynamics".  They  then  proceeded 
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to  introduce  the  so-called  interstitial  working  into  the  work  integral.  Mindlin’s  derivation, 
however,  clearly  indicates  that  the  correct  work  term  is  the  surface  integral  of  (3.26),  and 
am  appropriately  modified  continuum  theory  would  lead  to  thermodynamic  compatibility. 
Similar  discrepancies  also  show  up  in  other  higher  order  gradient  theories  (see,  e.g.  Gurtin, 
1989). 


4.  Stress  Functions 

The  solution  to  the  displacement— equation  of  equilibrium  was  obtained  by  Mindlin 
in  terms  of  the  solutions  to  four  uncoupled  partial  differential  equations  of  order  six.  We 
give  a  more  direct  deduction  and  show  that  Mindlin’s  representation  can  be  further 
simplified.  Moreover,  our  deduction  is  constructive  in  that  it  lends  itself  to  the  actual 
construction  of  general  solutions.  The  well-known  Galerkin— Somigliana  and 
Papkovich— Neuber  representations  serve  as  the  starting  point,  and  the 
component -function  approach  employed  by  Doyle  (1967)  is  then  used  to  reduce  the  system 
to  a  set  of  thirteen  uncoupled  2nd  order  equations. 

The  displacement-equation  of  equilibrium  (3.29)  is  first  rewritten  in  the  form 


d2  D2  -  D2  D2  ]  ?V.u  +  D2  D2  V2u 

ll  12  21  22J  21  22 

In  case  of  zero  body  force,  the  divergence  of  the  above  yields 

D2  D2  V2V*u  =  0  . 

11  12 


+  f=0  . 


(4.1) 


(4.2) 


Applying  D^D^V  to  (4.1)  and  using  (4.2),  we  conclude  that  for  zero  body  force 


D2  D2  D2  D2  !74u  =  0  (4.3) 

11  12  21  22 

which,  for  the  case  e  =  0,  is  just  the  well-known  result  that  u  is  biharmonic.  We  are 
thus  led  to  the  Generalized  Galerkin— Somigliana  Representation: 


2 

K+ 1 


U 


D2  D2  V2G  -  Id2  D2  -£=l  D2  D2 

11  12  L  tl  12  *+1  21  22. 


VV*G  . 


(4.4) 
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The  result  of  substituting  (4.4)  into  (4.1)  is 


D2  D2  D2  D2  V4G  =  -4rf 
11  12  21  22  *+l 

which,  together  with  (4.4),  completes  the  desired  generalization. 
Let  the  first  term  of  (4.4)  be  denoted  by  B,  i.e. 


(4.5) 


B  =  D2  D2  V2G  .  (4.6) 

11  12 

Then,  by  (4.5), 


D2  D2  V2B  =  --Lf  .  (4.7) 

21  22  *+l  V  ’ 

If  the  objective  is  to  reexpress  (4.4)  in  terms  of  B,  the  solution  to  a  lower  order  equation, 

the  remaining  terms  in  (4.4)  must  be  appropriately  modifed. 

We  begin  by  expressing  the  solution  to  (4.7)  in  terms  of  three  component  functions 

as  folllows: 


B  =  bM  +  f2C 

B^  +  e2C 

B« 

(48) 

21  22 

D2  D2  B  ,  bW= 

D2  72B  ,  B^2 

)=  D2  V!B  , 

(4.9) 

22  21 

22 

21 

V2b(°) 

II 

1 

* 

+  t-J 
>-* 

(4.10) 

D2  B^ 
21 

-  2  f 
*+ 1 

(4.11) 

D2  B^ 
22 

=  2  f 
K+l 

(4.12) 

where 


=— A-  ■  c,b  =  4^t-  .  (4.13) 

a02-a0l  a02-a0i 

in  which  0  =  1,2;  and  C^,  C  are  defined  for  a  similar  equation  to  be  developed  in  the 
sequel.  The  validity  of  the  above  decomposition  may  be  easily  verified  by  direct 
substitution  (Doyle,  1967). 
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The  divergence  of  (4.6)  is 

D2  D2  V2V*G  =  V*B  (4.14) 

11  12 

which,  together  with  the  first  of  (4.9),  yields 

D2  D2  V2(D2  D2  V*G)  =  .  (4.15) 

11  12  21  22 

From  the  identities,  which  are  applicable  for  any  F, 

V2mF  =  *«V2F  +  27«F  ,  (4.16) 

D^*p  =  **Di(3f'-2<JvJ*p  ■  <4'I7> 

The  following  representation  is  conceived: 

D2  D2  V.G  =  l  a.B^  +  i  B  (4.18) 

21  22  *  ^  0 

where  B  is  a  new  scalar  function.  Substituting  the  above  info  (4.15),  we  find  that  B 
o  0 

must  satisfy 

D2  D2  V2B  =  ~  [d2  D2  s.f-  2e2(a2  +a2  D2  )V«f 
11  12  0  *+l  11  12  11  12  11 

=  "XT  I0'  o'  [■•*-  2e2(a2  +  a2  )V-fl  -  2eV  C  D2  V2V.f 
*+l  I  11  12  L  11  12  11  11  12 


—  2e4a2  C  D2V2V»f  .  (4.19) 

12  12  11  J 

It  follows  from  the  above  that  the  component  representation  for  B^  is: 

B  =  +  e2C  B^+  e2C  B^  (4.20) 

00  11  0  12  0 

B(o>  _  D2  D2  B  ,  =  D2  V2B  ,  B^  =  D2  V2B  ,  (4.21) 

0  11  12  0  0  12  O  0  11  0 

V2B^*  =  2  [*.f-  2e2(a  +  a  )V.fi  ,  (4.22) 

0  *+l  (.  11  12 

D*B(‘)=__4  2a2V#f  (4.23) 

11  0  11 

D2  B^  =  --Ke2a2  ' 

12  o  K+i  12 


We  have  thus  effectively  expressed  the  last  term  of  (4.4)  in  terms  of  and  B  . 


(4.24) 
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Substituting  (4.8)  and  (4.9)  into  (4.14),  we  get 


V2 

D2  D2  V«G  -  t2 

C  D2  +  C  D2 

V«B 

11  12 

.  21  22  22  21. 

. 

=  V*B 


(o) 


which,  in  view  of  (4.16),  permits  the  convenient  representation 

D2  D2  V*G  =  e2(C  D2  +  C  D2  )?*B  +  i  a»D'  '  -t-  a  d  - 
11  12  21  22  22  21  ° 


-a.B^°)  +  i  b(°) 


where  B^  satisfies 


V2b(°)  =  -a.V2B^0^  =-|Ta«f  . 

0  K+l 


(4.25) 


(4.26) 


(4.27) 


The  second  term  on  the  right-hand  side  of  (4.4)  is  now  expressed  in  terms  of  B  and  B 


(«) 


by  (4.26).  In  view  of  (4.22)  and  (4.27),  it  is  convenient  to  decompose  Bj  1  into  two 
parts  as  follows: 

b(°)*  =  d>  +  B^,  (4.28) 

o  o 

V2*  =  -  4  f2(n2  +q2  )v.£  (4.29) 

Y  K+l  X  11  12' 

The  final  form  of  the  Generalized  Papkovich-Neuber  Representation  is  obtained  by 
substituting  (4.6),  (4.18)  and  (4.26)  into  (4.4),  viz. 

>(.)  *  1 


2u  =  (ft-f-l)B  -  Vb«B'0'  +  5  V 


(«-i)b>i)bW 


-  £2(«.-i-1)(C  d!  +  c  d!  }VV*B 
21  22  22  21 


(4.30) 


where  the  last  term  may  be  further  simplified  by  (4.8)  —  (4.13).  The  most  explicit  form  of 
the  representation  is 

2u  =  jj«+l)B^  -  Va*B^  -  VB^ 

e2(«+l)  fc  B^  +  C  B^-(a2  +  a2  )VV*b(o) 

v  L  21  22  21  22 

+  J  (St]»  +  c„b<,,!  +  c,,b(,j)]] 

-  e4(*+l)VV«|c  a2  +  C  a2  B^l 

v  '  [  21  21  22  22  j 


(4.31) 
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where  e  2<p  is  the  particular  solution  of  (4.29),  as  the  homogeneous  solution,  which  is 

harmonic,  may  be  absorbed  by  B^.  The  elasticity  solution  is  recovered  as  e  tends  to 

o 

zero.  We  note  that  each  one  of  the  component  functions  is  governed  by  a  second  order 

equation.  The  functions  B^,  B^  and  tfi  are  associated  with  the  Laplacian  operator 

0 

and,  hence,  has  no  boundary-layer  phenomena.  The  functions  B^,  (a  =  1,2)  are 
governed  by  the  operators  defined  by  (3.30)  and  contribute  to  the  new 

boundary— layer  phenomena. 


For  antiplane  deformations  defined  by  u^(z  ,  z^)  we  have 

2u  =  (*+1)bM  +  e2(*+l)(C  B^  +  C  B^) 

3  3  21  2  22  3 


(4.32) 


where  B'a'  are  of  the  boundary-layer  type.  For  plane  deformations  defined  by  ua(z^, 


z  )  the  displacements  are 


2u  =  kB^-B^  -z„B^  +  e2(*+l)  [c  B^  +  C  B^  -  (cr  +q2  )B 
a  a  o,o  p  P,o  21  o  22  o  22  P.P01 

1  u-11  r  (i)  (*)  11 

+  ijvr  C  B  +  C  B 

L  (_k+1J  ^  11  0,0  12  o,Q  JJ 

4  ,  _(0  1 

-  e  (/i+l)  C  a2  B a  o  +  C  a  .5  o  n  (^ 

L  21  21  Apo  22  2.  ApoJ 


(4.33) 


(A)  (A) 


where  B  B  are  of  the  boundary-layer  type.  It  follows  that  the  terms  B  are 
ct  0  o>u 

expected  to  have  the  most  significant  contribution  to  the  new  boundary-layer  phenomena. 

Thus,  the  constants  a  and  a  are  also  expected  to  be  more  significant  than  a  and 

11  12  21 

a  .  We  shall  see  that  the  former  pair  appears  in  all  the  surface  energetic  expressions. 

22 


5.  Boundary  Layer  Along  A  Smooth  Boundary 

The  appearance  of  the  small  parameter  e  in  the  operator  defined  by  (3.30) 

indicates  that  the  main  difference  between  elasticity  and  cohesive  elasticity  lies  in  a 
boundary  layer.  An  elasticity  solution  is  but  the  interior  solution  of  an  associated  cohesive 
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elasticity  problem.  The  self— equilibrating  state  u  =  given  in  (3.27)  is  a 

boundary— layer  phenomenon  and  was  discussed  by  Mindlin. 

To  fix  ideas,  let  us  consider  a  body  with  a  smooth  boundary  defined  by  s  =  a0.  The 

following  traction  boundary-value  problem  is  considered  for  the  case  f  =  0: 

°o  , 

t  =  r(«°)  (5-1) 

t  =  eT^a^e),  (5-2) 

t  =  e2r(a0;e),  (5.3) 

where  the  functions  on  the  nght-hand  side  are  prescribed  generalized  tractions,  which  may 
be  specified  to  depend  on  e.  The  explicit  appearance  of  e  in  the  governing  equation 
(3.29)  and  surface— traction  expressions  (3.12)  -  (3.14)  suggests  the  appropriateness  of  a 
regular  perturbation  solution  in  the  form 

u(a,e)  u(°)(a)  +  eu^(a)  +  ...  .  (5  4) 

Since  (3.29)  contains  only  powers  of  e2,  it  is  clear  that  both  u^  and  satisfy  the 

ordinary  elasticity  equation,  i.e. 

[Elasticity  Operator]  and  u^  =  0  ,  (5-5) 

which  are  of  second  order.  The  solutions  to  (5.5)  cannot  be  adjusted  to  satisfy  the  three 
sets  of  boundary  conditions  (5.1)  —  (5.3).  The  expansion  (5.4)  is  thus  termed  an 
outer-expansion  (Cole,  1968,  and  Van  Dyke,  1975).  The  appropriate  boundary  conditions 
for  (5.5)  can  only  be  determined  from  matching.  To  facilitate  such  a  computation  it  is 
necessary  to  set  up  a  set  of  normal  coordinates  relative  to  the  boundary  surface  s  = 
Such  coordinate  systems  are  routinely  used  in  the  theory  of  shells. 

Let  (  (a=l,2)  denote  a  set  of  Gaussian  coordinates,  so  that  the  boundary  surface 

is  given  by 

*  =  ••  =  »•«„)  =  ■“((,.(,)  •  <5-6! 
The  associated  base  vectors  and  first  fundamental  tensor  are  denoted  by  a^  and  a^. 
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In  addition,  a  third  base  vector  a  is  identified  with  the  normal  n  to  the  surface,  so  that 


the  position  vector  2  becomes 


=  «0(O  +  (a.  ,  a  =  n 


;5.7) 


where  (  is  the  dimensionless  distance  from  the  surface  along  n.  The  interior  of  the  body 

is  therefore  defined  by  i  <0. 

3 

In  terms  of  the  normal  coordinates  a.,  the  displacement  vector  is 


u(a;e)  =  u  ({)  =  u^O  \l J  .  (5-8) 

* 

so  that  u-  are  the  physical  components.  The  outer-expansion  (5.4)  is  now  rewritten  as 


u  ~ 


+  <  '■‘W«a.o)<J  +  ••• 


a./v  a-- 

y  jj 


(5.9) 


which,  for  small  values  of  (  ,  may  be  approximated  by 


u 


uj'(')((Q,0)  +  u‘W(fa,0)(]+ 


a./V  a.- 
y  JJ 


(5.10) 


The  coordinate  (  is  now  stretched  to  become  a  boundary-layer  coordinate  (  defined  by 


(5.11) 


The  inner— expansion  of  the  outer-expansion  is  obtained  by  substituting  (5.11 )  into  (5.10), 
and  the  result  is 


u”  lu_/°^Q>0)+e 


u*(l)Ua-°)+K!  3)((q-0) 


+ 


a  ■ 
JJ 


(5.12) 


Let  t-^  ((j)>  n  =  0,1,  denote  the  physical  components  of  the  ordinary  elasticity 

*  *(n) 

stress  tensors  associated  with  the  ordinary  elasticity  solution  u- v  '((■),  i-e- 


*  *  '?_«•  * 
r.  ■  =  2e. .  +  — - y  6-  -e,  , 
ij  lj  ii-1  ij  kk 


(5.13) 


where  e-j  are  the  physical  components  of  the  ordinary  strain  tensor.  On  the  surface  ~ 
0,  we  have 
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e*  =u*  ((  ,0) 
3  3  3  ,3^Q  ' 


e  = 
aa 


—  ua,a^0’°'>  +  ^~u (a  =  l>  2>  no  sum) 


(5.14) 


aa 


aa 


*  _  1 
£  —  a 


a% 


aa  aa 


where  R  ,  R  are  the  normal  radii  of  curvature.  The  curvatures  are  positive  if  the 
u  22 

centers  of  curvature  are  on  tne  side  f  <0. 

3 

The  desired  inner-expansion,  consistent  with  (5.12),  is 

*(o) 


u  ~ 


Uj  fUa,0)  +  «Vj((a,0  +  ... 


(5.15) 


where  v-  are  again  physical  components.  The  governing  equations  deduced  from  (3.29) 


are 


1-a2  Jl 

fi  2^1 

1  -  a  — 

u  del 

12  dfi 

d(2 


3  _ 


=  0  , 


(5.16) 


1  —  a2 
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”  H2 J 


i  d2' 


dC 


=  0  . 


The  leading  terms  of  the  displacement  gradient  tensors  are: 


E  =  v  , 
33  3 .{ 


E*  -1 

3a  2 


a,{ 


+  2«*<.)  *<•> 


Ct3 


a, 3  J 


*  *(•)  - 
a0  ea0  ' 


E*  •  =  -  v-  „  ; 
331  e  i,U 


*  1 

E  =  —  v-  ,  ( t  ; 

3331  (2  l,((( 


(5.17) 


(5.18) 


where  all  components  are  physical  and  refer  to  the  base  a-. 

There  are  more  than  a  few  nonvanishing  stress  components  associated  with  the 
above  but  only  the  following  are  needed  for  the  boundary  conditions: 
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T*  -  *+1  V  4-  -YV  1  3~*  *U) 

33  K—  1  3>^  +  ^3  ,U(  +  /C-l  aa 

T*  =  T*  =v  >+2v  ttt  +  2e*W-u*(o)  ; 
3a  03  a,(  2  a3  al3 


T  =  eav  ,,  , 

333  3,44 


(5.19) 


T  =  c2 
3333 


’jj/3=S2(a  +a)vflu  ; 

(fl.  +  VA.<«>  +  7.  'o 


(.)' 


333^ 


=  e 


)  73  Va,(  +  2<^  +  *Vva,«(]  +  7)(2e«  ’  -  ua!’M 


Finally,  the  leading  terms  of  the  traction  vectors  are  obtained  from  (3.12)  —  (3.14).  They 
are 


0 

t  = 


*  1  *  1  * 

rp  A  rp  [  x  m 

3l  £  23X,(  f2  3331,  it) 


(5.20) 


\lf*u  ’ 

t  =  [t*  ■  -  -  T*  .  ,1  a .//7.  , 

l  331  e  3331, fj  1'  u 

1  /J~ai  -  . 

' y  u 

The  traction  boundary  conditions  (5.1)  -  (5.3)  are  first  rewritten  in  terms  of  (  and 


2  * 

t  =  T  .a/ 
3331 


a-,  i.e. 


*  =  Wf  J*  4  =  t  =  f2T*(f  ;e)  . 


'a"  "a 

The  final  form  of  the  boundary  conditions  are 


0>  +  T,i,(  +  ^3,(U  =  ljlUa^-\  ‘aa] 


(5.21) 

(5.22) 

(5.23) 


S=T \<  +  {2T-a)\,U(+l*UUU  =  **«.>  ■  Ela  )((a'°);  I5  24) 

*  ^  + 2<  Wa,<« = -  v2^/’  -  <:’>  ■  <5  25> 
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[a(Oj  +  a)-^7j]vAi({-2Ws  +  (3()vJjW{  =  {-*(^;0)  ,  (5.26) 

VA,( +  h  - 2(1 Va.K(« +  2(/3s +  VyKUU( 


=  ^v+“2IV0)-2^°)(^'0) 


(5.27) 


The  following  identities,  which  are  obtained  from  (3.31),  will  be  used  repeatedly: 


a2  +  a2  =  ^tt(<i-27)  ,  a2  a2  =  0  i 

tl  12  *+lv  "  11  12  *+l 

a2  +  a2  =  2(a  +  a  )  —  7  ,  a2  a2  =  2(0  +  0  )  . 

21  22  3  *'  3  21  22  v  5  6J 


(5.28) 


Using  these  identities,  we  find  that  (5.24),  (5.27)  are,  respectively,  the  first  integrals  of  the 
governing  equations  (5.16),  (5.17).  The  solutions  to  these  equations,  for  (  <  0,  are 


*(  \  *(  ^  C/“  { 

v  =uUJ(£  +  ,0)  +  A  e  11  +  B  e 

3  3  a  3.3  a  13  13 


*0=^0  «a'°>  +  fuA3  (fa'°>  +  k20e  n  +  B20e 

*(o)  0*  *(o)  *(  \ 

UA3^a>0)  =  ^^a)  +  uA3(^*0)-2ei0^a>0)  > 


(5.25) 


(5.26) 


(5.27) 


(5.28) 


where  the  last  two  conditions  are  the  results  of  matching,  i.e.  the  outer— expansion  of  the 
inner— expansion  must  match  the  inner-expansion  of  the  outer  expansion  given  by  (5.12). 
Using  (5.13)  and  (5.15),  we  find  from  the  above  that  the  boundary  conditions  for  the 
elasticity  solution  u-'i;(^)  are  just  the  traction  conditions: 

=  °’«a)  ,  (5  29) 

and,  hence,  the  first  term  of  the  outer-expansion  (5.9)  is  completely  defined. 

The  remaining  boundary  conditions  (5.22),  (5.23)  and  (5.25),  (5.26)  involve  only 
derivatives  of  v  and,  hence,  the  constant  terms  u- ^\^Q,0)  cannot  be  determined  from 
the  inner  expansion.  A  lengthy  calculation  yields  the  following  concise  results  for  the  other 
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two  sets  of  constants: 


a/3j- 


f  2  2-j2 

+  da  )  da  L-  ~ 


Ct  /v  +  da  d  a  L.  -  On^  +  CLa  da  L- 

0° _ Pi  Pi  h  PQ _ Pit  th  Ji 


2  2 

"b  Cla  I  d  n  L; 


*0r 


a&>  +  “a 


2  2  1 
+  a 


% 


[*00  +  4«]2  “ft 


(5.30) 


a0>  +  afl|  Q&  L\i  +  \a0 


2  2 

+  da 

?0  ft 


2  1 


laA>  +  aA 


’A 


iA> 


2  2  -)  2 

+  “ftj  “ft 


(5.31) 


where 


2  1  2  1 

a  =  ~ TT  7  ,  cr  =  *  7  . 

10  K+l  20  2  '3 

(5.32) 

L  =  a  a  T*  (/  ;0)  , 

31  11  12  ft+1  3  v  Q 

(5.33) 

L  =-*-1/3 

32  *+l  0 

-a2  u*^°\{  ,0)  +  ~4[t*(£  ,0)  -  7 

10  3,3  /t-fl  3  '  1  dd 

(5.34) 

la  =  • 

(5.35) 

L/3,  =  -  “j.  "ftf  <V>  +  -  TJ(2‘ft  ’  -  - 

(5.36) 

exception  of  the 

*f  ) 

terms  u-^1'  ((  ,0),  which  are  but  constants 

in  the 

inner-expansion  (5.15),  the  desired  two-term  inner-expansion,  (5.15),  is  now  complete. 

Had  we  actually  began  with  the  presumption  that  are  complex,  (3.34)  and 

(3.35),  then  Bp  j  should  be  replaced  by  j  and  (5.30)  and  (5.31)  still  hold,  as  the 
latter  is  merely  the  complex  conjugate  of  the  former. 


6.  Traction— Free  Surface  Free  Engery 

Consider  now  a  body,  of  volume  V  and  bounding  surface  dV,  that  is  completely 
free  of  external  forces,  i.e. 
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(6.1) 


It  follows  from  (5.229)  that  the  first  term  of  the  outer— expansion,  u*^(s),  is  identically 
zero  and,  hence, 


by  (5.33)  -  (5.36). 


where 


L»— HFT*.  ■L3i  =  LA  =  Lft  =  0' 

(6.2) 

The  inner— expansion  (5.5)  is  merely 

u  =  u^  ~  ev  ({)  a  //"a"  i 

3  3  33 

(6.3) 

£/<*  U* 

v=Ae  11  B  e  12  , 

3  13  13 

(6.4) 

A  =  *“!  a2  (a2  +a2)/?/D  , 

13  *+l  nv  10  i2y  V  1 

(6.5) 

B  =£=ia2(a2  +a2)/3/D  , 

13  *+l  12'  10  11  O'  1 

(6.6) 

D  =  (a2  +  a2  )  a  -  (a2  +  a2  )  a  , 

1  10  12  11  10  11  12 

(6.7) 

which  are  determined  by  (5.30),  (5.31). 

The  above  solution  is  a  self-equilibrating  state  and,  hence,  (3.27)  and  (3.28)  apply. 


Thus, 


U’  =  U’  =  AF  ,  F  =  i  e  0  (a2  -  a2  )/D 
e  s  *  s  2  *+1  'n  v  i2/# 


'e  *“  s  ’  *  s  2  ~  K+i  "o '  ii  12"  i 


(6.8) 


where  A  is  the  total  surface  area  and  F  the  free  surface  energy  associated  with  the 

s 

traction— free  surface  of  a  completely  unloaded  body.  The  dimensional  form  of  Fg  is 

(k-1)  b2  (a2  -  a2  ) 

Dim.  F  =- -  0  -H - 12- 

s  2(k+1)h13  D 


(/t-l)b2  a  +  a 

= - 2 - ii - 11 -  ,  (6.9) 

2(k+1  )fil3  (a2  -  a  a  ) 2  -  a  a  (a  +  a  )2 

V  10  11  12'  11  12V  11  12 

which  was  first  obtained  by  Mindln.  Thus,  to  the  first  order  of  magnitude,  as  e-»0,  the 


26 


surface  energy  associated  with  a  regular  traction— free  boundary  of  a  completely  unloaded 
body  is  a  material  constant. 

The  expression  (6.8)  remains  the  same  if  a  ^  is  the  complex  conjugate  of  a^, 
(3.34).  Since  F5  must  be  positive,  we  have 

(6.10) 


2  2 
(a2  -  a  a  )  -  a  a  (a  -ha  )  >0, 

10  11  12  11  12  11  12 


a  condition  that  must  be  satisfied  by  the  new  material  constants. 

The  surface  energy  given  by  (6.8)  is  the  energy,  per  unit  area,  associated  with  the 

formation  of  a  new  regular  surface.  The  free  end  of  a  thin  wire  is  not  a  regular  surface  as  a 

whole  unless  it  is  sufficiently  thick,  in  terms  of  the  material  scale  L ,  as  it  is  interacting 

with  the  free  cylindrical  surface.  As  to  how  thin  is  thin  or  how  thick  is  thick  the 

self— equilibrating  state  associated  with  a  free  cylinder  must  by  analyzed.  It  goes  without 

saying  that  the  surface  in  the  neighborhood  of  a  notch  (or  crack)  tip  is  not  regular  in  the 

sense  that  (6.8)  is  also  not  applicable.  The  conclusion  is  that  the  energy  required  in  the 

formation  of  new  surfaces  is  not  always  determined  by  (6.8),  and  many  relevant 

self— equilibrating  states  must  be  analyzed  to  understand  the  phenomenon.  In  short,  the 

energy  required  to  scoop  a  smooth  marble  out  of  a  chunk  of  solid,  or  slice  it  into  halfs,  may 

be  computed  from  (6.8).  The  so-called  toughness  may  not  be  directly  obtained  from  the 

material  constant  F  . 

s 

Let  us  now  consider  a  traction-free  surface  on  a  loaded  body.  To  avoid  the 
unnecessary  complication  of  dividing  up  the  bounding  surface  into  different  portions,  we 
assume  that  the  load  is  merely  a  body  force  so  that  the  bounding  surface  is  still  completely 
free.  The  outer— expansion  u*,0)(£),  (5.9),  is  now  no  longer  zero.  Using  the  fact  that  dV 
is  traction  free,  we  have 

<6n> 

where  T*^0) (( a>0)  is  the  surface  (elasticity)  stress  invariant  produced  by  the  nonzero 
body  force  f.  The  inner— expansion  associated  with  the  still  traction-free  dV  is  now  given 
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by  (5.15),  (5.25),  (5.26),  (5.30),  (5.31)  with 


L  =0,  L  =  +  a  JAr*'0|((rt 

31  32  /t+1  ^0  10  2(7— tl)  QCt  'a 


«„.°) ; 


=0'  ' 


(6.12) 


(6.13) 


The  total  potential  energy  of  the  system  may  be  computed  from  (3.26),  (3.25),  viz. 


U*  =  -  j  J  f«udi/  +  j  e2PQ  3n?«uda  , 


(6.14) 


as  the  surface  dV  is  still  free.  For  e  =  0,  we  recover  the  ordinary  elasticity  potential 


energy  associated  with  the  system,  i.e. 


u’L=  =~Uf,u 


<°>dt/ 


(6.15) 


(6.17) 


where  u,0) (■)  =  u*(0) ((). 

The  second  integral  of  (6.14)  is  now  computed  to  give 

U‘  =  UE1  +  AfVO  (616) 

where,  to  the  first  order  of  magnitude, 

FS(Taa>  =  2  e4 

=5''3op^-“.5r^T:r«a'°)]  ■  (6i7) 

which  is  to  be  compared  with  (6.8).  It  is  noted  that  the  contributions  of  u(0)  and  vq  in 

the  surface  integral  of  (6.14)  are  of  one  order  lower  than  that  of  v^  and,  hence,  do  not 

appear  in  (6.17).  The  contribution  of  the  modulus  of  cohesion  0  is  a  constant  triple 

stress,  the  last  term  of  (3.21),  which  gives  rise  to  F  (0)  =  F  of  (6.8).  It  follows  that  the 

s  s 

interaction  between  d  and  f,  which  gives  rise  to  r* 1 0) ,  is 

o  aa 

)  -  F  =  -  r'  << "  F,  (6.18) 


s  s'  aa 


K)P0  aa  s 


(6.18) 


which  is  negative  if  the  surface  stress  (not  to  be  confused  with  surface  tension,  a  term  that 
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is  deliberately  avoided  in  our  presentation)  invariant  T*ala0)  is  positive.  The  total 
potential  energy  of  the  system  now  becomes 

U*  =  U*x  +  AFs  +  A(AFS)  .  (6.19) 

Consider  now  an  infinite  slab  of  thickness  H  with  a  cylindrical  hole  of  radius  R. 
The  slab  is  uniformly  stretched  at  infinity  by  a  tensile  stress  t  =  aef  while  er  is  the  unit 
radial  vector.  The  elasticity  potential  energy  is 

UE1  =  U*-i^lo2R2H  (6.20) 


where  U*  is  the  energy  associated  with  the  slab  without  the  hole  and  the  second  term  is 

a> 

the  so-called  flaw  (hole)  energy  (Sih  and  Liebowitz,  1967).  The  cohesive  elasticity 
potential  energy  may  be  obtained  from  (6.16)  and  (6.17)  with 

=  2a  +  2i/a  =  a  (plane  strain)  .  (6.21)*9 

Thus, 


U*  =  U*1  +  2,RH  F., 


W, 


(6.23) 


and  the  change  of  U*  associated  with  a  hypothetical  hole  expansion  is 

5U* 


^  = -ii*+ilERo2  +  2xh  [i-(^h  J  Fs 

[no 


(6.24) 


where  the  first  term  is  the  elasticity  potential  energy  release  rate.  The  vanishing  of  (6.24) 
is  the  condition  for  the  hypothetical  expansion  and  the  elasticity  counterpart  of  such  a 
condition  is  the  Griffith  criterion: 


-  ll*+H  HRo2  +  2tHFs  =  0  .  (6.25) 

The  relevancy  or  irrelevancy  of  the  hypothetical  expansion  condition  is  not  the  purpose  of 
our  discussion.  Rather,  we  use  (6.24),  (6.25)  to  illustrate  the  main  difference  between 
cohesive  elasticity  and  elast  ;y.  The  surface  energy  term  is  a  part  of  the  totcil  potential 
dergy  computed  from  the  fie.u  solution  in  (6.24),  whereas  it  is  merely  an  appended  side 


♦ 

9  The  stress  concentration  factor  is  2. 
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condition  in  (6.25).  Moreover,  the  role  of  the  interaction  energy,  (6.18),  is  apparent  in 
(6.24). 


7.  Interfacial  Energies 

Let  us  now  assume  that  the  upper  and  lower  half  spaces  V±  are  taken  from  two 
different  solids.  We  propose  to  glue  them  together,  by  a  layer  of  adhesive  of  vanishing 
thickness,  to  form  a  composite  body  V.  Such  a  process  cannot  be  nontrivially  defined  by 
the  ordinary  theory  of  elasticity.  Noting  that  the  most  striking  difference  between 
elasticity  and  cohesive  elasticity  is  the  existence  of  a  constant  triple  stress  =  e2/3Q, 
(3.21),  we  assume  that  the  desired  adhesion  is  the  result  of  a  sigular  variation  in  (3q 
Moreover,  the  two  half  spaces  V*  are  assumed  to  be  identical  in  every  respect  other  than 
the  possibility  of  a  difference  in  0Q.  Before  proceeding,  it  is  convenient  to  summarize  the 
relevant  one— dimensional  theory  as  follows: 


d2l 

df2 

fl-a’  -£) 

dfv  <i-l 
d(* 

K+l 

I 

J!1 

LI  “  d(!J 

1  -  a2  — 
12  d(2 

dv 

0  , 


(7.1) 

(7.2) 


1  /  \  d2v  n  d^v 

J  =  (o-7)  — /3  —  , 
3  d{2  d{4 


i  _  a  ,  dv  ^ d3v 
3  0  d{  d{3 


(7.3) 

(7.4) 


where  v(0  is  the  v^  of  (5.16),  (5.22)  -  (5.24),  and  the  possibility  of  a  spatially  varying 

j dQ  is  retained  in  the  equation  of  equilibrium  (7.1).  The  above  equations  are  to  be  applied 

to  V*  defined  by  (  ^  0  in  which  the  modulus  of  cohesion  0Q  is  assumed  to  be  of  different 

constant  values  p.  In  addition,  8  is  assumed  to  have  certain  singular  behaviors  at  (  = 
o  o 

0.  Specifically,  we  assume  that 

W)  =  i+(i-i30)t0(a  +  if0nc, 


(7.5) 
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where  $Q(()  and  fi(()  are,  respectively,  the  unit  step  function  and  delta  function,  such 
that 

aff0«)  =  ««>• 

•/ 

—  Q 

The  constant  /3^  is  a  new  parameter  characterizing  the  adhesion.  For  convenience,  it  is 
defined  in  terms  of  yet  another  parameter  the  modulus  of  cohesion  of  the  adhesive 

material,  by  the  expression 


Hi)  d{  =  l 


(7.6) 


a  +  a  )  a2  a1 

_  11 _ 1  2  1  11  12 _ 

a2  -  a  a  1  2  -  a  a  (a  +  a2  ] 
10  11  12  J  n  12  l  11  1 2  J 


(7.7) 


where  the  dimension  of  the  dimensional  form  of  the  factor,  which  is  introduced  for 
convenience  (c.f.  (6.9),  (6.10),  (3.34)),  is  l  and 

Dim.  ^  ,  (7  8) 

so  that  the  interpretation  of  (7.5)  is  dimensionally  correct. 

o  ± 

The  solution  satisfying  (7.1)  and  t^(0)  =  0  is 


.  *(/a  *l/a 

v^(0=A*e  »  +  B*e  12  (^0)  (7.9) 

where  the  four  constants  A*,  B*  must  now  be  determined  to  meet  the  implication  of 
(7.5).  It  follows  from  all  the  above  relations  and  additional  continuity  requirements  that 

^tuw--:u((o)=-^(0o-0o) ,  (7.io) 

v*«(0!"v;«(0)  =  '3^  '  (7'U) 

v, ^(0)  =  v,^(0)  ,  (7.12) 

r;(0)  =  -  0)  ,  (7.13) 

T*(0)  =  rj(0)  (identity),  (714) 


where  the  last  condition  is  identically  satisfied  if  (7.10)  and  (7.12)  are  met.  It  is  also  clear 
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from  (7.13)  and  (7.2)  that,  in  addition  to  (7.10)  and  (7.11),  the  4th  and  5th  derivatives  of 
v  are  also  discontinuous,  even  though  the  generalized  tractions  are  all  continuous.  The 
condition  (7.12)  conforms  with  the  ordinary  elasticity  requirement,  as  there  is  no 
concentrated  force  in  the  problem. 


The  four  constants  obtained  from  the  application  of  (7.10)  -  (7.13)  are 


(7.15) 


where  is  given  by  (6.7).  The  combination  of  (7.9)  is  a  self-equilibrating  state,  as  far 
as  the  composite  body  V  is  concerned,  and  (3.26)  gives  the  potential  energy  for  the 
system 

U*  =  ATj  ,  (7.17) 


where  Tj  is  the  interfacial  free  energy.  Since  the  interfacial  free  energy  must  be 


nonnegative,  we  have 


(7.19) 


where  the  right-hand  side  is  positive  by  (6.10). 

The  energy  of  adhesion,  E^,  is  defined  as  the  amount  of  work  necessary  to 
increase  the  separation  of  two  surfaces  from  zero  to  infinite  distance,  and  is  defined  by 
(Murr,  1975) 
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E .  ,  =  F  +  F  -  rT 
Ad  s  s  I 


(7.20) 


where  F  are  computed  by  (6.8).  The  energy  of  adhesion  must  also  be  non-negative. 


Thus, 


'  2  2 

a  &  ♦  -  a  -  a  a 

fr  <  fr =  8  +  8  +  J—ls - Li,  i.lL  . 

Ko  -  poM 


0-0 
Ho  Ho 


+  y  +  e, 


o  o 


(7.21) 


In  summary,  we  have 


fh  <  eh  <  (fir. 

Hom  ~  Ho  -  HoM 


(7.22) 


Perfect  Adhesion 


^  =  ^m'ri  =  °.E^  =  Fl  +  F;  ■ 


Ad  s  s 


(7.23) 


Adhesionless  Adhesion 


^  =  1  rr  =  +  Fs '  EAd  =  0 


(7.24) 


For  the  case  of  identical  V  ,  i.e.  0  =  0  —  0  ,  we  have 

o  o  o 


—  ph  -  oA  -A  _  « . 
-^om-^o  -  voM  ~  ^ 


(7.25) 


It  is  seen  from  (7.9),  (7.15)  and  (7.16)  that 


v,  >  e  0  for  0  =  0  =  0  ,0^=0 
i  Ho  Ho  Ho  '  Ko 


w,t  =  v 


i  ~ 


of  (6.4),  for  i  =  0o=f>, „  ,  /£  =  20, 


(7.26) 


(7.27) 


The  last  condition  indicates  that  the  traction— free  solution  is  recovered  via  the  singularity 

representation  (7.5)  without  the  specific  stipulation  that  the  generalized  tractions  r  and 
2 

t ^  be  zero.  In  this  connection,  it  is  more  convenient  to  interpret  the  singular  0Q— 
variation  of  (7.1)  as  the  body  force,  i.e. 
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Thus,  the  solution  to  two  disjointed  half  spaces,  in  the  sense  of  (7.27),  may  be  obtained 
from  the  solution  to  a  full  space  subjected  to  a  suitably  defined  singular  bodv  force.  This 
observation  opens  up  the  possibility  of  defining  a  microcrack  via  the  introduction  of 
singular  body  force  systems.  A  crack  is  termed  a  microcrack  if  its  dimensionless  length  is 
of  the  order  of  e  and,  as  a  result,  every  term  of  the  governing  equation  (4.1)  is  of  the  same 
order,  as  the  independent  variable  s  must  be  rescaled  by  e.  The  explicit 

component— function  representation  obtained  in  Section  4,  together  with  the  above 
body-force  observation,  lead  us  to  believe  that  an  explicit  solution  for  the  microcrack 
problem  is  an  achievable  goal.  This  problem  is  being  pursued  by  us. 

In  terms  of  cohesive  elasticity,  which  is  still  a  crude  continuum  model  for  solids,  a 

grain  boundary  or  an  interphase  interface  is  but  a  special  adhesive  boundary.  Thus, 

different  physical  interpretations  could  be  assigned  to  the  quantities  P^,  Tj  and  EAd 
When  a  grain  boundary  meets  a  free  surface,  a  thermal  groove  profile  for  grain 

boundary— surface  equilibrium  may  Oe  induced.  If  the  angle  sustained  by  the  groove  is 
then  the  grain  boundary  free  energy  which  may  be  identified  with  Tj  of  (7.8),  is 

Tgb  =  2FsCos  (ns/2)  .  (7.29) 

The  groove  configuration  is  a  singular  boundary  layer  phenomenon,  in  .erms  of  our 
analytical  treatement,  and  must  be  analyzed  separately.  It  is  now  also  possible  to  refine 
the  mathematical  definition  of  an  interfacial  crack  -  another  class  of  singular  problems  to 
be  solved. 

8.  Thin  Film  Energy  and  The  Concept  of  An  Interface— Phase 

A  layer  of  solid  is  called  a  plate  if  its  thickness  is  much  larger  than  the  material 
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length  scale  L.  It  is  termed  a  film  if  its  dimensionless  thickness  is  of  the  order  of  f.  The 
governing  equations  for  a  free  film  are  just  (7.1)  —  (7.4),  and  the  boundary  conditions  are 


0  t  2 

r  =  r  =  T  =  0  at  £  =  ±  h 

3  3  3  ^ 


where  2h  is  the  dimensionless  thickness.  The  solution  is 

v(0  =  ^  f-  a’  [o’  +  c!  1  Sink  —  Sia i  J- 

(*+l)Dj(h)  “l  •»  “J  a,, 


where 


+  a  a  +  a  Sinh  — —  Sinh  — ^ 
12 L  io  nj  a  a 


2  2  12 


D,(h)  =  a  a  +a  Cosh  —  Sinh - 

1  it  L  io  1 2 J  a  a 


2  2  l2  ~  h  h 

—  a  a  -+■  a  Cosh - Sinh - 

12I  io  nj  a  a 

12  11 


which  is  to  be  compared  with  of  (6.3)  -  (6.7). 

The  above  solution  is  self-equilibrating,  and  (3.27)  still  applies.  The  total  potential 
energy  for  the  system  is 

U*=AFp,  (8.21 

ff  =  e  Sr  $,[an  ~  «,] Sinh  zr Sinh  /TVh)  1  (8  3) 

11  12 

where  Fp  is  the  film  free  energy,  per  unit  area  of  the  film,  which  is  to  be  compared  with 
(6.8).  We  have 

1  im  Fp(h)  =  2F  (8.4) 

h-t  <d 

where  Fs  is  the  free-surface  free  energy  (6.8).  The  limit  simply  states  the  fact  that  the 

free  surfaces  of  a  (thick)  plate  do  not  interact.  As  h  tends  to  zero,  we  have 

1  im  Fp(h)  =  e  h/  a  -a  a  ]  ,  (8.5) 

h-*o  0  l  io  n  1 2 J 

which  is  to  be  interpreted  as  the  mathematical  approximation  for  Fp(h),  for  small  values 
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of  h.  Since  Fp  must  be  non-negative,  we  have  (c.f.  (3.34)) 


a  >  a  a 
10  11  n 


(8.6) 


It  follows  from  the  above  that  (6.10)  may  be  replaced  by 

a  -a  a  >fa  -fa  ]  fa  a  W2  ,  (8.7) 

10  11  12  l  11  12  J  L  11  12  J 

Here  we  see  that  material-parameter  properties  are  not  solely  governed  by  the  positive 
definiteness  of  the  strain  energy  density  (c.f.  the  paragraph  before  (3.34)).  It  is  anticipated 
that  many,  many  more  conditions  of  the  kind  will  eventually  come  out  of  many  different 
physical  situations.  It  can  be  easily  checked  that  Fp  is  a  monotonically  increasing 
function  of  h. 

Gluing  two  half  spaces  and  a  layer  together,  the  layer  in  between  is  termed  an 
interface— phase  if  its  dimensionless  thickness  is  of  the  order  of  e.  Let  us  assume  that  the 
two  half  spaces  are  identical,  so  that  only  half  of  the  symmetric  problem  needs  to  be 
analyzed.  The  set  of  equations  (7.1)  —  (7.4)  still  applies,  and  the  lower  half  of  the  problem 
is  defined  by 


r'«>  a+5# 

=  .  <».«> = 
v(0  [P0 


-  l 


(— h<  (  <0) 
(  (  <— h) 


(8.8) 


We  still  keep  the  simple  assumption  that  all  the  other  material  properties  are  the  same  for 
both  the  half  spaces  and  the  layer. 

The  associated  self-equilibrating  solution  is 

v.AO  =  ~  Cosh  +  a  Cosh  +  ^>/k2  (— h<£<0)  (8.9) 

11  11  12  12 


where 


A  (£+h)/a  Bn  U+h)/a 
v>((0  =  ire  n  +  ^e  12  «<-; 


h) 


u 


12 


A  =  — ASinh-^-  ,  B  =  -  B  Sinh  —  , 
o  a  ’  o  a 

11  12 


(8.10) 


(8.11) 
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A  =  A* /(Cosh.  —  +  Sinh  — )  ,  B  =  B*/(Cosh  —  +  Sinh  -^-), 


11 


u 


12 


12 


A*  = 


11 


a2  -  a2 


*— i  i  B*  _ 

^77  >  B  -~ 


12 


a  -  a 
12  11 


k — i  i 

h2 


2  H+l  l2 


(8.12) 


(8.13) 


1 2  11 

The  total  potential  energy  of  the  system  is  again  given  by  (3.27)  except  that  PQ{i) 
is  now  the  function  defined  by  (8.9).  We  have 


rn  =  ‘ 


u*  =  Arn  , 
r*=il 


(8.14) 


^ +1 J  h4[a22  -  a2J 


12 


[hCosh  h  -  a 

Sinh  — 1 

2 

a 

hCosh - a 

Sinh  — — 

a 

i  2  a 

ii 

a 

i  i  a 

l  1  2 

u. 

1 1 

lij 

Sinh  — —  +  Cosh 

a  a 

12  1  2 


Sinh  —  +  Cosh  — 
a  a 

u  li 


(8.15) 


where  Tjj  is  the  interface-phase  free  energy,  per  unit  area  of  the  interface,  which  is  to  be 
compared  with  Tj  of  (7.18)  where  the  interface— phase  was  taken  as  a  idealized  dividing 
surface.  We  have 


1  im  rTT(h)  =  e 

h-0  11 


r*-i] 

U+ij 

N 

/3a  a  (a  +  a  ) 
Al  u  12'  n  i2y 


(8.16) 


which  is  to  be  taken  as  the  mathematic  limit  of  for  small  values  of  h.  It  is  noted 

from  (8.15)  and  (8.16)  that  r j j  is  always  positive,  regardless  of  the  sign  in  front  of  in 
(8.8).  In  other  words,  Tjj  remains  the  same  for 

4>«>  “Vj  4  [(£]’- ']  <817> 

The  difference  between  the  above  and  (8.9)  lies  in  the  fact  that 


37 


f>oit)<f>0  .  for(8-9> 


(8.18) 

4(0  >  0O  .  *30<°)  =  0O  +  5  4  for  (8.17) 

The  strength  of  such  a  joint  is  dictated  by  the  minimum  value  of  PQ{()- 

The  PQ(()  variation  defined  by  (8.8)  is  arbitrary.  It  is  used  to  illustrate  the 
meaning  of  an  interface-phase.  However,  it  is  clear  that  0Q(()  may  be  required  to  yield  a 
minimum  for  Tjj.  This  proposition  is  not  pursued  further  in  this  paper.  It  is  also 
apparent  that  the  energy  associated  with  a  film  deposited  on  top  of  a  substrate  can  also  be 
meaningfully  defined. 


9.  Stretching  of  Thin  Films. 

Consider  now  a  (thin)  film  occupying  the  region 


<  eh  or 


(\<  h 


where  (  =  z  /e.  The  film  is  biaxially  stretched,  so  that 
3 

u,  =  e  z  ,u=ez  ,  u  =  rv(f) 

1  11  1  2  22  2  3 


(9.1) 


(9.2) 


where  e  and  e  are  the  constant  elasticity  strains.  For  convenience,  we  define  r  and 
11  22  J  11 

r  by 
22 


K+l  ,  3-* 
t  =  — T  e  H - r  e  , 

U  K—  1  i,  22 


K+l  3—fi 

r  =  - r  €  H - r  e 

22  K—  1  22  11 


(9.3) 


which  are  obtained  by  substituting  e  ,  e  ,  e  =0  into  the  ordinary  stress-strain 

11  22  33 

relation  with  k  being  the  first  expression  (not  the  plane— stress  version)  of  (3.18). 

In  terms  of  cohesive  elasticity,  the  nonzero  components  of  the  strains  are 

E  =e  ,E  =  f  ,  E  =  v, /  ; 

11  11  22  22  33  { 


^111  E  V,U  '  ^313  (2V‘U( 


(9.4) 
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The  nonzero  components  of  the  stresses  are 


Tv,(+7t  v,^  (a=l,2;  no  sum)  , 

(9.5) 

+  e  )  +  v,>  +  Tv,,  t ,  ; 

11  22'  K—  1  ’{  ' 

(96) 

T  =  eav,f.  , 

333  ’(( 

(9?) 

T  =  T  =  e(a  +  i  a  )v,  >  >  ; 

322  232  v  1  2  2'  u 

(9.8) 

+  7V,{  +  +  7i(«m  +  «M)1  , 

(9.9) 

'{3lT1('„  +  <!1)^+ Vafll 

•(  +  <«,  +  +  ^  v'((£l  U  • 

(9.10) 

13|7.(<U  +  ',A0+ Yafll 

,)  v,(  +  (20  +  20t  +  ^v.^l  ia0}  , 

(9.11) 

TaW  -  'Is  +  'n)Safr6  +  3  *7<  +  ‘<nS06  + 

+  3  +  ‘al’PS  +  6tfaS> 


+  3  K  +  V'f  +  ■  <912> 

The  governing  equation  for  v  is  just  (7.1),  and  the  traction— free  conditions  on  (  = 


±  h  are 


5«  f(i  -  a 1  j!)(i  _  a!  A)  Jxl  +  3qs 

* 1  L  11  d/2  12  (it!  Kr-l 


d(2  12d{2  Q 

(a-7)i!l-^  =  o  , 

df 2  d(< 


e  =  0  , 
aa  ’ 


V«  +  V^  +  ^-o. 


(9.13) 

(9.14) 

(9.15) 


The  solution  is 
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v  =  —  eaQi  +  A  Sinh  +  B  Sinh 

11  12 


(9.16) 


where 


A  =  -  f+f  “.(“I,' +“,>  Sinh  IT  [^o  +  (7.  "  Iff  ll'aal/'V1')’  (9'17) 

12 

B  “  Si  «  ,(“,+•.>  Sinh  7T  K  +  <  V -  Si  1>£aa]  /Dl<h>’  <918> 

11 

in  which  D^(h)  is  just  the  expression  given  by  (8.3).  The  solution  reduces  to  (8.2),  which 

was  obtained  for  a  film  with  a  fixed  edge,  when  e  =0. 

°  aa 

On  the  cross-section  z  =  constant,  the  stress-traction  is 

i 


t  =  t  e  ,  t  (/)  =  T  -  2T  +  3T 

11  1  11  11  311 i 3  3311 >33 


(9.19) 


which  follows  from  (3.12).  It  is  noted  that  the  stress— traction  is  no  longer  independent  of 


the  thickness  coordinate.  We  have 


=  T„ +  7(0  ■ 


Similarly,  the  stress-traction  on  the  cross-section  z^  =  constant  is 


(9.20) 


<0  -  If  T-< +  v2 Vai)v'«{  +  {20,+0,+0>)','UU(-  (9'21) 


t  =  t  e,t  (0  =  T  -  2T  +  3T 

22  2  22v  '  22  322.3  3322.33 


t  U)  =  T  +  T{()  . 
22  22 


The  average  of  r(f)  over  the  thickness  is 


r(«d( 


(9.22) 


(9.23) 


“  H  [If  *<h>  +  (27  +T,-2a &&  + 

=  ~^-faa"[,3o  +  (7,"lil7)'aa]  F(h) 


(9.24) 
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where 

p<h>  -  Sr  [sir  *  -  (21+ v2  Va2>  -  (20 +20 +0, )(°-7 )/0 

x(a  -a  )  Sinh  — —  Sinh  — —  /  hD,  (h)  . 
v  n  12  a  a  '  lv  ’ 

11  12 

The  averaged  forms  of  (9.20)  and  (9.23)  are 

(t  )  =  r  +  (r)  ,  (t  )  =  t  +  (r)  . 

11  11  22  22 

A  necessary  condition  for  a  film  of  finite  size  to  be  completely  free  is  (t  )  = 

follows  that  there  exists  an  e°  such  that 

e  =e  =  e°  ,  r  =  r  =  e°  , 

11  22  11  22  k-1 


2(  7- 
“s+T 


Trlft-[0o  +  (7,-|+f'')  2<t]  F(h)  =  0  , 


(9.25) 


(9.26) 


(t  )  =  0.  It 
22 


(9.27) 


(9.28) 


where  the  last  condition  is  merely  the  requirement  that  (t  )  =  (t  )  =  0.  The  associated 

ll'  22 

solution  is  therefore  self-equilibrating,  and  (3.27)  again  applies.  The  total  potential  energy 
of  the  system  is 

U‘  =  AF  *  ,  ^  [p0  +  ( Tj  -  -r)2«»]  Fp  ,  (9.29) 

where  Fp  is  the  film  free  energy  associated  with  a  film  with  a  fixed  edge,  (8.3).  The  new 
film  free  energy  Fp  is  associated  with  a  film  with  a  free  edge.  It  must  be  nonnegative  and 
also  less  than  Fp,  and  hence 


-4,  <<vl+T  i)2'0  <0 


(9.30) 


It  can  be  shown  that  for  positive  0Q,  which  is  presumed  throughout  this  paper,  a  material 
volume  contracts  in  the  presence  of  a  traction— free  surface.  For  example,  v,^  of  (8.2)  is 
negative.  It  follows  that  e°  <  0  which,  together  with  the  second  inequality  of  (9.30), 


implies  that 


3— K  2(  K-l)  3—K/  ,  w  n 

\  -  s+1 1  =  -isr1  \  -  s+T  (■ V"»,)  > 0 


(9.31) 


The  first  of  (9.30)  and  (9.28)  leads  to  the  conclusion  that 

Ffh)  =  2 7  +  1—2 a-a\-{ 00 +70  +  ^ )( 0-7)  /  /*]  Fp/  e/9?  <  0, 

if 

1  <  k  <  3 


(9.32) 

(9.33) 


which  is  a  condition  for  the  ordinary  theory  of  elasticity.  Since  Fp  is  positive,  we  have 

7  -  (27 +7r2a-a  )  -  (20  +20  +0^0-1)10  <  0  .  (9.34) 

A  unidirectional  tension  test,  in  the  z  -direction,  performed  on  the  free  film  is 
defined  by 

e  =  e°  +  e  ,  e  =  e°  +  e  ,  (9.35) 

11  11  22  22 


(t  +^e  +  (7_3=*  )F(h)le  =0,  (9.36) 

'  22'  «— 1  22  *-l  ii  ^2—1  wi  *+T"  ^  y  aa  ’  K  ' 


<t  )  «  2(e  -e  )  ,  (9.37) 

11  1!  22' 

where  (9.36)  has  been  used  in  defining  (t  ). 

The  vanishing  of  (t  ),  (9.36),  yields 

e!!  =  -2rG(h)  (9-38) 

where 

G(h>  =  [1-5ffe)(7l-^fl)FH/[l-!^i(7.“fel)F(h)]  (939) 

It  follows  from  (9.31),  (9.32),  (9.33)  and  (9.25)  that 


G(h)  >  1  ,  lim  G(h)  =  1  . 

(9.40) 

h-t  ® 

The  tension-test  result  finally  becomes 

(tii>=2[1+^G(h)]eu  . 

(9.41) 

If  the  film  is  actually  a  plate,  i.e.  h-t®,  the  elasticity  result  is  recovered,  i.e. 


(t^)  =  ^en  =  2(l+i/)  e  (Elasticity) 


(9.42) 
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where  2(l+i/)  is  to  be  identified  with  Young’s  modulus,  as  Dim.(t^)  =  M(ln)-  The 
cohesive  elasticity  solution  is 

(t  )  =  (  2(l+i/l  +  ^  |G(h)-l]|  e  .  (9.43) 

It  follows  from  (9.40)  that  the  apparent  modulus  for  a  film  is  greater  than  Young’s 
modulus  and,  in  fact,  size  dependent.  It  is  widely  known  that  the  apparent  strengths  of 
some  materials  are  effected  by  strain  gradients.  Using  a  couple— stress  theory,  which 
employs  only  a  portion  of  the  strain  gradients,  Mindlin  (1962)  was  able  to  show  a  reduction 
in  stress— concentration  factor  around  a  small  hole.  The  above  result  is  consistent  with  the 
general  belief  that  increasing  strain  gradients  appear  to  make  some  materials  stronger. 
Finally,  the  analysis  of  this  section  has  essentially  established  the  groundwork  for  a  new 
plane— stress  theory  suitable  for  films. 

10.  Concluding  Remarks 

The  inclusion  of  higher  order  gradients  of  displacements  in  a  continuum  theory  leads 
to  boundary-layer  phenomena  that  are  absent  in  all  grade— 1  theories  of  which  the  ordinary 
theory  of  elasticity  is  a  special  case.  This  character  may  be  exploited  to  facilitate  the 
introduction  of  nontrival  self— equilibrating  states,  so  that  surface  phenomena  of  solids  may 
be  defined  by  special  displacement  fields.  If  material  isotropy  is  presumed,  a  the.  :y  must 
include  at  least  the  first  three  gradients  of  the  displacements  before  a  material  constant  can 
be  introduced  to  capture  the  desired  surface  phenomena  (Mindlin  1965).  This  is  the 
grade— 3  theory  of  Mindlin  and  the  constant  that  gives  rise  to  all  the  surface  phenomena 
was  termed  the  modulus  of  cohesion  by  him.  It  is  for  this  reasc  that  we  have  coined  the 
name  cohesive  elasticity. 

The  most  striking  approach  in  Mindlin’s  formulatic  .  is  the  inclusion  of  the  linear 
term  e2/3E...  in  the  strain  energy-density  W  of  (3.17).  It  c  'es  rise  to  a  constant  triple 
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stress  and  leads  to  the  concise  energy  expression  U*  of  (3.27).  This  practice,  however,  is 
a  matter  of  convenience  in  the  sense  that  one  could  have  left  the  term  out  of  W  but 
introduced  the  constant  triple  stress  as  an  initial  or  residual  state.  In  fact,  one  could  have 
interpreted  the  triple  stress  as  a  special  body  force  in  the  sense  of  (7.28).  The  important 
point  is  that  however  we  interpret  it,  it  is  possible  to  introduce  a  new  constant  to  generate 
a  self-equilibrating  state. 

It  is,  of  course,  completely  legitimate  to  include  the  linear  term  in  the 

elasticity  version  of  strain  energy— density.  Such  a  term  may  also  be  interpreted  as  a  initial 
or  residual  stress.  The  stress,  however,  must  be  completely  relieved  in  the  presence  of  a 
free  boundary,  as  elasticity  equations  do  not  exhibit  boundary— layer  behavior  and,  as  a 
result,  the  exercise  leads  to  no  useful  results.  It  is  obviously  for  this  reason  that  a  grain 
with  an  eigenstrain  must  interact  with  a  nontrivial  surrounding  (Mura,  1982,  Eshelby, 
1957).  A  nontrivial  displacement  field  is  associated  with  a  free  grain  according  to  cohesive 
elasticity.  We  have  analytically  defined  the  meaning  of  grain  boundary  energy  in  Section 
7.  Unlike  ordinary  elasticity,  cohesive  elasticity  requires  the  implementation  of  comer 
(boundary)  conditions  in  case  the  boundary  is  not  a  smooth  surface.  In  this  connection,  it 
is  instructive  to  recall  the  presence  of  a  corner  force  in  a  plate  theory  (a  grade— 2  theory). 
A  notch—  or  crack-tip  and  a  corner  on  a  grain  are  among  many  meaningful  problems  to  be 
analyzed.  Even  the  meaning  of  a  concentrated  force  may  be  re— examined.  We  have 
already  found  that  the  displacement  under  a  point  load  is  finite. 

A  grade-2  theory  contains  couple  stress  and  double  stress.  It  is  mathematically 
more  tractable  than  a  grade— 3  theory  and  also  exhibits  boundary— layer  behavior. 


However,  since  it  is  not  possible  to  include  a  term  linear  in 


in  W,  it  is  not  possible 


to  introduce  a  material  constant  to  induce  a  self-equilibrating  state  for  defining,  say, 


surface  free  energy.  On  the  other  hand,  one  could  introduce  an  initial  or  residual  double 


stress,  instead  of  triple  stress,  to  induce  a  self-equilibrating  system.  The  constant  double 


stress  is  not  a  material  constant  any  longer,  but  is  still  associated  with  a  special  chunk  of  a 
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material.  This  possibility  may  turn  out  to  be  a  simpler  way  of  generating  the  desired 
surface  phenomena. 

Higher  gradient  theories  were  the  popular  topics  of  research  of  the  sixties.  Of  the 
many  known  published  results,  Mindlin’s  work  consistently  placed  the  most  emphasis  on 
relating  to  real  microscopic  phenomena  of  solids.  If  the  goal  of  today’s  micromechanics 
initiative  is  actually  aiming  at  establishing  physically  relevant  microcontinuum  mechanics 
field  theories,  then  Mindlin’s  work  are  indeed  of  a  pioneer  nature  and  should  be  carefully 
reexamined. 

Every  time  the  word  continuum  is  mentioned,  the  term  constitutive  relation  follows 

One  of  the  commonly  asked  questions  during  the  sixties  was  how  to  determine  the  so  many 

constants  involved  in  a  grade— n  theory.  It  is  our  belief  now  that  the  negative  implication 

of  the  question  was  a  result  of  our  preoccupation  with  the  term  constitutive  relation  which, 

in  one  way  or  another,  is  perceived  as  a  response  to  a  uniform  sample.  The  existence  of  a 

uniform  sample,  say,  the  tension  specimen,  is  a  consequence  of  the  assumption  that  the 

behavior  of  a  solid  is  completely  characterized  by  E-j.  While  the  concept  of  Cauchy  stress 

and  stress  vector  is  not  dependent  on  the  former  presumption,  they  are  actually  consistent 

in  that  the  components  of  the  surface  traction  vector  would  turn  out  to  be  just  the 

appropriate  components  of  the  stress  tensor  if  the  latter  were  defined  in  terms  of  a  strain 

energy— density  via  the  stationary  potential  energy  approach  of  Section  2.  It  is  clear  from 

Section  2  that  once  higher  gradients  are  included  the  exact  meaning  of  a  Cauchy  stress 

o 

becomes  vague  and  the  direct  connection  between  stress  and  stress  vector  t,  is  gone.  In 
general,  there  is  no  uniform  sample  to  speak  of  anymore.  Many  new  nontrivial  problems  of 
a  more  fundamental  nature  must  first  be  solved  before  new  test  specimens  can  be  devised. 
The  surface  free  energy,  unidirectional  test  solution,  and  point  load  solution  (not  included 
in  this  paper)  are  but  a  few  of  the  anticipated  test  configurations. 

No  continuum  field  theories  can  be  a  hundred  percent  physical.  For  example,  it 
could  be  argued  that  Young's  modulus  should  be  derived  from  more  fundamental  physical 


45 


constants,  but  the  "correct"  Young’s  modulus  is  always  measured  from  a  tension  specimen 
which  was  designed  in  accordance  with  the  mathematical  solution  that  a  =  P/A  and  e  = 
A/L.  it  is  our  oeiief  that  many,  many  more  P/A  formulas  must  and  will  be  discovered  in 
the  future.  Finally,  in  connection  with  failure,  the  all  important  Eshelby  tensor  must  be 
modified.  The  tensor  itself  could  be  straightforwardly  defined  by  casting  the  spatial 
derivatives  of  W  into  a  divergence  free  form,  but  the  associated  conservation  laws  remain 
to  be  examined.  It  is  even  possible  to  couple  the  current  grade— 3  theory  with  a  higher 
order  diffusion  equation,  the  Cahn-Hilliard  equation  (Cahn  and  Hilliard,  1957,  Gurtin, 
1989),  so  that  the  modulus  of  cohesion  will  depend  on  concentration  and  one  of  the 
diffusion  coefficients  will  depend  on  the  triple  stress.  It  is  perhaps  time  again  to  reexamine 
the  role  of  higher  gradients  in  our  pursuit  for  the  understanding  of  failure  of  solids. 
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